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The  work,  on  OMR  Contract  #N00014-84-K-0027\ has  progressed 
steadily  during  the  first  year  of  the  contract,  '^A  major  focus 
of  the  first  year  was  on  the  mixed-mode  fracture  problem. 
Theoretical  work,  on  the  computational  methodology  for  the 
calculation  of  mixed-mode  stress  intensity  factors  was 
performed.  The  convergence  properties  of  several  algorithms 
were  delineated  and  guidelines  for  the  accurate  calculation  of 
mixed-mode  stress  intensity  factors  were  established.  This 
study  also  investigated  the  effects  and  influence  of  loading 
holes  on  the  calculations,  (loading  holes  are  usually  ignored  in 
theoretical  studies,  if  they  are  modeled,  they  are  usually 
modeled  poorly).  While  previous  studies  have  investigated  the 
accurate  modeling  of  loading  holes,  their  influence  on  fracture 
specimens  had  not  been  sufficiently  explored.  In  addition  to 
these  fundamental  studies,  a  proposed  mixed-mode  specimen  was 
carefully  analyzed.  This  specimen  (originally  designed  by 
Professor  H,  Richard)  has  been  proposed  for  studies  near  the 
range  of  pure  Mode  II  loading  (so-called  "pure  shear").  The 
calculations  demonstrated  the  utility  of  the  specimen  for 
measuring  critical  shear  conditions  for  pure  Mode  II  failure. 
This  study  was  performed  in  conjunction  with  an  exchange 
student  to  The  George  Washington  University  (Mr.  Peter 

Bauerle),  His  final  report  carefully  documents  the  work, 
performed  and  is  included  as  Appendix  A  of  this  report. 

Concurrent  with  the  theoretical  work.,  independent 
experimental  verification  of  this  fracture  specimen  was  per- 


formed.  Several  refinements  were  made  to  the  original  design. 
This  design  (which  involves  a  very  specialized  grip 
configuration)  was  constructed  and  tested.  PMMA  specimens  were 
tested  to  determine  the  failure  envelope  under  the  full  range 
of  mixed-mode  loading.  In  addition,  uncracked  specimens  were 
tested  to  determine  their  resistance  to  shear  loading.  The 
results  demonstrate  the  ability  of  this  loading  grip  system  to 
accurately  induce  the  full  range  (from  pure  Mode  I  to  pure  Mode 
II)  of  planer  mixed-mode  loading  conditions.  This  study  was 
also  performed  in  conjunction  with  an  exchange  student  to  The 
George  Washington  University  (Mr.  Roland  Gerstner).  His  final 
report  completely  covers  the  work,  performed  and  is  included  as 
Appendix  B  of  this  report. 

Most  of  the  previous  recent  fracture  research  sponsored  by 
ONR  focused  on  studying  the  nature  of  plastic  deformation  in 
the  vicinity  of  three-dimensional  crack,  fronts.  These  studies 
were  completed  during  this  contract  year.  Experimental  surface 
contractions  were  measured  and  compared  with  finite  element 
predictions.  For  the  worst  cases,  the  results  agreed  within 
4%.  The  plasticity  model  employed  and  the  finite  element  grid 
convergence  were,  therefore  demonstrated.  In  addition,  several 
global  response  parameters  were  compared  with  an  experiment. 

It  was  shown  that  global  parameter  agreement  does  not  guarantee 
accurate  local  modeling.  The  results  of  this  study  are 
included  in  a  paper  entitled  "Prediction  of  Plasticity 
Characteristics  for  Three-Dimensional  Fracture  Specimens  - 


Comparison  with  Experiment"  which  has  been  accepted  for 
publication  in  Engineering  Fracture  Mechanics.  This  paper  is 
included  as  Appendix  C  of  this  report - 

The  three-dimensional  plasticity  studies  have  also 
demonstrated  many  of  the  local  response  characteristics  near  a 
crack,  front.  Having  established  the  accuracy  of  the  modeling, 
the  response  of  the  interior  of  the  specimen  can  now  be 
accurately  predicted.  The  influence  of  hardening,  specimen 
thickness,  local  relaxation,  crack  growth,  etc.,  was  accurately 
predicted  and  investigated  for  straight  crack.  Mode  I 
specimens.  While  no  conclusions  regarding  appropriate  failure 
criteria  have  been  drawn,  it  was  demonstrated  that  the  failure 
was  definitely  governed  by  a  local  parameter  and  it  should  be 
related  to  the  local  deformation  response.  Previous  studies 
have  indicated  the  inadequacy  of  a  global  failure  parameter, 
however,  the  local  deformation  characteristics  have  only  now 
been  delineated  and  verified.  The  results  were  presented  as  an 
invited  paper  at  the  American  Society  of  Metals  meeting  held  in 
Salt  Lake  City,  Utah  (December,  1985).  This  paper  is  included 
as  Appendix  D  of  this  report.  This  paper  will  subsequently  be 
published  in  the  proceedings  of  this  conference  by  the  American 
Society  of  Metals. 

The  major  focus  of  this  contract  has  been  to  develop  the 
methodology  and  to  begin  a  systematic  study  into  the  fracture 
behavior  of  specimens  subjected  to  high  temperatures.  To  this 
end,  a  test  system  for  the  high  temperature  and  mechanical 
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loading  of  specimens  was  procured,  installed  and  tested.  After 
verifying  the  instrumentation,  preliminary  studies  were 
performed  to  delineate  the  constitutive  behavior  of  Inconel  718 
when  loaded  at  elevated  temperatures  (temperatures  above 
SOO^C).  Uniaxial  tests  were  performed  and  the  data  was  "fit" 
with  several  constitutive  models.  Upon  the  completion  of 
uniaxial  testing,  the  experimental  work,  focused  on  establishing 
the  instrumentation  necessary  to  study  crack,  growth  phenomena. 
Attention  is  focusing  on  accurate  measurement  of  crack,  mouth 
opening  behavior  for  modeling  purposes.  In  addition, 
preliminary  tests  are  being  run  to  determine  the  requisite 
loading  to  produce  different  crack,  growth  rates  of  importance 
for  application.  The  experimental  work,  to  date  is  summarized 
in  Appendix  E  of  this  report.  This  Appendix  also  covers  the 
constitutive  modeling  being  done. 

Concurrent  with  the  experimental  and  constitutive 
modeling,  work,  has  focused  on  developing  and  testing  a  finite 
element  code  to  perform  fracture  analyses.  The  preliminary 
code  using  a  simple  constitutive  model  has  been  completed  and 
tested.  Work,  is  underway  to  build  in  crack,  growth 
capabilities,  large  strain  modeling  and  alternative 
constitutive  models  with  mixed  hardening  properties.  The 
formulations  being  employed  are  discussed  in  Appendix  F  of  this 
report . 

In  addition  to  the  work,  discussed  above,  two  other  publi¬ 
cations  were  produced  which  summarized  and  concluded  previous 


work,  performed  under  ONR  support.  These  papers:  "Finite 
Element  Methodology  for  Elastic-Plastic  Fracture  Problems  in 
3-Dimensions"  and  "Effect  of  Specimen  Thickness  on  Crack.  Front 
Plasticity  Characteristics  in  Three-Dimensions"  are  included  as 
Appendices  G  and  H  of  this  report. 

As  is  evident  from  this  report  and  the  publications 
included,  much  has  been  accomplished  in  the  first  year  of  this 
contract.  The  next  year  will  be  spent  establishing  the  basic 
experimental  and  numerical  approaches  for  accurate  modeling  and 
investigation  of  creep  fracture  processes.  While  our  efforts 
are  producing  useful  and  informative  results,  we  have  only 
scratched  the  surface  of  this  problem. 
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ABSTRACT 


Stress  intensity  factors  are  material  independent 
parameters  which  measure  the  stress  and  strain  distribution  at 
the  tip  of  a  crack,  in  a  loaded  solid  body.  These  stress 
intensity  factors  contain  the  geometry  of  the  crack..  They 
enable  to  predict  if  a  prescribed  load  can  be  applied  to  a 
damaged  material  without  losing  the  required  safety. 

In  this  research  all  calculations  were  done  with  the  finite 
element  program  APES  which  uses  12-node  quadrilateral 
isoparametric  elements.  In  order  to  obtain  the  singularities 
at  the  crack,  tip,  the  following  two  crack,  tip  elements  were 
utilized; 

enriched  elements 

collapsed  cubic  isoparametric  elements 

(1,  §  4  point  elements) 

9  9 

Since  the  stiffness  matrix  of  the  enriched  elements  already 
contains  the  singularity  functions,  Kj  and  Kjj  can  be 
calculated  together  with  the  displacements  by  the  computer 
program.  Using  the  collapsed  cubic  isoparametric  elements,  the 
stress  intensity  factors  Kj  and  Kjj  have  to  be  determinated 
with  the  aid  of  crack,  border  displacement  functions  after  the 
run  of  the  program. 

These  calculations  were  done  for  a  special  specimen 
developed  by  H,  A.  Richard  for  mixed  mode  problems.  Also,  the 
influence  of  the  holes  where  the  forces  are  applied  was  studied. 
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1 .  INTRODUCTION 


During  the  last  20  years  many  investigations  have  been  done 
to  determine  the  stress  intensity  factor  (Mode  I)  where 
the  forces  are  applied  perpendicular  to  the  crack..  However, 
most  cracks  lead  to  mixed  mode  problems  where  also  the  stress 
intensity  factor  Kjj  (Mode  II,  see  Figure  1)  has  to  be  taken 
into  consideration.  For  this  case  it  is  important  to  develop 
methods  which  allow  to  determine  both,  Kj  and  more 

accurately.  Only  dependable  stress  intensity  factors  allow  an 
engineer  to  decide  whether  a  damaged  part  of  a  machine  or  a 
building  can  be  used  without  losing  the  required  safety. 

In  this  paper,  the  finite  element  program  APES  will  be 
utilized  to  determine  the  two-dimensional  stress  intensity 
factors  K|  and  Kp.  The  12-node  quadrilateral 
isoparametric  elements  of  this  program  already  lead  to  a  high 
accuracy  with  a  relative  small  number  of  elements,  compared  to 
lower  order  elements  of  other  finite  element  programs. 

On  the  following  pages,  only  plane  strain  will  be  taken 
into  consideration.  Plane  strain  means  that  along  the  z-axis 
(third  dimension)  no  strains  are  allowed  according  to  the 
Poisson’s  ratio.  Therefore,  in  addition  to  the  stresses  in  x- 
and  y-direction  also  stresses  in  z-direction  have  to  come  into 
being.  In  reality,  this  state  of  stress  can  just  be  realized 
in  the  middle  of  the  crack  tip  of  a  Compact  Tension  Shear 
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specimen  if  the  specimen  is  thick,  enough.  Plane  strain 
represents  the  crucial  load  case  of  a  CTS  specimen  which  leads 
to  cleavage  fracture  without  any  prior  plastic  deformation. 


2.  INTRODUCTION  TO  THE  12-NODE  QUADRILATERAL  ISOPARAMETRIC 

ELEMENT 

The  APES  program,  with  which  this  research  was  done, 
utilizes  the  12-node  cubic  isoparametric  element  (shown  in 
Figure  2).  One  12-node  element,  given  in  the  x,y  coordinate 
system,  can  be  transformed  into  the  natural  n  coordinate 
system  where  the  element  is  mapped  to  a  square  and  the 
coordinates  are  dimensionless  [1]. 

The  displacement  components  can  be  written  in  the  form 

12 

u=  z  N. (C,n)  u.  (1) 

i=l  ^  ^ 

12 

v=  I  N. CC.n)  V.  (2) 

i=l  ^  ^ 

u,v  are  the  components  of  the  displacements  of  a  point  with  the 
coordinates  C,  Ht  whereas  u^^,  are  the  displacement 
components  of  the  node  i.  Nj^(5,  n)  represent  the  shape 
functions  or  interpolation  functions  at  the  node  i  which  are 
given  by 

NiCC,n)  «  -j^Cl  +  CCi)Cl  +  10  +  9(5"  +  n")]- 

[-  10  +  9(5?  +  n?)] 

C3) 

+  TS^Cl  +  CCi)(l  +  9nni)(l  -  n")(l  -  n?) 

*  +  955i)Cl  -  C")(l  -  C?) 


where  ri^  are  the  coordinates  of  the  point  i.  It  has  to  be 
mentioned  as  the  important  feature  of  the  cubic  isoparametric 
elements,  that  the  same  shape  functions  are  utilized  to 
describe  the  geometry  of  the  12-node  elements  in  the  x,  y 
coordinate  system. 


12 

E 

i=l 

N^CC.n)  X. 

(4) 

12 

E 

i=l 

N^CC.n)  X. 

(5) 

where  x^,  represent  the  nodes  of  the  element  in  the  x,  y 
coordinate  system.  Evaluating  the  shape  functions  at  each  node 
i.  N(  n  )  becomes: 


where  [B]  is  given  by 


[e]  -  [B][u] 


[e] :  strain  field 

[u] :  displacement  field 


[B]  = 


[  •  •  •  »B  j 

[B^] 


.] 


3N^ 


3N. 


3N 


w~ 


1 
3x 


[D]  is  the  material  stiffness  matrix  for  plane  stress 


[D] 


l-v" 


V 

1 


0 

0 

Cl-v) 


for  plane  strain  [2] ; 


*  C1T0)X1T?v) 


(1-v) 

V 

V 

n 


V 

(1-v) 

V 


n 


Cl^v)Cl-2v) 
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[J]  represents  the  Jacobian  matrix  (transformation  matrix) 


If  det  I J I  is  zero,  stress  and  strains  become  singular. 

The  advantages  of  the  12-node  quadrilateral  isoparametric 
elements  are  [3]: 

the  edges  of  the  elements  can  take  the  shape  of  a 
cubic  function.  Therefore,  the  12-node  elements  can 
be  easily  adapted  to  many  practical  problems, 
strains  and  stresses  vary  cubically  over  the  element, 
which  means  that  only  a  few  elements  are  sufficient  to 
simulate  the  required  stress  and  strain  distribution 


in  a  solid. 
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3.  CRACK-TIP  ELEMENTS 


In  order  to  simulate  the  nodal  crack,  tip  singularity  in  a 
proper  way,  the  finite  elements  around  the  crack,  tip  have  to 
reproduce  the  strain  singularity.  The  finite  element  program 
APES  already  contains  enriched  elements  which  tak.e  the  strain 
singularity  into  account  with  the  aid  of  specific  singular 
assumptions. 


3.1  Enriched  Elements 


3.1.1  Introduction  to  the  Enriched  Elements 
The  enriched  elements  have  the  same  shape  lik.e  the  ordinary 
cubic  isoparametric  elements  and  they  can  be  adapted  to  the 
required  geometry  of  the  specimen  in  the  same  way.  In  order  to 
simulate  the  strain  singularity  at  the  crack,  tip  special 
singular  functions  are  utilized. 

S.  E.  Benzley  describes  the  displacement  assumptions  for 
enriched  elements  with  4  nodes  [4].  For  12-node  elements  these 
displacement  functions  can  be  written  as  follows: 


“l  ■  *k  “ik  *■  fjt  ’Jiiic) 


*  "llWa  ^k  '?2ik” 


^  -•••  .•»  Jr  - - -  ■ 


vrhere  =  nodal  displacements 

f  a  shape  functions 

^li’^2i  “  special  singular  assumptions 

^lik.’^2ik.  “  value  of  Qi£.Q2i  evaluated 
at  node  k 

R  a  correction  function. 

The  correction  function  R  provides  for  a  disappearance  of 
the  singular  assumptions  at  the  borders  where  both,  enriched 
and  12-node  cubic  isoparametric  elements,  have  to  fit 
together.  However,  these  compatibility  conditions  are  just 
satisfied  directly  at  the  nodes.  Between  the  nodes,  the 
displacements  of  these  two  different  elements  are  not 
coincident . 

Strains  are  obtained  by  the  derivation  of  the  displacements 
and  they  are  not  continuous  at  the  border  of  enriched  and 
ordinary  12-node  elements.  Therefore,  strains  and  stresses  are 
not  correct  close  to  the  crack  tip  and  they  can  not  be  utilized 
for  the  determination  of  the  stress  intensity  factors.  The 
stiffness  matrix  of  the  enriched  elements  becomes  [4]: 
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.V. 

•  *  •  ' 


rll 

I 

,22 


,12 


where  u  «  element  nodal  displacements 
regular 

stiffness  matrix  from  enriched  terms 
coupled  stiffness  matrix  from  regular  and 
enriched  terms 

F  »  singular  load  vector  which  becomes  a  null  vector 
if  no  enriched  element  is  on  a  loaded  boundary. 

The  equation  system  (11)  points  out  that  all  singular 
assumptions  of  the  enriched  elements  are  contained  in  the 
stiffness  matrix  and  that  the  stress  intensity  factors  Kj  and 
Kjj  can  be  evaluated  directly  by  the  program  with  the 
displacements . 

Furthermore,  it  is  advantageous  that  the  size  of  the 
enriched  elements  can  be  quite  large  without  losing 
considerable  accuracy.  Therefore,  complete  problems  can  be 
solved  quite  easily  within  a  short  data  preparation  time. 
Usually,  enriched  elements  are  working  within  an  accuracy  of 
about  5\. 

However,  it  has  to  be  noted  that  smaller  enriched  elements 

need  not  necessarily  lead  to  a  higher  accuracy  which  is  the 

case  for  i  §  4.  point  elements. 

9  9 

3.1.2  Application  of  the  Enriched  Elements 
The  determination  of  the  stress  intensity  factors  Kj  and 
Kjj  was  done  using  a  special  specimen  developed  by 
H.  A.  Richard  for  mixed  mode  problems  [5],  shown  in  Figure  3. 
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First,  a  modified  compact  tension  shear  specimen  was 

studied,  where  the  influence  of  the  holes  was  neglected  (Figure 

4).  Equally  distributed  stresses  were  applied. 

With  the  aid  of  the  sinus/cosinus  relationship  each  applied 

force  F  under  the  load  angle  a  can  be  decomposed  to  the  forces 

F,  and  F^  where 
z  s 

F2  leads  to  pure  Mode  I 
Fs  leads  to  pure  Mode  II 

It  should  be  mentioned  that,  in  order  to  obtain  pure  shear 
stress  on  the  elongated  crack,  line  for  case  Mode  II,  also  a 
moment  M  with 

M  X  Fs  •  c 

has  to  be  applied  in  opposite  direction  of  F^.  Applied 
tensile  and  shear  stresses: 

for  Mode  I :  o, 

z 

for  Mode  II: 

The  shear  stress  distribution  was  chosen  in  the  form  of  a 
parabola  where 

W  =  1-5  ^o  (Figure  5) 


wt 

wt  ^bmax  ~  v  ~ 
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Figure  6  and  7  point  out  the  finite  element  grid  of  the  compact 
tension  shear  specimen  with  the  ideal  force  application. 
Altogether,  20  finite  elements  are  used,  whereas  4  enriched 
elements  are  arranged  around  *:he  crack,  tip. 

In  Figure  6,  the  boundary  conditions  for  pure  shear  are 
used:  all  nodes  on  the  x-axis  are  constrained  in  x-  and 

y-direction.  It  should  be  added  that  for  pure  shear,  according 
to  the  Moor's  strain  circle  (Figure  8),  e  vanishes  in  shear 
direction,  whereas  y  reaches  its  maximum.  However,  because  of 
the  bending  moment  M  there  exist  tensile  and  compressive 
stresses  along  the  x-axis  (Figure  4)  which  should  allow  some 
displacements  of  the  nodes  at  the  boundaries  in  x-direction. 
This  is  not  possible.  Since  the  maximal  tensile  and 
compressive  stresses  exist  at  the  corners  and  since  the 
boundaries  are  far  away  from  the  crack,  line,  these  effects  can 
be  neglected. 

In  Figure  7,  the  boundary  conditions  for  pure  tension  are 
shown.  Only  one  node  on  the  x-axis  has  to  be  constrained  in  x- 
and  y-direction,  whereas  all  other  nodes  on  the  boundaries  can 
move  in  x-direction.  This  takes  the  Poisson's  effect  into 
consideration. 

Chart  1,  Figure  9  and  10  point  out  the  dimensionless  stress 

intensity  factors  Kj  and  Kjj  for  load  application  angles 

changing  from  0  to  90  degrees,  a  was  0.6. 

w 
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3.2  i  §  4  Point  Elements 
9  9 

3.2.1  Introduction  to  the  i  S  4  Point  Elements 

9  9 

In  order  to  achieve  the  stress  and  strain  singularity  at 
the  crack,  tip,  collapsed  cubic  isoparametric  12-node  elements 
can  be  utilized  (Figure  11)[71.  These  crack.-tip  elements  have 
to  be  distorted  as  follows: 

all  nodes  lying  on  the  side  of  the  crack,  tip  have  to 
overlap  (nodes  1,10,11,12). 

the  nodes  lying  on  the  sides  close  to  the  crack,  tip 

are  placed  on  the  points  with  the  distances  i  I, 

9 

4  i  and  i  away  from  the  tip  (therefore  the  name  i  &  4 
9  9  9 

point  elements). 

The  mathematical  proof  of  the  singularity  for  these  elements  is 
given  in  [7]. 

The  advantages  if  these  i  §  4_  point  elements  are: 

9  9 

one  does  not  need  special  crack  tip  elements. 

Ordinary  12-node  elements  can  be  utilized. 

smaller  i  S  4  point  elements  will  lead  to  a  higher 
9  9 

accuracy,  which  is  not  necessarily  the  case  for 


enriched  elements. 


However,  these  have  to  be  noted  as  disadvantages: 


because  of  the  higher  number  of  elements,  the  1.  S  4 

9  9 

point  element  method  is  more  time  consuming  than  the 
enriched  element  method. 

the  stress  intensity  factors  Ki  and  Kji  have  to  be 
calculated  in  post-process. 

The  stress  intensity  factors  Kj  and  Kjj  can  be  obtained 
by  using  the  following  displacement  equations  at  the  crack  tip 
(Figure  12)[8]: 


-  IT  /It  I  tjCk-l)  +  sin^  j] 

/It  7  [7^^^^^"  7^ 

7  ^7^^'^^  7^ 
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3.2.2  Optimum  Size  of  the  1.  8  4  Point  Elements 

9  9 

In  order  to  find  out  the  optimum  size  of  the  184  point 

9  9 

elements,  the  problem  shown  in  Figure  13  was  studied  with  the 

appropriate  boundary  conditions. 

Figure  14  points  out  the  finite  element  grid  for  shear  with 

which  the  best  results  could  be  obtained.  Altogether,  30 

finite  elements  were  used,  6  crack,  tip  elements  were  arranged 

around  the  crack,  tip.  The  length  i  of  the  184  point  elements 

9  9 

was  1.23%  of  the  length  of  the  crack  a. 


16 

For  shear,  the  stress  intensity  factors  Kj  and  Kjj  were 

evaluated  along  the  6  -  180*  line,  since  this  is  the  only  line 

where  the  tensile  stresses  vanish.  The  equations  (17)  and  (18) 

were  utilized.  The  results  can  be  seen  in  Figure  15. 

Because  the  nodes  of  the  crack,  tip  are  constrained  in  both 

directions  the  stress  intensity  factors  calculated  very  close 

to  the  tip  have  to  be  neglected  (about  first  two  elements). 

The  stress  intensity  factors  Kj  and  Kjj  directly  at  the 

crack  tip  were  obtained  graphically. 

Figure  16  shows  the  results  of  a  12  element  mesh  using  4 

crack-tip  elements  (shear).  The  length  of  the  184  point 

9  9 

elements  was  1%  of  the  crack  length.  This  grid  turns  out  to  be 
too  rigid. 

Figure  17  points  out  the  results  of  a  36  element  grid  using 

6184.  point  elements  around  the  crack  tip.  It  turns  out  that 
9  9 

just  increasing  the  elements  without  using  more  crack  tip 
elements  need  not  necessarily  lead  to  a  higher  accuracy  (Length 
of  the  crack-tip  elements;  0.667%  of  the  length  of  the  crack.) 

For  pure  tension,  the  best  results  could  be  obtained  along 
the  6  a  90*  line  where  the  equations  (15)  and  (16)  were  used 
(Figure  18). 

The  dimensionless  Ki  and  Kji  factors  are  (a  =  0.6): 

w 


The  dimensionless  Ki  and  Kji  factors  are  (a  -  0.6): 

w 


«  [“] 

0 

90 

Kj 

It 

3.99  (4.12) 

0  (-  0.01) 

(-)  0.11  ((-)0.05) 

1.57  (1.63) 

Fcr  a  equal  to  zero  degree,  the  dimensionless  Kj  factor 
becomes  3.99,  whereas  Kjj  turns  out  to  be  -  0.11.  Usually 
K'jj  should  vanish  for  this  loading  case.  However,  this 
discrepancy  should  not  be  overevaluated,  since  for  Mode  I,  all 
nodes  at  the  crack  tip  are  constrained  in  both  directions. 

This  restriction  can  lead  to  different  results. 

For  a  equal  to  90  degrees,  the  dimensionless  Kjj  factor 
turns  out  to  be  1.57,  whereas  Kj  vanishes. 

The  values  in  parentheses  show  the  results  of  the  enriched 
element  method. 

3.2.3  Application  of  the  i  S  4  Point  Elements  to  the 

9  9 

GTS  Specimen 

Now,  the  1  §  4  point  elements  are  applied  to  the  compact 
9  9 

tension  shear  specimen.  The  stress  intensity  factors  of  the 

GTS  specimen  are  determined  with  the  aid  of  the  finite  element 

grid  of  Figure  19.  Twelve  154  point  elements  are  arranged 

9  9 


around  the  crack  tip.  This  grid  is  utilizing  the  boundary 


conditions  already  shown  in  Figure  6  and  7.  Ideal  distributed 
forces  are  applied  to  the  specimen. 

There  are  several  reasons  why  the  whole  finite  element  grid 
of  the  CTS  spcimen  (Figure  19)  has  to  be  used  and  why  it  can 
not  be  replaced  by  the  grid  shown  in  Figure  14. 

First,  in  Figure  14  all  nodes  at  the  crack,  tip  are 
constrained  in  x-  and  y-direction.  This  leads  to  restrictions 
which  need  not  represent  the  reality.  All  nodes  at  the  crack, 
tip  should  be  able  to  move  in  all  directions  in  order  to  adapt 
to  the  required  situation.  This  restriction  can  be  the  reason 
for  the  inaccurate  stress  intensity  factors  close  to  the  crack, 
tip  in  Figures  15-18. 

Secondly,  the  boundary  conditions  of  Figure  14  (shear)  do 
not  simulate  the  boundary  conditions  for  pure  shear.  According 
to  equation  (14),  all  nodes  along  the  boundaries  should  be  able 
to  move  in  y-direction  (0  =  0®,  Kj  *  0,  Kjj  »  0)  except  for 
the  nodes  at  the  crack,  tip.  In  Figure  14,  this  would  lead  to 
instabililty  and  could  not  be  realized. 

The  finite  element  grid  of  Figure  19  does  not  use  any 
boundary  conditions  at  and  around  the  crack,  tip. 

In  Figure  20  and  21,  it  can  be  seen  how  the  stress 
intensity  factors  are  evaluated.  For  pure  tension  (Figure 
20a, b),  the  stress  intensity  factors  were  determined  along 
the  9  »  90®  line  (Figure  12,  equations  15  and  16),  whereas  for 
pure  shear  (Figure  21a, b),  the  0  »  180®  line  was  utilized  along 
the  upper  crack,  surface  (equations  17  and  18),  The  calculations 


along  the  upper  and  lower  crack,  surface  lead  to  identical 

results.  In  Figure  20  and  21,  the  k.  factors  of  all  nodes  along 

the  9  »  90®  or  9  a  180®  line  could  be  used  in  order  to 

determine  the  stress  intensity  factors  at  the  crack,  tip.  For 

all  calculations,  the  nodal  displacements  relative  to  the  crack. 

tip  were  taken  into  consideration-  Since  all  nodes  of  the 

crack  tip  are  allowed  to  move  separately,  only  that  crack  tip 

node  was  chosen  as  the  reference  point  which  was  laying  on 

the  6  =  90°  or  9  =  180®  line  of  the  respective  1.  8  4.  point 

9  9 

element . 

In  Chart  3,  Figure  22  and  23,  the  dimensionless  stress 

intensity  factors  are  shown  for  a  »  0.6.  For  pure  tension,  the 

w 

dimensionless  Kj  factor  turns  out  to  be  3.98,  whereas  Kjj 

vanishes.  For  pure  shear,  the  dimensionless  stress  intensity 

factor  Kjj  becomes  1.56,  whereas  Kj  is  zero. 

The  results  of  the  184  point  elements  are  compared  with 

9  9 

those  of  the  enriched  elements  and  of  Richard  [9].  Richard 

used  a  superelement  in  order  to  simulate  the  crack  tip 

singularity.  In  general,  the  stress  intensity  factor  of  the 

184  point  elements  and  the  results  of  Richard  are  laying  very 
9  9 

close  together.  The  differences  are  equal  to  or  smaller  than 
1%  and  can  be  neglected.  The  enriched  elements  are  a  little 
bit  more  inaccurate. 

Figure  24  and  25  point  out  the  location  of  the  crack-tip 


20 

nodes  after  applied  shearing  and  tensile  stresses.  For  shear, 
the  nodes  are  more  distributed  In  direction  of  the  x-axls, 
whereas  for  Mode  I.  the  crack  tip  nodes  are  more  arranged  along 
the  y-axls.  Unfortunately,  It  Is  not  possible  to  find  out  In 
which  direction  the  crack  will  propagate.  The  nodes  of  the 
crack- tip  elements  are  randomly  distributed. 

Figure  26  and  27  show  how  the  crack  will  behave  for  the 
loading  cases  tension  and  shear.  For  pure  tension  (Figure  26), 
the  nodes  on  the  upper  crack  surface  undergo  larger 
displacements  than  the  nodes  of  the  lower  crack  surface. 
However,  in  the  vicinity  of  the  crack  tip,  the  displacement  of 
the  corresponding  nodes  on  both  surfaces  are  identical. 

Figure  27  shows  the  behavior  of  both  crack  surfaces  for  pure 
shear.  Along  the  whole  crack  line,  both  surfaces  are  sliding 
against  each  other  as  expected.  Away  from  the  crack  tip,  both 
crack  surfaces  are  penetrating  slightly.  This  penetration  is 
very  small  compared  with  the  displacements  in  x-direction  or 
the  displacements  in  Figure  26.  It  lies  within  the  accuracy  of 
the  program  and  should  not  be  overevaluated.  It  can  be  said 
that  for  pure  shear,  no  penetration  occurs  along  the  crack 


surface. 


LOADED  HOLES  WITHIN  A  PLATE 


4.1  Modelling  of  the  Holes 

Up  to  now,  just  the  ideal  force  application  to  the  compact 
tension  shear  specimen  has  been  taken  into  consideration.  In 
reality,  the  forces  are  applied  with  the  aid  of  bolts  which 
could  lead  to  a  different  stress  and  strain  distribution  around 
the  crack  tip. 

In  most  of  the  previous  examinations  just  a  nodal  force  was 
assumed  in  the  middle  of  the  holes.  Also  Richard  [10]  was 
working  with  single  nodal  forces  when  he  was  studying  the 
stress  intensity  factors  of  his  CTS  specimen. 

However,  there  are  required  so  huge  forces  in  order  to  be 
able  to  tear  apart  a  CTS  specimen  that  a  plastic  deformation 
comes  into  being.  This  plastic  deformation  leads  to  a  force 
application  over  a  broader  cross  section  at  the  holes.  No  one 
knows  how  broad  this  plastic  zone  will  be. 

In  this  research  another  approach  was  done.  It  was  assumed 
that  the  bolts  fit  exactly  into  the  holes  of  the  CTS  specimen 
and  that  the  forces  can  be  applied  over  the  hole  cross  section 
of  the  holes.  All  calculations  were  done  for  the  pure  elastic 
case. 

Figure  28  and  29  point  out  the  finite  element  grids  of  a 


arbitrary  plate  and  a  bolt.  The  boundary  conditions  of  the 
bolt  are  imitating  the  hole  of  the  plate:  the  points  1-7  of 
the  bolt  can  just  move  along  the  tangent  of  the  hole.  For 


both,  hole  and  plate,  the  symmetry  could  be  used. 

When  the  bolt  is  loaded,  a  nodal  force  F  is  applied  in  the 
center  of  the  bolt  in  negative  y-direction  which  leads  to  an 
elastic  deformation  and  to  reaction  forces  in  the  nodes  1-7  in 
y-direction.  These  forces  (Fj^  ~  ^7^  applied  at  the 
nodes  1*  -  7*  of  the  plate.  The  displacements  of  these  nodes 
(1*  -  7*)  can  be  taken  as  prescribed  displacements  for  the 
nodes  1-7  of  the  bolt.  After  9  iterations,  the  changes  of  the 
forces  in  the  nodes  1-7  of  the  bolt  were  less  than  1%.  The 
obtained  forces  and  displacements  are  given  in  Chart  2. 

It  should  be  mentioned  that 
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IF.*  1.3873  -  1.3873  F 
i*l  ^ 

The  reason  is  that  the  finite  element  method  is  based  on  the 
principle  of  virtual  displacements  which  leads  to 

SWg:  work,  done  by  external  forces 

6We  =  6Wi  C19) 

•SW^:  work,  done  by  internal  forces 


a  ■ 


stress  matrix 


e  -  strain  matrix 
b  -  body  force  matrix 
p  -  exterior  surface  force  matrix 
u  >  displacement  matrix 
A  -  surface  area 
V  »  volume 

Because  of  the  separation  of  bolt  and  plate,  for  the  plate, 
just  the  nodal  forces  have  to  be  taken  into 

consideration.  Therefore,  the  resultive  force  which  is  applied 
to  the  hole  of  a  plate  is 

13 

Fres  ”  1.3873  F  -  F.  (21) 

Figure  30  shows  the  deformation  of  the  bolt  and  the  plate 
under  the  influence  of  the  applied  force  F. 

Figure  31  and  32  point  out  the  dimensionless  stress 
distribution  around  the  bolt  and  the  hole  of  the  plate. 

4 . 2  Application  of  the  Holes  to  the  Simplified  Compact 
Tension  Shear  Specimen 

First,  the  holes  are  applied  to  the  simplified  CTS  specimen 

of  Figure  12  which  leads  to  Figure  33  (a  »  0.6). 

w 

In  Figure  34  it  can  be  seen  how  the  stress  intensity 
factors  Kj  and  Kjj  are  evaluated  for  shear.  It  turns  out 


that  the  difference  between  the  stress  intensity  factors  using 
a  CTS  specimen  with  and  without  holes  is  less  than  0.2%.  For 


both  cases  Kj  is  pretty  exactly  zero  (shear). 
The  dimensionless  factors  are: 


including  the  holes: 


without  the  holes: 


0.191 


1.573 


1.576 


which  lies  within  the  accuracy  of  the  procedure  and  can  be 
neglected.  Since  for  pure  shear,  the  shear  force  is  applied  to 
a  single  hole  and  the  other  two  holes  are  utilized  to  apply  the 
bending  moment,  this  case  can  be  considered  as  the  extreme 
loading  case.  For  all  other  loading  cases  the  difference  of 
the  stress  intensity  factors  using  specimen  with  and  without 
holes  will  not  be  more  than  0.2%.  Therefore,  for  the  following 
problems  the  influence  of  the  holes  can  be  neglected  and  only 
CTS  specimen  with  the  application  of  equally  distributed  forces 
will  be  used. 


4.3  Application  of  the  Holes  to  a  CTS  Specimen 


Figure  35  shows  the  CTS  specimen  with  enriched  elements 
including  the  holes.  In  this  grid  the  ideal  boundary 
conditions  could  be  utilized  (statically  stable): 

at  the  upper  three  holes  the  required  forces  are 
applied 

the  nodes  of  the  lower  holes  are  fixed  in  the  shown 
way:  they  can  move  in  one  direction.  In  the 
direction  perpendicular  to  it  prescribed  displacements 
are  used  in  order  to  achieve  the  same  force 
distribution  which  is  applied  to  the  upper  three  holes. 

With  this  grid,  the  boundary  conditions  of  Figure  6  should  be 
tested  which  are  not  100%  correct.  Unfortunately,  this  mesh 
does  not  work,  because  the  allowed  bandwidth  of  the  program  was 
exceeded. 


5.  EXAMINATION  OF  DIFFERENT  CRACK  LENGTHS 


Up  to  now,  only  one  crack,  length  (a  =  0,6)  has  been  studied 

w 

In  this  chapter,  also  two  other  crack,  lengths  (a  =  0.5  and 

w 

a  =  0.7)  will  be  examined  with  the  aid  of  a  finite  element  grid 
w 

similar  to  that  one  of  Figure  19.  Equally  distributed  forces 
are  applied  to  the  specimen. 

In  Figure  36a-39b,  it  can  be  seen  how  the  stress  intensity 

factors  are  evaluated  for  pure  tension  and  pure  shear.  All 

results  are  put  together  in  Chart  4,  Figure  40,  41  and  42. 

For  pure  shear,  the  dimensionless  Kjj  factors  are 

l,37(a  *  0.5),  1.56(a  =  0.6)  and  1.74(a  =  0.7),  Except  for 
w  w  w 

a  =  0.5,  all  Ki  factors  vanish.  For  a  *  0.5,  Ki  is 
w  w 

slightly  negative  {-  0.02),  which  would  lead  to  a  penetration 
of  the  crack,  surfaces.  For  this  loading  case. 


is  equal  to  1.4%  which  lies  within  the  accuracy  of 


the  program  and  should  not  be  overevaluated. 

All  obtained  stress  intensity  factors  are  compared  with  the 

results  of  Richard  [9],  The  differences  between  the  results  of 

the  i  S  4  point  element  method  and  the  results  of  Richard  are 
9  9 


small,  except  for  the  crack,  length  a  =  0.5  (a  =  90°)  where 


Richard  obtains  a  Kt  factor  of  +  0.08.  The  1.  §  £  point 

9  9 

elements  lead  to  1  Kj  ]  =  +  0.02  which  is  smaller. 

For  pure  tension,  the  dimensionless  Kj  factors  become 

2.81(a  =  0.5),  3.98(a  =  0.6)  and  6.19(a  =  0.7).  All  stress 
w  w  w 

intensity  factors  Ku  are  zero. 


6 .  SUMMARY 


In  this  research,  the  stress  intensity  factors  Kj  and 

Kji  were  evaluated  for  a  special  compact  tension  shear 

specimen  developed  by  Richard.  Two  different  cracR-tip 

elements  were  utilized:  enriched  elements  and  1^4  point 

9  9 

elements.  The  finite  element  grid  using  enriched  elements 

leads  to  quite  accurate  results  with  a  relative  low  number  of 

elements.  In  order  to  obtain  more  accurate  results,  i  &  4 

9  9 

point  elements  have  to  be  used.  For  this  case,  the  size  of  the 

finite  elements  arranged  around  the  crack,  tip  has  to  be  very 

small.  The  optimum  results  could  be  achieved  for  i  §  4  point 

9  9 

elements  with  an  element  length  of  about  1.23%  of  the  crack 
length. 

Furthermore,  the  holes  were  examined  where  the  forces  are 
applied  to  the  compact  tension  shear  specimen.  Two  different 
finite  element  grids  were  compared;  one  mesh  with  a  force 
concentration  around  the  holes  and  another  finite  element  grid 
with  equally  distributed  forces  at  the  whole  cross  section. 

Both  finite  element  grids  led  to  the  same  results  which  means 
that  the  CTS  specimen  is  so  well  designed  that  the  holes  do  not 
have  an  influence  on  the  stress  intensity  factors. 

Altogether,  the  stress  intensity  factors  were  evaluated  for 
three  different  crack  lengths.  For  pure  shear,  the 
dimensionless  stress  intensity  factors  K,,  become 


1.37 (a  »  0.5),  1.56 (a  =  0.6)  and  1.74 (a  0.7).  Kj  is  pretty 
w  w  w 

exactly  zero  for  all  three  crack,  lengths. 

For  pure  tension,  the  dimensionless  stress  intensity 

factors  Ki  turn  out  to  be  2.81(a  »  0.5),  3.98(a  »  0.6)  and 

w  w 

6.19(a  >  0.7),  whereas  the  Kxj  factors  vanish. 


^  J  J  > 


30 


REFERENCES 


[  1]  S.  L.  Pu,  M.  A.  Hussain  and  W.  E.  Lorensen,  "The  Collapsed 
Cubic  Isoparametric  Element  as  a  Singular  Element 
for  Crack  Problems,"  in  International  Journal  for 
Numerical  Methods  in  Engineering.  Vol.  12..  PP. 
1728-1730,  1978. 

[  2]  K.  Kussmaul,  "Festigkeitslehre  I,"  Vorlesunesmanuskr lot . 
MPA  Stuttgart,  p.  169,  1982. 

[  3]  L.  N.  Gifford,  "APES  -  Second  Generation  Two-Dimensional 
Fracture  Mechanics  and  Stress  Analysis  by  Finite 
Elements,"  DTNSRDC  Report  #4799,  p.  12,  December 
1975.  - 

[4]  S.  E.  Benzley,  "Representation  of  Singularities  with 

Isoparametric  Finite  Elements,"  in  International 
Journal  for  Numerical  Methods  in  Engineering.  Vol. 
8,  pp.  538--542,  1974. 

[5]  H.  A.  Richard,  "Einleitung  von  iiberlagerter  Normal-  und 

Scherbeanspruchung  in  Bruchmechanikproben," 
University  of  Kaiserslautern,  1981. 

[6]  M.  M.  Filonenko-Boroditsch,  "Festigkeitslehre  Band  I," 

VEB  Verlag  Technik  Berlin,  p.  191,  1954. 

[  7]  M.  M.  Pu,  M.  A.  Hussain  and  W.  E.  Lorensen,  "The 
Collapsed  Cubic  Isoparametric  Element  as  a 
Singular  Element  for  Crack  Problems,"  in 
International  Journal  For  Numerical  Methods  in 
Engineering.  Vol.  12,  pp.  1730-1734,  1978. 

[  8]  J.  Eftis,  N.  Subramonian  and  H.  Liebowitz,  "Crack  Border 
Stress  and  Displacement  Equations  Revisited,"  in 
Journal  Fracture  Mechanics.  Vol.  9,  p.  189,  1977. 

[  9]  H.  A.  Richard,  Habilitationsschrift ,  University  of 
Kaiserslautern,  p.  213,  Chart  15,  1984. 

[10]  H.  A.  Richard,  Habilitationsschrift,  University  of 
Kaiserslautern,  p.  168,  1984. 


MEASUREMENTS 
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Stress 

Young’s 

Poisson 

Stress 


1  in.  s 
1  lb.  . 


INCH 

kilo  pound 

<?,T)  pound/inch2 

Modulus  (E) 
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Intensity  Factor  (k) 
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2.54cm  »  2.54  x  10“2  m 
453. 6g  =  0.4536  k.g 
3  JL 


1  psi  *  6.8948  X  10" 


-  334.5 

-  278.3 

-  216.9 

-  131.7 


-  375.8 

-  334.5 

-  278.3 

-  216.9 

-  131.7 


18.6 

52.1 

60.2 


<|B> 


Mode  I 
Opening  Mode 


Mode  II 
Sliding  Mode 


L  - 

y  =  0:Tw«=* 


3Q  ,  la 

bh  ■  2  F 


b;  thickness 
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Figure  19b:  Detail  A 
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Figure  22:  Dimensionless  K,  Factors  for  -  =  0 
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Figure  23;  Dimensionless  Kj,  Factors  for  =  0.6 


Figure  24:  Location  of  the  Crack-Tip 
Applied  Shearing  Stresses 


Magnification  Factor  for  Displacements:  i  x  10 

0*1  psi 
a  =  3  in. 
w  =  5  in. 
t  =  2.5  in. 

Figure  26;  Crack  Near  Crack  Tip  After 
Applied  Tensile  Force 


Figure  30:  Deformation  of  the  Bolt  and  of  the 
Plate  After  Applied  Forces 
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Figure  37b:  K, 
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Abstract 


In  a  literature  survey  it  was  found  that  the  CTS 
specimen  proposed  by  Richard  is  the  most  suitable  specimen 
for  mixed  mode  experiments,  which  were  conducted  with  PM/1A. 
As  it  was  found  in  tests  with  CT  specimens  that  the  critical 
stress  intensity  IJi^of  PMMA  depends  on  the  stress  intensity 
rate  K,  the  testing  was  done  at  a  constant  K.  The  fracture 
cr i ter  ion 


gives  a  reasonable  fit  to  the  experimental  data.  Some 
problems  in  mode  II  testing  with  the  CTS  specimen  were  found 
that  need  further  consideration.  In  additional  tests  the 


modulus  of  elasticity  and  the  Poisson  ratio  were  obtained. 


« 

4 

i-- 

•> 

'y. 

^^•v 

S 

s 

s 

> 

\ 

CONTENTS 

W-  ' 

"j*  V* 

i&S 

1 

1  INTRODUCTION 

PAGE  C 

1 

• 

II  NOTATION 

2 

III  LITERATURE  REVIEW 

mF 

1 

1.  MIXED  MODE  SPECIMENS 

3  -S':- 

•! 

2.  FRACTURE  CRITERIA 

9 

i 

3.  EARLIER  EXPERIMENTAL  WORK 

U.  SUMMARY 

13  vS; 

5.  PLAN  OF  INVESTIGATION 

20 

IV  EXPERIMENTAL  PROCEDURE 

1 

1.  MATERIAL 

21 

2.  TESTING  PROCEDURE 

21 

i 

V  RESULT  AND  DISCUSSION 

■ 

1.  MECHANICAL  PROPERTIES 

26 

2.  FRACTURE  TOUGHNESS  AT  VARYING  STRESS 
RATES 

INTENSITY 

27 

i 

3.  FRACTURE  TOUGHNESS  UNDER  MIXED  MODE 

LOADUNG 

29  h 

VI  CONCLUSIONS  ^ 

35 

*.■ 

■ 

y. 

i 

• 

k  •*. 

# 

V 

i'  1 

I  Introduction 


While  mode  I  fracture  mechanics  are  well  established  and 
commonly  used,  mixed  mode  fracture  mechanics  did  not  reach 
more  than  a  prelimary  stage  yet.  There  is  no  standardized 
testing  procedure.  There  are  many  different  failure  criteria 
but  none  of  them  can  be  judged  as  true  or  false  because  of 
the  lack  of  reliable  data. 

On  the  other  hand,  mixed  mode  fracture  is  of  great 
practical  intrerest  as  almost  no  real  structure  is  so  simple 
that  a  crack  is  subjected  to  pure  opening  mode.  Mixed  mode 
loading  occurs  at  curved  and  inclined  cracks,  branched 
cracks,  welded  or  glued  joints,  and  in  many  other  cases. 

This  paper  starts  with  a  literature  review  of  mixed  mode 
fracture  mechanics  with  three  points  of  interest:  mixed  mode 
specimens,  fracture  criteria,  and  earlier  experimental  work. 
As  the  values  for  mechanical  properties  of  plexiglas  found 
in  literature  differed  considerably,  testing  of  these 
properties  was  done.  A  series  of  tests  with  CT  specimens  {1} 
was  conducted  to  find  out  the  dependancy  of  the  critical 
stress  intensity  K,^  on  varying  stress  intensity  rates. 
Finally,  mixed  mode  experiments  with  the  compact  tension 
shear  (CTS)  specimen  proposed  by  Richard  {2-8}  were 


performed . 


The  following  literature  review  consists  of  the  three 
parts  mixed  mode  specimens,  fracture  criteria,  and  earlier 
experimental  work.  From  this  review  a  plan  of  investigation 
is  concluded. 


1.  Mixed  Mode  Specimens 

To  date,  there  is  no  standadized  specimen  for  mixed  mode 
experiments.  During  the  literature  review,  10  different 
specimens  of  general  use  were  found  which  will  be  discussed 
briefly.  These  specimens  are  (see  figure  1) 

51  Panel  with  an  angled  center  crack  under  tension  (ACCT)  , 
sometimes  with  notch  (ACNT)  . 

52  Panel  with  a  single  angled  edge  crack  under  tension 
(SAECT) 

53  Panel  with  a  curvilinear  center  crack  under  tension 
(CCCT) 

SA  Panel  with  a  center  crack  under  biaxial  ’’ension  (CCBT) 

55  Single  edge  cracked  specimen  under  j-point  bending 
(SEC3B) 

56  Single  edge  cracked  specimen  under  l4-point  bending 
(SECAB) 

57  Tube  under  torsion  with  an  oblique  crack 

58  Braz i 1 i an  disk  (60) 

59  Specimen  proposed  by  Banks-Sills 

SIO  Compact  tension  shear  specimen  (CTS) 
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Beside  these,  spherical  shells  under  pressure  {10,11}  and 
X-,  H-,  and  T-shaped  specimens  were  used  {12}.  The 

abbreviations  in  parentheses  are  used  in  the  following. 

As  with  the  most  often  used  ACCT  specimen  no  pure  mode  II 
condition  can  be  acieved,  the  following  specimens  were 
employed  for  mode  II  testing 

511  Compact  shear  specimen  (CS) 

512  Inplane  shear  specimen  (IPS) 

These  specimens  will  be  discussed  briefly  as  well. 

The  panel  with  an  angled  crack  under  uniaxial  tenstion 
(figure  la)  is  the  most  often  used  specimen  {13~23}t  but 
there  are  some  variations:  Erdogan  and  Sih  {13}  used  a 

large  plate  (a/b  =  0.2)  subjected  to  "uniform  tension  at 
infitnity"  and  gave  for  that  the  stress- i ntens i ty  factors 

Kj  =  cfV^sin  p,  Kjj  r  <J\/rrasinpcosp  (3) 

which  are  used  by  other  resesarchers  as  well. 

Shah  {18}  and  Wilson  {23}  used  a  narrow  panel  (a/b  ••  0.5) 
with  only  one  loading  hole  on  each  end.  The  stress  intsiisity 
factors  for  this  configuration  are  not  explicitly  given,  but 
it  is  obvious  that  they  are  different  from  those  of  the 
large  panel.  Wu  et  al.  {34,35)  used  a  large  panel  (a/b  « 
0.17)  with  an  angled  elliptic  notch  loaded  by  three  pins  on 
each  side.  Many  researcher  do  not  state  at  all  which  kind  of 
ACCT  specimen  they  used. 

The  obvious  problem  is  that  the  results  cannot  be 
compared  simply  by  the  crack  angle  as  the  mixed  mode  ratio  n 
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is  different.  Equation  (2)  is  strictly  valid  only  for  an 
infinite  panel  loaded  at  infinity. 

The  biggest  limitation  of  this  specimen  is  that  no  pure 
mode  II  condition  can  be  achieved.  Therefore  researcher 
could  not  perform  experiments  over  the  complete  mixed  mode 
range  or  they  had  to  use  another  specimen  for  mode  II 
experiments.  The  loads  increase  considerably  for  decreasing 
angle  and  in  some  cases  the  specimens  broke  in  the  loading 
holes  for  small  angles  {21}. 

Another  problem  is  the  fatigue  precracking  as  the 
application  of  the  load  in  testing  direction  would  result  in 
curved  cracks  which  are  of  no  use.  Several  researchers 
prevented  the  problem  by  using  oniy  thin  machined  slits 
{21,25.33}  or  by  introducing  the  crack  with  a  razor  blade 
{15}  which  is  possible  in  the  case  of  PMMA.  If  fatiguing  is 
necessary  as  in  the  case  of  steel  or 
aluminum  alloys,  the  following  procedure  has  to  be  done 
{12}:  A  central  not  angled  notch  is  machined  in  a  larger 
plate.  Then  this  plate  is  precracked  in  the  usual  way.  The 
final  specimen  is  cut  out  of  this  lager  plate  afterwards,  as 
shown  in  figure  2.  This  procedure  wastes  a  lot  of  time  and 
mater i al . 

The  plate  with  an  angled  edge  crack  subjected  to  uniaxial 
tension  (SAECT)  (figure  lb)  is  used  less  often  {12,36.37} • 
The  problems  in  production  and  application  of  this  specimen 
are  similar  to  those  just  stated  for  the  ACCT  specimen. 

Hussain  et  al.  {38}  suggested  a  plate  with  a  curvilinear 


central  crack  subjected  to  uniaxial  tension  (figure  1c)  for 
mixed  mode  and  mode  II  testing  and  performed  some  mode  II 
experiments.  The  specimen  did  not  get  much  attention  as  it 
is  quite  difficult  to  produce  a  curvilinear  crack.  The  mode 
II  conditon  is  at  77*0“  (not  at  79-6*as  stated).  A  pure 
mode  I  condition  cannot  be  achieved  as  K,  gets  zero  at  'n  ■ 
137*5  and  there  is  K,  <  0. 

The  cracked  plate  under  biaxial  tension  (figure  Id)  is 
used  by  several  workers  {25.39“^i}*  There  are  some 
variations:  Some  used  a  straight  crack  and  just  changed  the 
ratio  d',/d'y  (39, AO},  others  rotated  the  crack  {25} »  and 
Radon  et  a1.  {41}  used  curved  cracks.  A  testing  machine 
suitable  to  produce  biaxial  loading  is  necessary.  Therefore 
the  effort  in  machinery  is  higher  than  for  those  specimens 
which  can  be  tested  in  uniaxial  tension. 

The  tube  under  torsion  with  an  oblique  crack  (figure  le) 
got  some  attention  {11,13}.  but  it  is  not  very  easy  to 
handle.  Also,  not  every  material  is  available  as  tube 
mater i a1 . 

The  '  hr or  four — point  bend  specimens  have  some 
advantages:  They  are  easy  to  produce  and  fatigue.  They  can 
be  tested  from  pure  mode  I  to  pure  mode  11,  although  a  pure 
mode  II  condition  may  only  be  achieved  with  a  notch,  as 
pointed  out  by  Richard  {42}.  The  bend  specimens  got  some 
attention  in  experimental  work  {37,43»44}. 

The  Brazilian  disk  specimen  which  also  can  cover  the 
complete  range  from  mode  I  to  mode  II  was  used  by  several 
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researchers  [k5,^^) .  As  the  crack  is  subjected  to 
compression  in  this  specimen,  friction  between  the  crack 
flanks  may  occur  {42}.  The  fracture  mechanism  for  a  crack 
under  compression  may  be  different  from  that  under  tension 
{47}.  Hence  the  results  between  tension  specimens  and  the 
BD  specimens  may  not  be  comparable: 

Banks-Sills  et  al.  (50)  recently  proposed  a  new  specimen 
for  mode  II  dominant  mixed  mode  experiments.  This  specimen 
is  not  suitable  for  mode  II  experiments.  No  experiments  with 
this  specimen  have  been  reported  yet. 

The  CTS  specimen  proposed  by  Richard  {2-8}  is  similar  to 
the  Banks-Sills  specimen.  But  instead  of  gluing  the  specimen 
to  the  grips,  Richard  uses  a  statically  determinant 
connection  by  means  of  pins.  This  is  easier  to  handle  in 
practical  application.  The  CTS  specimen  can  be  used  from 
pure  mode  I  to  pure  mode  II.  Richard  et  al.  {42}  claim  that 
the  CTS  specimen  is  superior  to  the  Banks-Sills  specimen  and 
all  other  specimens  in  mode  II  testing.  A  problem  with  the 
CTS  specimen  is  that  the  influence  of  geometrical  deviations 
which  are  inevitable  in  machining  the  specimens  is  not  known 
yet.  No  other  experiments  than  those  performed  by  Richard 
himself  are  reported  in  literature  so  that  the  basis  for 
direct  comparison  is  tiny. 

The  K|,g  values  which  are  necessary  to  confirm  the 
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different  mixed  mode  fracture  criteria  are  evaluated  with 
two  different  mode  II  specimens,  the  compact  shear  (CS) 
specimen  {49}  and  the  inplane  shear  (IPS)  specimen.  Although 


Chrisholm  and  Jones  {^9}  claim  to  achieve  virtually  pure 
mode  II  condition  at  the  crack  tip  of  the  CS  specimen,  later 
researchers  got  different  results:  Riddle  (50,51}  found  that 
K|  is  much  larger  than  anticipated  and  that  already  small 
variations  in  the  boundary  conditions  change  the  ratio  K|,/K, 
extremely.  And  Richard  et  al.  {Ii2}  came  to  the  conclusion 
that  "this  specimen  is  unsuitable  for  the  evaluation  of 
The  discussion  about  this  specimen  is  not  settled  yet.  A  new 
FEW  calculation  by  Elenz  (52)  shows  that  the  state  of  stress 
at  the  crack  tip  is  essentially  mode  II.  The  CS  specimen 
was  used  by  {21}  in  their  experimental  work. 

The  inplane  shear  specimen  did  not  get  that  much 
theoretical  attention.  It  was  used  by  {13,21}  for  mode  II 
testing,  but  the  introduction  of  the  forces  imposes  some 
uncertainties. 

In  general,  it  is  very  difficult  to  achieve  pure  mode  II 
condition.  Paris  and  Sih  (53)  found  that  only  internal  or 

very  deep  external  cracks  show  significant  mode  II 

displacements.  Several  researchers  (9,^3, 54}  found  in 
numerical  analyses  that  the  plastic  zone  i-ncre.?pos  with 
increasing  K||.  Beside  some  other  problems  this  imposes  that 
the  thickness  requirement  for  mode  II  specimens  may  be 
greater  than  that  for  mode  I  specimens  (54) .  A  further 

discussion  of  the  problems  and  some  other  mode  II  specimens 

are  found  in  {42}. 

The  question  remains  what  the  influence  is  on  the  value 
of  K„  if  the  state  of  stress  is  not  pure  mode  II  but 
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established  criterion  is  Sih's  strain  energy  density 
function  {l6,17f59}*  The  predictions  of  this  criterion 
depend  on  the  Poisson  ratio.  While  the  difference  to  the 
other  two  criteria  is  within  10%  for  a  Poisson  ratio  around 
0.3f  the  predictions  are  very  different  for  extreme  values 
of  V.  Therefore,  experiments  with  materials  with  an  extreme 
Poisson  ratio  are  best  to  prove  the  validity  of  this 
cr i ter  ion. 

Many  criteria  {13. 17t55»59.66)  base  on  the  assumption 
that  fracture  is  a  local  phenomen  and  that  hence  the  stress 
distribution  in  the  vicinity  of  the  crack  tip  can  be 
described  just  by  the  singular  terms  of  the  William'^  stress 
funtion  (67)  •  For  a  definition  of  the  stresses  around  the 
crack  tip  see  figure  3. 

Inspired  by  a  paper  of  Cotterel 1  {68},  Williams  and  Ewing 
{15}  pointed  out  that  the  stress  (S’,  parallel  to  the  crack  tip 
(see  figure  3)  may  not  be  neglegible.  Hence  they  added  the 
first  nonsingular  term  of  the  William's  stress  function  {67} 
in  their  analysis.  Later  Eftis,  Subromanian  and  Liebowitz 
{69.70}  confirmed  this  appiosch.  Seversl  author:-  {15.71-73} 
based  their  fracture  criteria  on  this  enlarged  function.  The 
problem  is  that  all  these  criteria  cannot  be  described  in 
terms  of  K,  and  K,  only,  but  need  an  additional  parameter. 
This  parameter  must  be  determined  experimentally  which 
usually  is  done  in  a  way  that  the  resulting  function  gives 
the  best  fit  to  the  experimental  data.  As  the  calculation  of 
possible  without  knowing  this  parameter,  these 
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criteria  are  not  listed  in  table  2. 


Host  criteria  are  for  brittle  fracture.  But  there  are 
also  four  criteria  that  take  plastic  deformations  into 
account:  the  criterion  by  Shih  {9},  the  crack  separation 
energy  rate  criterion  {k3,73},  the  T-criterion  {31-33,71*), 
and  the  J-integral  {63*65}  •  It  is  especially  interesting 
that  the  T-c iter  ion  predicts  for  pure  mode  II  an  angle  of 
fracture  initiation  of  about  -97  decree,  while  the  J-intgral 
predicts  0  degree  which  are  extreme  values. 

A  second  group  of  fracture  criteria  are  the  empirical 
criteria.  They  are  evaluated  as  best  fit  to  the 
experimental  data.  Normally  they  require  the  knowlegde  of 
K^.  The  most  general  form  of  them  is 


Awaji  and  Sato  {1*6}  found  u  -  v  -  1.6  for  graphite  and 
plaster  and  u  ■  v  -2  for  marble.  Shah  {18}  got  u  ■  v  ■  1 
for  steel  which  was  confirmed  by  Chiu  and  Liu  {39}  for 
aluminum  and  was  found  by  Leicester  {75}  for  wood.  Wu  {76} 
got  u  ■  1  and  v  -  2  for  wood  and  fiberglass  -  reinforced 
plastics  which  later  was  found  by  Richard  {1*,8}  for  PMHA. 
Lee  and  Advani  {77)  confirmed  these  exponents  in  an  analysis 
based  on  the  Griffith  theory  {78}. 

The  mixed  mode  problem  is  dealt  with  in  official 
regulations  too  {79*80} .  The  ASHE  Boiler  and  Vessel  Code 
{79}  suggests  to  project  the  inclined  crack  in  the  direction 
of  the  maximum  stress  and  subject  then  the  projected  crack 
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to  mode  1  loading.  The  fracture  criterion  is 

s  <J\/ic(asinp)  (5) 

After  some  algebraic  transformations  (see  appendix  1),  this 
equation  becomes 


''ic  • 

~2 

(6) 


This  criterion  is  dicussed  in  several  papers  {8,57>8i}>  As 
can  be  seen  in  figures  4  to  6  in  {57).  the  ASME  criterion 
makes  in  some  cases  prediction  of  higher  fracture  loads  than 
other  criteria,  especially  In  the  case  of  pure  shear.  But 
this  does  not  mean  that  the  ASHE  criterion  is  unsafe,  as 
this  is  just  a  relative  difference  and  the  other  criteria 
are  not  convincingly  verified  by  experiments  either. 

This  short  review  of  the  mixed  mode  fracture  criteria  is 
far  from  being  exhaustive.  Already  in  1977,  Gilles  {82} 
rep,',  ted  of  35  different  mixed  mode  fracture  criteria.  This 
paper  also  is  mainly  restricted  to  homogeneous  materials 
subjected  to  tension.  For  orthotropic  materials  (wood, 
composite  materials)  or  for  cracks  under  compression,  other 
criteria  can  be  found  in  literature.  For  a  further 


discussion  of  fracture  criteria  and  also  comparison  between 
them,  the  reader  is  referred  to  {8,19,57,58,81,83-85}. 


3.  Earlier  Experimental  Work 


There  is  no  standardized  testing  procedure  for  mixed  mode 
testing.  Most  researchers  took  ASTM  E399  (1)  as  a  guideline. 
But  as  already  mentioned,  the  thickness  requirement  stated 
there  may  not  ensure  linear  elastic  conditions  as  the 
plastic  zone  increases  in  mixed  mode  loading.  Pook  {14} 
introduced  the  criterion  that  the  net  section  stress  at 
fracture  initiation  should  be  less  than  80%  of  the  yield 
strength  of  the  material  to  ensure  brittle  fracture. 

Many  researchers  already  conducted  mixed  mode 
experiments.  Most  research  was  done  with  PMMA,  but  steel 
and  aluminum  alloys  got  considerable  attention  too. 

The  problems  and  results  with  PMMA  will  be  discussed 
later,  now  the  other  materials  will  be  considered.  Tables  3 
to  -5  give  the  results  found  in  literature  for  steel 
{18,19,25.33.43}.  aluminum  alloys  {1-2,14,19,20,23,28,33,86}, 
and  several  other  materials  {4,27,29,31,33,46}.  As  can  be 
seen  in  column  3  of  each  table,  most  researcher  used  the 
ACCT  specimen.  Column  5  gives  the  range  of  the  experiments 
by  the  mixed  mode  index  n.  Many  researchers  did  not  test  the 
complete  range  from  mode  1  to  mode  II.  Another  important 
information  is,  how  many  specimens  were  tested,  which  is 
given  in  column  6.  It  is  felt  that  that  3  or  4  specimens 
are  not  enough  to  make  any  conclusion  about  the  tested 
material.  Columns  7  to  9  finally  give  K|j,K,^,  and  as 
far  as  they  are  stated  in  the  ci ted' papers.  A  comparison  of 
the  experimentally  found  K„^/K,jWith  those  predicted  by  the 
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different  fracture  criteria  (table  2)  reveals  that  only  the 
J-integral  {63-65}  and  the  criterion  by  Irwin  {61,62} 
predict  K,[^  “  K,^  and  are  therefore  close  to  the  experimental 
results.  Nost  researchers  who  perfomed  these  experiments 
therefore  used  empirical  fracture  criteria. 

The  angle  of  fracture  initiation  predicted  by  the  various 
criteria  differ  considerably.  Hence  experimental 
verification  is  necessary.  An  interesting  experiment  to 
check  the  validity  of  Sih's  strain  energy  density  criterion 
{16,17,59}  was  performed  by  Finnie  and  Weiss  {1*9}.  They  used 
cross-rolled  beryllium  (V  =  0.0)  and  found  that  the  measured 
fracture  angle  is  much  higher  than  predicted  by  the  strain 
energy  density  criterion  but  is  close  to  the  prediction  of 
the  Erdogan  and  Sih  criterion  {13}.  Awa j  i  and  Sato  {146} 
measured  the  fracture  angle  in  mode  II  for  graphite  SA  Ci) 
(V  ■  0.07)  and  graphite  7^77  =  0.20)  and  found  i  t  to  be 
around  67  degrees  in  both  cases.  All  these  experiments  give 
evidence  that  the  angle  of  fracture  Initiation  is 
essentially  independent  of  the  Poisson  ratio. 

Although  the  J-integral  and  the  T-criterion  predict 
totally  different  fracture  angles  for  mode  II,  both  are 
verified  experimentally.  Riddle  {50}  performed  experiments 
with  aluminum  and  found  the  angle  to  be  around  0  degree, 
while  Theocaris  {31}  used  Polycarbonate  of  Bisphenol  A 
(PCBA)  and  found  "angles  absolutely  greater  than  90*". 

Plexiglas  got  the  largest  attention  in  mixed  mode 
fracture  research  (I* ,  1  1  ,  1 3 . 15.2 1 , 25-27 , 33 . 31* . 37 . **0, It  1 , l*l»l  . 
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The  results  are  given  in  table  6.  In  addition  to  the 
information  included  in  the  other  tables,  this  table  states 
the  way  in  which  the  crack  was  produced  (column  k) ,  the 
displacement  rate  used  for  testing  (column  5),  and  whether 
slow  crack  growth  was  reported  in  the  paper  (column  6).  Most 
researchers  found  Kji^K^  to  be  around  0.9  which  is  close  to 
the  values  predicted  by  several  fracture  criteria  (table  2). 
This  is  not  that  surprising  as  most  of  the  fracture  criteria 
were  constructed  to  explain  the  behavior  of  PMMA. 

As  it  was  planned  to  perform  the  new  experiments  with 
PMMA  too,  this  material  got  further  attention  in  the 
1 i terature  review. 

Many  mixed  mode  exper iments-  requi re  the  knowledge  of  the 
Poisson  ratio.  In  addition,  it  is  useful  1  to  know  the 
modulus  of  elasiticity  and  the  strength  of  the  material.  The 
values  found  in  literature  {^.7.11.33.37.87-9^}  are  given  in 
table  7.  The  range  for  the  modulus  of  elasticity  is  from 
2.0  to  5-0  GPa  and  for  the  Poisson  ratio  values  between  O.30 
and  0.43  were  stated.  Often  It  was  not  very  clear  whether 
the  stated  value  was  obtained  experimentally  or  whether  it 
was  chosen  as  the  best  fit  to  the  respective  theory. 

The  modulus  of  elasticity  is  dependent  on  the  straining 
rate  (95).  but  most  sources  do  not  state  this  rate.  Williams 
{95}  gives  the  equation 

E  =  [0.655  +  2.74(e)°-^°9]GPa  (7) 


where  e  is  the  straining  rate  with  theunit  1/s.  Williams 


applied  this  formula  for  e  ■  5x10'*  1/s  to  10*  1/s.  He 
therefore  got  values  for  E  from  1.6  to  5-2  GPa. 

The  fracture  toughness  testing  of  PMMA  involves  some 
problems.  In  the  following,  those  found  in  the  literature 
are  stated  and  discussed: 

1.  Most  testing  machines  are  designed  for  metals.  Therefore, 
the  lowest  load  range  is  around  5000  N.  Even  with  a  high 
accurancy  of  0.1?!,  the  error  because  of  the  small  fracture 
loads  for  PMMA  may  be  high.  This  problem  is  reported  by 
Phadke  et  al.  {21}.  They  used  a  load  range  of  500  kg  (I49OO 
N)  with  a  least  count  of  0.5  kg  (4.9  N)  .  For  the  smallest 
reported  fracture  load  of  7.5  kg  (73.6  N) ,  the  error  is 

6.7?:. 

It  is  also  not  clear  whether  ASTM  E399  {1}  is  applicable 
to  PMMA  or  not,  as  this  norm  is  designed  for  metallic 
materials.  Despite  this  uncertainty,  that  standard  was 
applied  in  many  papers,  as  no  other  standard  is  known. 

2.  Many  researchers  consider  PMMA  as  ideal  isotropic  and 
brittle,  while  other  scientists  proved  that  socalled  crazes 
establish  at  the  crack  tio  {89,96-99}.  The  crack  tip  is 
supposed  to  have  atomar  dimensions  and  the  fracture  process 
is  a  molecular  process.  As  PMMA  has  large  molecular  chains, 
the  orientation  of  these  chains  influences  the  state  of 
stress  in  a  way  that  it  is  inhomogeneous  in  molecular 
dimensions.  Crazes  only  develop  in  an  inhomogeneous  state  of 
stress  which  shows  that  this  materia]  is  not  that  ideal 
isotropic  in  a  microscopic  range.  Furthermore,  crazes  are 


kind  of  plastic  deformations  which  may  cause  some  error  in 
applying  LEFM  theories. 

It  was  suggestd  to  anneal  the  the  specimens  to  reduce  the 
crazes  {100}  which  was  done  by  (3^) •  In  addition,  Hai  {100} 
suggested  to  minimize  the  occurence  of  crazes  "by  a  coating 
of  silicone  oil  onto  the  specimen  surfaces".  The  influence 
of  this  coating  is  not  readily  understood  and  no  other 
reference  to  it  was  found  in  literature. 

3.  The  fracture  toughness  of  PMMA  depends  on  the 
temperature,  the  relative  humidity  of  the  air  or  the  water 
content  of  the  material  respectively,  and  the  presence  of 
chemicals.  Therefore,  some  researchers  stored  the  specimens 
under  specific  conditions  for.  1500  hours  {S**} .  or  did  the 
testing  in  climatized  rooms  {100}.  But  most  researchers  did 
not  have  the  facilities  and  did  tne  testing  under  whatever 
conditions  there  were  in  their  laboratories. 

I*.  As  was  discussed  in  several  papers  {95>  lOI-lOA}  ,  the 
fracture  toughness  of  plexiglas  depends  on  the  crack 
propagation  speed.  A  group  of  researchers  {95il02}  state 
that  there  exists  a  "unique  relationship  between  critical 
stress  intensity  factor  and  crack  speed  which  is  independent 
of  the  test  piece  geometry  used"  {95}.  while  Cotterell  {101} 
claims  that  "although  the  fracture  toughness  does  increase 
with  velocity  of  fracture,  ...  it  is  not  a  function  of 
velocity".  The  results  of  the  first  group  show  a  high  data 
scatter  for  crack  speeds  greater  than  1  mm/s  which  is 
probably  due  to  an  inaccurate  method  of  determining  the 


U.  Summary 


1.  The  CTS  specimen  is  the  best  specimen  for  mixed  mode 
experiments  found  in  literature.  It  can  be  used  for  the 
complete  mixed  mode  range  from  mode  I  to  mode  II  and  Is 
easy  to  handle. 

2.  There  are  innumerable  fracture  criteria  in  literature 
which  are  partly  contradicting. 

3.  A  considerable  amount  of  experimental  mixed  mode  work  was 
found  in  literature.  The  problem  is  not  that  there  are  no 
experimental  results  available,  but  that  the  available 
results  are  not  reliable. 

V.  The  mechanical  properties  for  plexiglas  stated  in  the 
literature  differ  considerably.  Especially  the  influence  of 
the  straining  rate  on  the  modulus  of  elasticity  is  not 
clearly  understood. 

5.  The  critical  stress  intensity  of  plexiglas  depends  on  the 
crack  propagation  speed. 
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The  material  used  for  the  tests  was  Plexiglas  G  Acrylic 
Safety  Glazing,  produced  by  Rohm  and  Haas,  Philadelphia 
(PA).  All  the  specimens  were  cut  out  of  a  single  half  inch 
thick  plate.  As  the  plate  was  cast,  the  thickness  varied 
cons i derab I y . 


2.  Testing  Procedure 

Table  8  gives  the  dimensions  of  the  tensile  specimens. 
All  specimens  were  instrumented  with  an  extensometer  between 
point  2  and  3  to  measure  the  elongation.  Specimens  T1  and  T3 
were  additionally  instrumented  with  a  slightly  modified 
extensometer  to  measure  the  contraction  of  the  specimen  for 
the  evaluation  of  the  Poisson  rario.  The  results  were 
recorded  on  a  X  -  Y  plotter. 

For  the  measurement  of  the  modulus  of  elasticity  and  the 
Poisson  ratio,  the  specimens  were  stressed  to  a  level  of  11 
MPa  and  then  unloaded  with  a  specific  displacement  rate 
(table  9).  The  time  of  the  unloading  process  was  measured  to 
calculate  the  straining  rate.  Then  the  specimens  were 
loaded  until  they  ruptured.  The  location  of  frature  was 
noted. 

The  6  CT  specimens  were  prepared  and  tested  in  the 
following  way:  A  razor  blade  was  forced  into  the  machined 


chevron  notch.  Then  the  specimens  were  fatigue  precracked  at 
5  Hz  with  a  controlled  stroke  amplitude  of  0.15  nw"  and  an 
initial  minimum  load  of  about  65  The  final  crack  length 
typically  was  reached  after  12000  cycles.  The  crack  front 
was  straight.  The  dimensions  of  the  CT  specimens  are  given 
in  table  10.  The  specimens  were  instrumented  with  a  clip 
gage  according  to  ASTM  E399  (I)  and  a  load  -  COD  curve  was 
plotted  on  a  X  -  Y  plotter. 

The  machine  (MTS)  in  use  for  the  experiments  had  as 
lowest  load  range  1000  Ibf  with  a  least  count  of  1  Ibf.  In 
previous  tests  with  this  machine  it  was  found  that  the  load 
control  in  the  lowest  range  does  not  work  properly. 
Therefore,  displacement  control  had  to  be  used  for  all 
tests.  Displacement  control  was  used  by  most  other 
researchers  too  {15*21 ,26,3i»f  37)  ♦ 

The  testing  was  done  with  speeds  between  5«1  and  63*5 
mm/min.  In  displacement  control,  tbe  displacement  recorded 
on  the  plot  is  proportional  to  the  time.  The  time  until 
fracture  initiation  found  there  and  the  critical  stress 
intensity  K^was  used  to  calculate  the  stress  intensity  rate 
k.  The  exact  way  of  the  calculation  can  be  found  in  appendix 
2. 

As  all  the  machines  and  gages  in  use  were  in  US  customary 
units,  the  CTS  specimens  and  the  grips  were  slightly 
modified  to  meet  these  units.  The  used  dimensions  are  shown 
in  figures  h  and  5> 

Altogether,  16  CTS  specimens  were  produced  and  tested.  A 
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razor  blade  was  forced  into  the  machined  chevron  notch  to 
start  the  crack.  Then  the  specimens  were  fatigue  precracked 
at  5  Hz  with  a.  controlled  stroke  amplitude  of  0.I5  mm  and  an 
initial  minimum  load  of  about  125  N.  “The  specimens  were 
loaded  through  the  holes  at  the  side  of  the  notch  and  not 
through  the  middle  holes  as  suggested  by  Richard  {8}.  The 
crack  length  a  -  1.5in  which  equals  a/W  ■  0.5  was  typically 
reached  after  10000  cycles.  Two  specimens  were  further 
fatigued  to  a/W  ■  0.55  and  two  others  to  a/W  ■  O.65.  The 
crack  front  was  straight,  but  the  crack  tip  was  only  by 
chance  exactly  on  the  center  line  of  the  specimen.  Most 
samples  showed  deviations  of  up  to  0.1  mm,  some  up  to  0.2 
mm.  The  worst  case  was  specimen  CTS1  with  1.5  mm  deviation  _ 
(^.8%).  In  addition,  the  crack  of  this  specimen  was  slanted 
in  through  thickness  direction  by  about  U  degrees. 

The  testing  of  these  specimens  under  mixed  mode 
conditions  imposes  some  problems.  Richard  {8}  suggests  in 
accordance  with  ASTM  E399  {1}  to  use  load  control  and  choose 
F  such  that  the  stress  intensity  rate  K  in  mode  I  is  between 
0.55  and  and  2.75  Mpa/nT  per  second.  Richard  himself  used  K  “ 
3.3  MPa/irT/s.  He  further  suggests  that  the  loading  rate  F 

used  for  mode  I  should  be  contained  for  mode  II  and  mixed 

•  • 

mode  loading.  The  problem  now  is  that  K  changes  for  F  ■ 
constant  at  different  angles.  It  was  shown  by  earlier  works 
{95,101,102}  that  the  fracture  toughness  of  PMMA  is  rate 
sensitive.  Therefore  it  is  felt  that  the  suggested 
procedure  gives  no  clear  distinction  between  the  influences 
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of  mixed  mode  loading  and  those  of  the  rate  sensitivity. 

There  are  three  principally  different  ways  of  loading  the 
specimens  under  different  mixed  mode  conditions: 

1.  displacement  rate  d  -  constant 

« 

2.  loading  rate  F  -  constant 

3.  stress  intensity  rate  K  •  constant 

Naturally,  there  are  arbitrari ly  many  ways  of  loading 
without  any  constant  magnitude.  But  these  do  not  seem 
logical . 

Host  researchers  did  their  mixed  mode  testing  at  a 
constant  displacement  rate  {21,25.26,3^.37}.  As  they  used 
many  different  specimens,  the  influence  of  the  changing 
stress  intensity  rate  is  unknown.  But  for  a  CTS  specimen, 
the  use  of  a  constant  displacement  rate  is  questionable,  as 
in  this  case  the  stress  intensity  rate  under  mode  11  loading 
is  only  2S%  and  17%  of  that  under  mode  I  loading  for  a/W  ■ 
0.5  and  0.7  respectively.  This  can  be  seen  in  table  11.  The 
way  how  these  figures  were  calculated  is  given  in  appendix 
2. 

Richard  {8}  used  for  his  experiments  a  constant  loading 
rate.  This  is  already  better,  but  the  stress  intensity  rate 
under  mode  II  loading  is  still  only  28%  of  that  under  mode  I 
loading  for  a/W  *  0.7.  as  can  be  seen  in  table  12. 

Hence  in  the  new  experiments  it  was  tried  to  keep  the 
stress  intensity  rate  K  constant.  The  problem  is  that  the 
stress  Intensity  rate  is  no  property  of  the  testing  machine 
as  it  is  in  the  cases  of  the  displacement  or  loading  rate. 


To  relate  the  mixed  mode  and  mode  II  stress  intensity  rate 
to  that  of  mode  I,  a  fracture  criterion  is  necessary. 
Because  of  its  simplicity  and  its  good  agreement  to  previous 
results,  the  fracture  criterion  {8,76,77} 


was  used  where  it  was  supposed  that  K,,^  *  K,e  .  Then  the 
necessary  displacement  rate  was  calculated  as  stated  in 
appendix  2.  The  results  are  shown  in  table  I3. 

The  constant  stress  intensity  rate  condition  is  equal  to 
the  condition  of  constant  time  of  loading  which  was  used  by 

ni). 

The  tests  were  done  with  displacement  control.  A  load  - 
LLD  curve  was  plotted.  This  is  sufficient  in  the  case  of 
Plexiglas  which  has  such  a  low  modulus  of  elasticity  (3  GPa 
versus  210  GPa)  that  all  deformations  of  the  machine  and  the 
grips  can  be  neglected. 

After  the  testing,  pictures  of  the  crack  tip  with 
magnification  were  taken  to  measure  the  angle  of  crack 
initiation.  In  addition,  the  fracture  surface  was  examined. 

The  calculation  of  the  stress  intensities  K,  and  K„  was 
done  using  the  calibration  factors  given  by  Richard  {8},  as 
stated  in  appendix  3* 
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V  Results  and  Discussion 

1.  Mechanical  Properties 

The  plots  for  the  Poisson  ratio  and  the  modulus  of 
elasticity  were  straight.  For  the  calculation  of  the  modulus 
of  elasticity,  the  average  of  the  cross  -  section  area  at 
point  2  and  3  was  used.  In  addition  to  the  specific  tests 
of  the  modulus  of  elasticity,  this  property  was  evaluated 
out  of  the  plots  of  the  rupture  tests.  The  results  agree 
with  a  deviation  of  3%  or  less. 

The  specimens  broke  at  different  locations  and  this  may 
be  because  of  stress  concentrations  due  to  machining.  It  was 
not  possible  to  get  consistent  values  for  the  ultimate 
stress  d^.  Hence  the  yield  strength  was  evaluated  which 
shows  little  data  scatter.  The  results  are  shown  in  table 
1i».  The  average  of  the  Poisson  ratio  was  calculated  to  be 
0.36. 

A  comparison  with  table  7  shows  that  the  new  results 
agree  well  to  those  stated  in  the  literature.  The  new 
results  show  slightly  increasing  values  for  the  properties 
with  increasing  straining  rate.  The  growth  is  about  11%  in  E 
and  for  an  eight  times  higher  strain  rate,  as  shown  in 
figure  6.  In  many  cases,  this  might  not  be  neglegible.  Hence 
the  reported  increase  of  E  with  increasing  strain  rate  could 
be  confirmed. 

Equation  (7)  given  by  Williams  (95)  yields  much  lower 
value's  for  E  than  found  in  the  new  experiments,  but  this  is 


due  to  the  definition.  Williams  defined  the  modulus  of 
elasticity  as  the  secant  between  the  origin  and  the  point  of 
3.5^  strain,  while  usually  the  initial  slope  of  the 

straining  curve  is  used.  Equation  (7)  predicts  a  15% 
increase  in  E  for  the  range  covered  by  the  new  experiments. 

2.  Fracture  Toughness  at  Varying  Stress  Intensity  Rates 

The  load  -  COD  plots  for  the  CT  specimens  were  almost 
linear  and  all  did  fulfill  the  requirements  of  ASTM  E399. 
Pre i ns  tab i 1 i ty  crack  growth  was  not  observed,  but  up  to  the 
stress  intensity  rate  K  «  2.15  MPa\^/s,  the  crack 

propagation  speed  was  low  (range  of  millimeters  per  second). 
The  fracture  surface  therefore  showed  the  river  markings 

which  are  typical  for  low  crack  speeds  {9I4}  .  Only  specimen 
CT6  (K  ■  3*78  MPa/m7s)  showed  the  smooth  fast-fracture 
surface. 

Table  15  shows  the  results  for  the  CT  specimens.  The 
results  are  consistent  with  Marshall  et  al.  {I03}  who  gave  a 
K|j  of  about  1.2  MPa/m  for  crack  speeds  in  the  range  of 
millimeter  per  second. 

Then  k  CTS  specimens  were  tested  under  mode  I  loading  at 
different  stress  intensity  rates.  Again,  the  load 
displacement  curves  were  almost  linear  and  met  the 

requirements  of  ASTM  E399*  Only  specimen  CTS7  (K  ■  0.23 

MPa/tn/s)  showed  the  rough  slow  -  fracture  surface,  the  other 
specimens  had  a  smooth  fracture  surface. 
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The  results  are  shown  in  table  16.  All  K,^  values  from 
tests  with  CIS  specimens  are  significantly  greater  than 
those  from  CT  specimens  at  the  same  stress  intensity  rate  K. 
This  can  be  seen  in  figure  7>  Also  the  change  from  the  rough 
to  the  smooth  fracture  surface  happened  at  a  much  lower  K 
for  the  CTS  specimens  than  for  the  CT  specimens.  The  reason 
for  that  behavior  is  not  clearly  understood.  It  might  be 
that  the  fracture  process  happens  faster  in  a  CTS  specimen 
as  this  specimen  is  stiffer  and  more  el ast i c_ energy  is 


stored  at  the  onset  of  fracture.  The  crack  separation  rate 
for  the  two  specimens  might  be  different.  Chiang  and  Miller 
{1*3}  calculated  the  crack  separation  energy  rate  for  ACCT 


and  SENS  specimens  and  found  it  to  be  different.  They  - 


attributed  the  difference  to  the  effect  of  hydrostatic 
pressure  which  influences  the  plastic  zone  around  the  crack 


tip.  Whatever  the  reason  may  be,  it  can  be  concluded  that 
there  is  no  unique  relationship  between  the  stress  intensity 
rate  and  the  critical  stress  intensity. 


As  the  crack  speed  was  not  measured  in  the  new 
experiments,  the  statement  that  the  relation  between  the 


crack  speed  and  the  critical  stress  intensity  is  unique  and 
independent  of^  the  specimen  geometry  {95}  could  not  be 
confirmed.  But  an  analysis  of  the  slope  of  the  test  plot 
after  crack  initiation  showed  that  a  higher  K,  was 


correlated  to  a  steeper  falling  curve  which  means  a  faster 
crack  speed. 

Although  the  values  for  the  critical  stress  intensity 


-■  -'-I 

I  I 


obtained  with  CTS  specimens  are  generally  higher  than  those 
found  with  CT  specimens,  they  are  still  within  the  range 
found  in  literature.  The  question  remains  why  Richard 
got  a  30%  higher  value  with  the  same  specimen  for  the  same 
material. 

The  first  reason  for  this  difference  is  that  he  used  a 
much  higher  stress  intensity  rate.  A  second  reason  may  be 
that  Richard  had  curved  crack  fronts,  as  can  be  seen  in 
figures  86  and  87  in  {8}.  The  curved  crack  front  is  due  to 
the  higher  frequency  used  for  fatigue  precracking  and  was 
found  in  earlier  tests  by  Pou'lose  too.  The  higher  frequency 
causes  a  heating  of  the  specimen  which  is  probably 
responsible  for  the  effect.  As’ was  found  by  Towers  and  Smith 
{III4}  in  a  numerical  analysis  of  the  CT  specimen,  a  curved 
crack  front  which  still  meets  ASTM  E399  requirements  can 
have  an  up  to  28%  higher  fracture  l-oad  than  an  equivalent 
straight  crack  front. 

In  summary  it  can  be  said  that  the  critcal  stress 
intensity  of  PMMA  increases  slightly  with  increasing  stress 
intensity  rate.  For  small  variations  of  the  stress  intensity 
rate,  the  variation  of  the  critical  stress  intensity  is  not 
outstanding  of  natural  data  scatter,  but  a  change  of  K  by  a 
factor  2  or  more  is  not  neglegible. 

3.  Fracture  Toughness  under  Mixed  Mode  Loading 

Beside  two  specimens  in  mode  I  which  were  already 


mentioned  in  the  previous  chapter,  all  specimens  were  loaded 
at  a  stress  intensity  rate  of  1.0  MPa\/m/s  with  a  maximum 
deviation  of  15%,  as  can  be  seen  in  table  I7.  The  deviations 
are  due  to  inaccurancy  in  the  calculation  of  the 
displacement  rate  and  scatter  of  the  stiffness  in  different 
specimens. 

Out  of  the  16  tested  specimens,  three  (CTS1,  CTS12, 
CTSI3)  had  to  be  rejected.  All  of  them  showed  the  same 
effect:  The  crack  growth  started  slowly  at  one  side  of  the 
crack  and  moved  up  to  3  inm  before  the  instability  condition 
was  reached.  The  fracture  surface  was  slanted  in  thickness 
direction.  In  all  cases,  some  geometrical  deviations  could 
be  found.  As  two  of  the  rejected  specimens  were  tested  under 
mode  II  loading,  it  seems  that  this  condition  is  especially 
sensitive  to  deviations  of  the  specimen  geomety. 

The  fracture  surface  of  all  other  specimens  was  smooth 
(mirror  surface)  and  showed  none  or.  very  little  (<  0.5  mm) 
slow  crack  growth  which  then  was  neglected.  The  results  of 
the  mixed  mode  experiments  are  given  in  table  I8, 

K,j  was  calculated  as  the  average  of  rhe  results  of 
specimens  CTS2,  CTS9,  and  CTSll  and  had  the  value  1.26 
ttPa/m.  The  evaluation  of  is  a  little  vague.  First  of  all, 
out  of  the  5  tested  specimens  2  had  to  be  rejected.  The 
result  of  the  remaining  3  is  not  that  consistent  either.  The 
mean  value  of  the  two  specimens  with  a/W  =0.5  is  1.39  MPa/nT 
which  yields  ■  1.10.  Taking  the  specimen  with  a/W  = 
0.65  in  account,  the  average  of  the  three  specimens  is  K,.  ■ 
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1.35  MPa/in  which  gives  K|,j/K,j  *1.07.  In  any  case,  K,|£/K^^>  1. 
This  is  much  more  than  the  values  found  in  literature  for 
PnttA.  As  already  mentioned,  most  researchers  found  K,|^/K|jto 
be  around  O.9  (table  6).  Only  Phadke  et  al.  {21}  found  Ki^/K, 
to  be  up  to  1.0,  depending  on  which  mode  II  specimen  they 
used  for  the  evaluation  of  Kjj^.  But  the  new  results 
correlate  well  to  results  found  for  steels  (table  3), 
aluminum  alloys  (table  4)  and  different  other  materials 
(table  5)-  For  many  materials  around  1.1  was  found. 

The  reason  for  this  discrepancy  between  the  results  for 
PMMA  and  other  materials  might  be  that  most  researchers 
neglected  the  influence  of  the  stress  intensity  rate  on  the 
critical  stress  intensity  of  PMMA.  Richard  {4,8}  found  Kn/K,^ 
■  0.93,  but  as  pointed  out  ealier,  his  stress  intensity  rate 
for  mode  II  was  just  about  30^  of  that  for  mode  I.  Assuming 
that  the  dependency  of  the  critical  stress  intensity  on  the 
stress  ■  intens i ty  rate  is  the  same  for  mode  I  I  as  it  is  for 
mode  1  and  also  assuming  that  the  dependency  can  be 
described  by  a  line  in  the  semi  log  -  plot  figure  7,  the 
stated  value  for  K,|g  is  about  1%  too  low.  A  correction  for 
this  would  result  in  K|,/K|j  of  about  1. 

A  problem  that  could  have  influenced  the  new  results  is 
crack  closure.  Especially  for  a/W  <  0.55.  there  is  some 
negative  stress  intensity  Kj  for  cx *  90°  (see  figure  5  in  {2} 
or  figure  78  in  {8}).  This  might  have  had  an  influence  on 
the  results  with  a/W  *  0.5* 

For  a/W  »  O.65  and  c*  *  90°,  K  is  virtually  zero.  But  as 
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the  tested  specimens  had  an  extremely  long  fatigue  crack 
(16.5  ™ji)  .  even  a  very  small  negative  K,  may  influence  the 
results.  This  is  confirmed  by  Bauerle  (115)  who  made  a  new 
FEH  analysis  of  the  CTS  specimen.  Although  he  found  for  a  » 
90”  K|  to  be  almost  .  zero  at  the  crack  tip,  K|  got 
increasingly  negative  along  the  crack  flanks.  As  the 

surfaces  of  the  fatigue  crack  are  in  touch  with  each  other, 
they  get  pressed  together  due  to  the  negative  K, ,  causing  an 
unknown  error  in  the  results.  The  crack  closure  problem  may 
be  increased  by  the  inevitable  out  of  middle  location  of  the 
crack  tip.  It  is  therefore  felt  that  a  mode  II  specimen 
should'have  a  small  amount  of  positive  K,  to  reduce  crack 
closure  and  friction  problems.  This  requirement  is  met  by 
the  modified  CTS  specimen  {1*2,116}  (see  figure  8). 

As  can  be  seen  in  figure  9»  the  mixed  mode  fracture 
cr i ter  ion 


suggested  by  several  authors  (8,76,77)  gives  a  good  fit  to 
the  experimental  results.  This  criterion  can  be  transformed 


into  {8} 


.  2C  Kj  +  1/  kI  +  4(o,Kjj)2'  ] 


where  is  the  comparative  stress  intensity  factor  and  CK= 
(K^/K,J  .  The  safe  -  design  condition  is  then 


V  Ic 


A  look  at  table  2  reveals  that  just  the  criterion  by 
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and  Nuismer  (55) i  while  the  other  criteria  make  more  or  less 
different  predictions.  The  measured  fracture  angle  with  the 
Erdogan  and  Sih  criterion  and  the  empirical  criterion  by 
Richard  {8}  are  shown  in  figure  12.  The  measured  fracture 
angles  of  beryllium  (V  -  0.0)  {1»7}  are  given  in  the  same 
figure.  It  is  remarkable  that  they  are  very  close  to  the  new 
results  despite  the  big  difference  in  the  Poisson  ratio.  It 
seems  that  the  angle  of  fracture  initiation  is  independent 
of  the  Poisson  ratio. 


Uf.'. 
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VI  Conclusions 


1.  The  thickness  of  the  CTS  specimen  should  be  symbolized  by 
B  and  not  by  t.  because  t  can  easily  be  held  for  a  symbol 
for  the  time.  B  is  used  to  describe  the  thickness  of  the  CT 
specimen  {!},  the  CS  specimen  arid  the  new  compact 
shear  (NCS)  specimen  {11?}.  and  there  is  no  reason  why  this 
should  be  changed. 

2.  The  testing  of  the  specimens  should  be  done  in  a  way  that 
the  comparative  stress  intensity  rate  Ki,  is  constant  for  ait 
mixed  mode  conditions  and  not  with  a  constant  loading  rate. 
To  calculate  K„ ,  the  criterion 


^  -  iCKi  » 


as  a  first  approximation  is  sufficient.  Even  if  it  turns  out 
that  the  implied  assumption  K,^  »  K,^  is  not  accurate,  the 
difference  in  the  stress  intensity  rate  over  the  mixed  mode 
range  will  be  much  less  than  with  'just  taking  a  constant 
loading  rate. 

3.  As  there  are  some  crack  closure  problems  in  pure  mode  II, 
the  crack  should  be  as  short  as  possible.  There  should  be 
set  not  only  a  lower  limit  for  the  crack  length  {8}  but  also 
an  upper.  In  order  to  cut  the  length  of  the  touching 
surfaces,  it  has  to  be  considered  whether  a  straight  through 
notch  should  be  used  instead  of  a  chevron  notch.  Maybe  the 
original  CTS  specimen  should  be  replaced  at  all  by  the 
modified  CTS  specimen  {1*2}  for  mode  II  testing.  The  modified 
CTS  specimen  has  a  slightly  positive  K,  for  cx  ■  90®  (116). 
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U.  More  research  is  needed  to  find  out  the  influences  of 


deviations  of  the  specimen  geometry.  It  is  not  enough  just 
to  keep  the  machining  tolerances,  as  the  fatigue  crack  tends 
to  grow  away  from  the  center  line.  Relatively  small 
deviations  already  showed  a  significant  effect. 

5.  Both  the  mechanical  properties  and  the  fracture  toughness 
of  PMMA  are  slightly  sensitive  to  the  loading  rate.  For 
changes  of  more  than  factor  2  in  the  loading  rate,  the 
deviations  are  not  neglegible. 

6.  The  fracture  criterion 


K 


Ic 


K 


T  2 


II 


K 


lie 


=  1 


gives  a  good  fit  to  the  experimental  found  mixed  mode 
fracture  initiation. 

7.  The  ASME  Boiler  and  Vessel  Code  gives  a  reasonable  fit  to 
the  experimental  data  up  to  a  crack  inclination  of  1*5°.  The 
predictions  for  steeper  inclinations  are  not  accurate  but 
they  are  on  the  safe  side. 

8.  The  criterion  of  maximum  tangential  stress  proposed  by 
Erdogan  and  Sih  gives  the  best  fit  to  the  experimentally 
found  angle  of  fracture  initiation.  It  seems  that  this  angle 
is  independent  of  the  Poisson  ratio. 
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Fracture  17,  RIO5-RIO7  (I98I) 


0.0 

0.268 

0.577 
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1.732 

3.732 

CO 


mode  I 


mixed 

mode 


mode  II 


Table  1  Relation  between  p,  m,  and  n 


Poisson 


mode  II 


Erdogan  and  Sih  [13 J 
Nuismer  [551 
Hussain  et  al.  [381 
Ametoy  et  al.  [5^1 
Radaj  and  Heib  [57 » 581 
Sih  [16,17,591 
Di  Leonardo  [60I 
Irwin  [61,621 
J  -  integral  [6 3-651 


legend t  i  =  independent 

-  =  no  prediction 

Table  2  Prediction  of  different  fracture  criteria  for 
mode  II 


Ueda  C25] 
Shah  [18] 


Chiang  [4]] 


2 

Material 

3 

Spec 

4 

T["C] 

M  SM41 

ACNT 

RT 

M  SM41 

CNBT 

RT 

H  4340 

CCT 

RT 

H  4340 

Tube 

RT 

H  4340 

ACCT 

-140 

H  l8MnMoV 

ACCT 

RT 

H  GC-4 

SECB 

RT 

S 

ACNT 

RT 

.92  110 1  - 


7 

^Ic 

8 

''lie 

- 

- 

- 

80.4 

- 

- 

88.3 

45.7 

49.7 

- 

- 

49.6 

- 

- 

1)  doubtful  method  of  determining 
Table  3  Mixed  mode  results  for  several  steels 
Legend  to  tables  3  "to  5 


column  2  M  =  mild  steel 

H  =  high  strength  steel 
S  =  spring  steel 

column  3  specimens,  see  chapter  III  1. 
column  4  RT  =  room  temperature 
column  6  A  =  number  of  the  tested  specimens 
column  7  +  8  values  in  MPa/m 

♦  =  value  found  in  a  figure 
-  =  not  stated  in  the  paper 


■m 


V  •V 


1)  shear  panel 

2)  with  a  CCT  specimen, 

3)  with  a  CCT  specimen, 


Kordisch  [27]] 

Richard 

Theocaris  C3l3 

Theocaris  [331 

Lemant  [2$! 

Awaji  and 
Sato  C46] 


Material 

Glass 

PU 

Araldit  B 

PCBA 
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Bronze 

Zinc 

G^SA24  // 
G  SA24 
G  7477 
G  SMl-24 
Plaster 
Marble 


Spec  T(‘G]  n^i^  n^^ 


ACCT  -  0  0.92* 

AGCT  -  0  0.85‘ 

GTS  RT  0  1 

AGNT  RT  0  0,92 

AGNT  RT  0  0.92 

AGGT  -196  0  0.63 

BD  RT  0  1 

BD  RT  0  1 

BD  RT  0  1 

BD  RT  0  1 

BD  RT  0  1 

BD  RT  0  1 


7 

8 

9 

^Ic 

Ktt 

lie 

Ktt 

lie 

Kt 

Ic 

0.59 

1.46 

2.48 

7.0® 

tm 

— 

— 

0.72 

0.83 

1.15 

0.69 

0.76 

1.11 

0.94 

1.09 

1.16 

0.81 

0.89 

1.09 

0.13 

0.15 

1.14 

0.93 

1.05 

1.13 

1)  calculated  from  given  data 

2)  three  different  types  of  bronze  ■ 

3)  G  a  graphite 

legend  see  table  3 


Table  5  Mixed  mode  results  for  various  materials 


Table  6  Mixed  mode  fracture  data  for  PM.ViA  stated  in  literature 


Footnotes  to  table  6 

1)  load  control  F  =  4kN/s 

2)  spherical  shells  Sf  178  -  270  mm  under  water  pressure 

3)  constant  time  of  loading  15  s 

4)  constant  time  of  loading  75  s 

5)  loaded  by  running  water 

6)  only  47  specimens  evaluated  for  Kj,  because  of 

pre instability  crack  growth 

7)  in  table  1  of  that  paper,  B  = is  stated 

8)  the  ligament  in  front  of  the  crack  tip  was  reduced  to 

2.5  mm.  The  results  are  probably  not  valid  because  of  the 
uncertain  state  of  stress  at  the  crack  tip  and  possible 
residual  stresses  due  to  grinding. 

9)  slov/  crack  growth  only  in  mode  I 

10)  constant  stress  intensity  rate  K=0.2  MPa\/m/s 

11)  6  specimens  at  Pjj/Py=0,  I7  at  p^^/Py  >  0 


Legend  to  table  6 

column  3  B  r  thickness  of  the  specimen  in  [ram] 
column  4  introducing  of  the  crack 
F  =  fatigue  precracking 
R  =  razor  blade  forced  into  the  notch 


N  =  machined  notch 


column  5  cl  =  displacement  rate  in 

column  6  did  pre instability  crack  growtth  occur? 

N  =  no,  Y=  yes,  F  =  in  some  specimens 
column  8  A  =  number  of  the  tested  specimens 


*  =  value  found  in  a  figure 
- =  not  stated  in  the  naner 


Reference 


t 


e 


Reference 

£ 

[GPa] 

V 

^ult 

[MPa] 

^^0.2 

[MPa] 

1 — 1 
+>  w 

• 

e 

Dubbel  [87] 

2.7  -3-7 

- 

50-77 

- 

- 

- 

Encyclopedia  [88] 

2.6  -  3*1 

0.35 

48  -  76 

- 

- 

Mascia  [89] 

3.12 

0.33 

- 

- 

- 

Ewing  [11] 

- 

0.35 

- 

- 

- 

- 

Richard  [4] 

- 

0.30 

- 

- 

- 

- 

Richard  [7] 

- 

0.36 

- 

- 

- 

- 

Swing  [37] 

2.0 

0.43 

- 

- 

- 

- 

Buhelt  [90] 

3.41 

0.33 

- 

- 

- 

Preece  [91] 

- 

0.35-0.37 

- 

- 

- 

- 

Zimmermann  [92] 

5.0 

0.35 

- 

- 

- 

high 

Xausch  [93] 

30 

- 

65 

- 

- 

- 

Leevers  [9^] 

2.9 

- 

80 

- 

- 

Rohm  i  Haas 

2.90 

0.35 

72 

- 

- 

Theocaris  [33] 

3.4 

0.34 

72.9 

64.3 

- 

- 

Table  7  Mechanical  properties  of  pr-".!A  stated  in  literature 


7.3"  (185.4) 
4.1"  (104.1) 


Spec 

width 

thickness 

1 

2 

3 

■■ 

1 

2 

3 

■■ 

T1 

12.5 

12.6 

12.6 

12.7 

12.0 

12.0 

12.1 

12.2 

T2 

12.5 

12.6 

12.5 

12,7 

12.1 

12.1 

12.0 

11.9 

?3 

12.7 

12.3 

12.9 

13.0 

12.1 

12.1 

12.0 

11.9 

TU 

12. i!- 

12.^ 

12.5 

12.5 

12.2 

12.1 

12.1 

12.0 

Table  8  Dimensions  of  the  tensile  specimens 


Spec 

W 

[mm] 

a 

Cmm] 

a/V/ 

B 

Cmm] 

CTl 

28.2 

0.555 

12.0 

CT2 

25.1 

0.495 

11.7 

CT3 

50.8 

24.9 

0.490 

11.7 

CT4 

(2.0") 

25.4 

0.500 

11.7 

CT5 

25.1 

0.495 

11.7 

CT6 

31.0' 

0.610 

11.7 

Table  10  Dimensions  of  the  CT  specimens 


a/'A’  =  0.5  I  a/W  =  0.7 


1.0 

0.96 

0.73 

0.67 

0.56 

0.56 

0.59 


1.0 

0.9^ 

0.71 

0.56 

0.40 

0.33 

0.29 


1.0 

0.93 

0.69 

0.51 

0.34 

0.24 

0.17 


Table  11  Relation  between  d,  F,  and  K  for  a 


constant  displacement  rate  d 


Stsec. 


K 


MPa/m 


Lmin  J 


s 


CTl 

CT2 

CT3 

CT4 

CT5 

CT6 


5.1 

0.32 

5.1 

0.33 

10.2 

0.71 

15.3 

1.09 

30.5 

2.15 

63.5 

3.78 

[Nl 

Kt 

Ic 

[  MPa/nT] 

244 

1.04 

290 

1.05 

315 

1.13- 

316 

1.16 

317 

1.15 

232 

1.25 

Table  15  Testing  conditions  and  results  for  the 
CT  specimens  (mode  I) 


Spso. 

■ 

F 

c 

[JcN] 

Kic 

[MPa/nT] 

CTS7 

CTS9 

CTS  2 

CTSll 

2.5 

10.9 

10.9 

12.7 

0.23 

0.96 

1.00 

1-25 

1.01 

1.12 

1.19 

1.17 

1.12 

.  1.23 

1.29 

1.27 

Table  16  Testing  conditions  and  results  for  the 
CTS  specimens  (mode  I) 


oc 

Spec. 

IMI 

[ MPa/m] 

^11 

[rJPa\/^ 

-0 

0°  . 

CTS9 

1.12 

1.23 

0 

0° 

CTS2 

1.19 

1.29 

0 

0° 

CTSll 

1.17 

1.27 

0 

0° 

15° 

CTS12 

invalid 

o 

o 

CTS14 

1.52 

Q| 

o 

o 

CTS5 

1.22 

WM 

32° 

45° 

CTSIO 

1.52 

1.14 

0.57 

34° 

60° 

GTSi: 

1.60 

0.85 

0.76 

49° 

75° 

CTS8 

2.20 
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Figure  8  Modified  CTS  specimen  C^2] 
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The  citerion  {8,76,77) 
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can  be  transformed  In  the  following  way:  Multiplying  both 
sides  by  and  adding  (1/4)  K,*,  one  obtains 
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Introducing  o,«  K^/Kj^  and  extracting  the  root,  this  finally* 
yields 


he  “  ) 


The  ASME  Boiler  and  Vessel  Code  (79)  gives  the  criterion 

s  (J'VTcasinp 

With  K|  ■  tfVica’  sin*p  ‘for  an  lnfinit«  plate,  one  obtains 

.  Kj(sinp)"^-^ 

With  p  *  arctan (K, /K|,)  and  applying  the  respective 
trigonometrical  formula,  this  yields 


Appendix  2;  Calculation  of  the  Stress  Intensity  Rate  K 


The  calculation  is  done  for  specimen  CT4  as  an  example. 

This  specimen  broke  at  the  critical  stress  intensity  > 

1.l6  MPa/nT  and  showed  there  a  COO  of  0.439  mm.  Neglecting 

the  deformations  of  the  crack  flanks,  the  point  of  measuring 

the  COD  has  the  distance  ■  1.625"  to  the  crack  tip,  while 

the  load  line  has  L2*  1.00".  LID  ar\d  COD  are  related  by  the 

simple  geometrical  relation 

LLP  _  ^ 

COD  '  L, 

which  yields  LLD  *  0.270  mm  at  the  initiation  of  fracture. 

In  displacement  control,  the  LL6  is  related  to  the 
displacement  rate  d  by  _ 

LLD  =  ht 


with  which  the  time  until  f 
■^c  =  ^ 


racture  can  be  calculated: 


The  average  stress  intensity  rate  K  then  is  calculated  by 

^IC  r-, 

K  =  =  1.09  MPa\^/s 

The  calculation  for  the  other  CT  specimens  is  done  in  the 
same  way. 

For  the  CTS  specimens,  the  LLD  was  recorded  and  therefore 


the  conversion  from  COO  to  LLD  is  not  necessary. 
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It  is  intented  to  keep  K  constant  for  the  mixed  mode 
loading.  The  problem  is  that  a  relation  between  the  machine 
parameter  d  and  K  has  to  be  found.  The  loading  rate  F  is 
related  to  d  by 

p  =  ai 


where  i  is  the  slope  of  the  load  ~  LLD  plot.  The  relation 

•  • 

between  K  and  F  is  given  by 
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where  W  is  the  width  of  the  specimen,  B  is  the  thickness, 
and  Yy  the  comparative  calibration  factor.  Richard  {8}  made 
the  FEtt  calculation  and  got  the  calibration  factors  Y,  and 
Y|  for  the  mode  I  and  mode  It  component  respectively  (see 
appendix  3) •  These  factors  have  to  be  related  by  a  fracture 
criterion.  The  criterion 
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II 
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lie 
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was  r.' osen  and  the  assumption  K,,^  =  K,^  was  made.  With  this 
criterion,  Yy  can  be  calculated  by 


Yy  =  ) 


This  has  been  done  for  a/W  ■  0.5  and  a/W  =  0.7* 

The  last  remaining  unknown  is  the  slope  i  which  has  to  be 
determined  experimentally.  For  this  reason,  one  specimen  was 
loaded  to  about  N  and  then  unloaded,  both  at  a 
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Appendix  3 

To  calculate  K,  and  K,,  Richard  {8}  gives  the  formulas 


If  —  y 

^  WB 


^II  ~  ^  VB  ^11 


where  Y,  and  Yg  are  the  calibration  factors  for  the  CTS 
specimens  evluated  by  a  FEM  calculation.  Yj  and  Yg  have  the 
following  values  {8} 
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ABSTRACT 


\ 


A  center-cracRed  panel  of  707 5 -Aluminum  alloy  which  has 
overaged  from  the  T651  condition  was  loaded  in  tension  and 
subsequently  unloaded  to  zero  applied  load.  The  permanent 
surface  deformation  was  measured  close  to  the  intersections  of 
the  crack,  front  with  the  free  surfaces.  The  permanent 
deformation  (being  a  good  indicator  of  the  extent  of  plastic 
deformation)  was  used  to  measure  the  accuracy  of  finite  element 
analyses .-  ^ 

The  same  specimen  was  modeled  using  20-node 
three-dimensional  isoparametric  elements.  A  fully  incremental 
elastic-plastic  formulation  was  employed  in  the  stress 
analysis.  The  residual  surface  deformations  after  unloading 
were  compared  to  the  experimental  results. 

The  average  experimental  results  compare  quite  favorably 
with  the  finite  element  predictions.  The  average  results  were 
employed  to  minimize  the  influence  of  material  inhomogeneity, 
load  misalignment,  and  lack,  of  symmetry  in  the  fatigue  crack.. 

The  scatter  in  the  results  from  measuring  the  different  sides  is 


INTRODUCTION 


Over  the  past  35  years,  the  field  of  fracture  mechanics  has 
evolved  and  developed  into  an  important  and  useful  tool  for  the 
design  of  engineering  components  and  structures.  Several  major 
problems  dealing  with  the  criticality  of  cracks  in  engineering 
components  and  structures  can  now  be  answered  with  great 
accuracy.  Specifically,  problems  involving  straight  cracks  in 
brittle  materials  undergoing  Mode  I  deformation  only  can  be 
accurately  predicted.  Most  problems  which  occur  in  practice, 
however,  involve  materials  which  are  ductile  in  the  loading 
applications  for  which  they  are  employed.  Many  cracks  are  also 
initiated  in  sites  which  involve  complicated  loading  which 
involve  more  than  a  single  fracture  mode.  Finally,  most 
problems  arising  in  application  involve  geometries  which  can  not 
be  accurately  approximated  two-dimensionally .  The  major 
research  in  fracture  mechanics  today,  therefore,  is  geared 
toward  addressing  the  issues  of  ductility,  mixed-mode  loading 
and  three-dimensionality. 

The  issue  of  plasticity  and  ductility  in  fracture  specimens 
has  long  been  a  concern  of  researchers.  Many  attempts  have  been 
made  to  propose  fracture  criteria  which  account  for  ductility 
and  to  develop  numerical  tools  "o  perform  stress  analyses.  The 
development  and  refinement  of  the  finite  element  method  has 


greatly  aided  the  progress  in  this  area.  Many  two-dimensional 
studies  have  been  performed  using  elastic-plastic  finite  element 
modeling.  While  much  of  the  early  work,  has  been  demonstrated 
erroneously,  several  accurate  computational  procedures  are  now 
available.  The  area  of  failure  prediction  has  not  been  as 
successful  as  the  area  of  stress  analysis.  No  viable  ductile 
fracture  criteria  have  been  proposed  which  pass  the  tests  of 
specimen  and  geometry  independence,  consistent  and  theoretically 
sound  formulation,  and  reproducibility.  The  best  that  can  be 
said  for  the  existing  criteria  is  that  for  limited  realms  of 
applicability  (usually  vary  within  10-15%  of  the  range  of 
brittle  criteria),  the  proposed  methods  offer  conservative 
estimates  for  failure  loads  which  are  not  as  strict  as  the 
brittle  predictions.  It  is  important  to  recognize  at  the  outset 
that  elastic-plastic  fracture  parameters  (e.g.,  J-integral, 

CTOD,  CMOD,  etc.)  can  either  be  shown  to  be  theoretically 
invalid  for  true  plasticity  problems,  or,  are  simply 
experimental  observations  which  not  do  pass  the  test  of  specimen 
and  geometry  independence. 

During  the  past  three  years,  the  authors  have  focused  their 
research  on  addressing  the  three-dimensional  aspects  of  ductile 
fracture.  A  major  first  step  has  been  the  development  of  an 
accurate  and  theoretically  consistent  computational  approach  to 
the  stress  analysis  of  three-dimensional  fracture  specimens.  In 


a  series  of  recent  papers  [1,2,3],  the  effect  of  specimen 
thick.ness,  material  hardening  characteristics  and  mesh 
characteristics  have  been  investigated.  These  results  give  much 
insight  into  the  necessary  properties  for  ductile  fracture 
criteria.  While  no  new  criteria  have  emerged  to  date  (either 
from  the  authors  or  others),  the  groundwork,  for  analyzing  and 
assessing  criteria  has  been  established. 

A  major  problem  with  ductile  fracture  problems  in 
three-dimensions  is  the  establishment  of  the  accuracy  of  the 
analysis.  Convergence  studies  are  extremely  costly  and  only 
show  the  consistency  of  the  approach.  They  in  no  way  guarantee 
agreement  with  the  behavior  of  real  materials.  To  address  this 
problem,  the  study  presented  in  this  paper  compares  the 
deformation  predicted  from  a  full  three-dimensional  incremental 
plasticity  finite  element  analysis  to  the  deformations  measured 
in  the  laboratory.  A  center-crack.ed  panel  was  chosen  for  the 
study  for  two  reasons:  first,  the  authors’  previous  studies 
have  been  performed  on  center-cracked  panels  and  second,  the 
specimen  is  easier  to  model  with  finite  elements  as  the  effects 
of  the  loading  holes  are  easier  to  account  for  (by  using  an 
accurate  gauge  length).  While  a  successful  comparison  does  not 
guarantee  the  accuracy  of  any  given  study  other  than  the 
present,  it  is  the  most  rigorous  way  of  establishing  the 
validity  of  the  approach  and  demonstrating  the  qualitative 
agreement  of  the  predictions  made  previously  with  the  behavior 
of  real  fracture  specimens. 
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ELASTIC -PLASTIC  FINITE  ELEMENT  FORMULATION 


The  stress  analysis  in  this  study  is  performed  utilizing  the 
finite  element  method  to  solve  the  basic  elastic-plastic 
governing  equations  for  the  deformation  of  continuum  solids. 

J2  flow  theory  plasticity  is  employed  with  the  standard 
associative  flow  law.  The  Newton-Raphson,  or  Tangent-Stiffness 
approach  is  employed  in  the  finite  element  formulation  to  handle 
nonlinearities.  The  Updated  Lagrangian  coordinate  system  is 
employed  to  handle  finite  strains.  The  formulation  of  all 
equations  is  outlined  in  this  section. 

The  J2  flow  theory  of  plasticity  assumes  that  the  material 
in  question  yield,  or  starts  deforming  plastically  when  the 
'•effective  stress”  (or  von-Mises  stress)  reaches  a  critical 
value  (called  the  yield  stress)*.  Prior  to  the  onset  of 
plasticity,  the  material  is  assumed  to  behave  linear 
elastically.  Subsequently,  the  deviatoric  stress  components  are 
related  to  the  deviatoric  strain  rate  through  the  tensor  relation 
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where  eij  are  the  deviatoric  strain  rates  given  by 
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Sij  are  the  current  deviatoric  stress  components  given  by 
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Sij  are  the  deviatoric  stress  components  measured  relative  to 
the  current  yield  surface  center  given  by 

^ij  “  ^ij  -  ^ij 

a^j  are  the  coordinates  in  stress  space  of  the  current  yield 
surface  center,  are  the  Cauchy  stress  components,  are 
the  "true"  strain  components  (discussed  in  a  subsequent 
section),  o_  is  the  effective  stress  given  by 
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e  7  ij  ij 


(5) 


and  Oq  is  the  effective  stress  measured  relative  to  the  current 
yield  surface  center 
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The  function  f(o_)  is  derived  from  the  uniaxial  stress-strain 
curve  and  is  consistent  with  the  Associated  Plasticity  Theory  (a 
complete  discussion  is  given  in  Reference  [4]),  Derivation  of 
f(a_)  for  a  multilinear  representation  of  the  stress-strain 
curve  will  be  discussed  subsequently. 

For  plastic  strains  which  are  incompressible,  the 
hydrostatic  plastic  strain  rate  is  zero.  The  total  hydrostatic 
strain  rate,  therefore,  is  related  to  the  related  to  the 
hydrostatic  stress  rate  by 

c  »  1  -  2v  a  (7) 

PP  E  pp 

Engineering  materials  exhibit  different  types  of  uniaxial 
hardening  behavior  when  subsequently  unloaded  after  being 
plastically  deformed.  Generally,  the  behavior  falls  between  two 
extremes  called  kinematic  and  isotropic  hardening.  The  uniaxial 
representation  of  these  behaviors  for  a  bilinear  material  are 
shown  in  Figure  1,  To  allow  for  various  hardening  behaviors  in 
the  multiaxial  formulation,  the  yield  surface  is  permitted  to 
move  and  expand  under  certain  constraints.  These  motions  are 
controlled  by  a  single  parameter,  e,  which  can  be  varied  from  0 
to  1.  A  value  of  zero  represents  isotropic  behavior  and  a  value 
of  1  represents  kinematic  behavior.  The  resulting  yield 
surfaces  in  a  three-dimensional  principle  stress  space  are  shown 
in  Figure  2.  The  yield  surface  center  moves  at  a  rate  governed 
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Figure  2:  Hardening  Models  in  Principle  Stress  Space. 
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To  allow  for  finite  strains  and  rotations.  Updated 
Lagrangian  approach  is  adopted  [5].  The  coordinate  system  is 
convected  with  the  deformation.  In  this  coordinate  system,  the 
"true  strain"  rate  is  related  to  the  determination  rate  (or 
velocity)  through 
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In  the  absence  of  rotation  the  stress  tensor  is  related  to  the 
strain  rate  tensor  in  the  classical  manner,  i.e.. 


This  stress  measure  is  the  Cauchy  stress.  Under  finite 
rotations,  the  stress  tensor  is  not  invariant.  At  zero 
strain  rate,  the  stress  rate  is  given  by  where  W  are  the 
rotation  rates 
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The  total  stress-deformation  relation  is,  therefore. 
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(13) 


Equations  (1)  -  (13)  form  a  complete  incremental 
representation  of  finite  plastic  deformation.  It  only  remains 
therefore,  to  quantify  the  uniaxial  behavior  through  the 
function  f(ag).  There  are  many  functional  ways  to  represent 
uniaxial  loading  behavior.  From  a  computation  standpoint,  a 
multilinear  representation  is  easily  implemented  and,  by 
allowing  for  enough  segments,  can  be  arbitrarily  accurate. 
Consider  the  multilinear  representation  of  a  true  stress-true 
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strain  curve  shown  in  Figure  3.  The  functional  relationship 
between  the  stress  and  strain  are  given  by 
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The  plastic  strain  rate,  therefore,  is  given  by 


c  *  a  a  /E  (15) 
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Using  equation  (1)  and  recognizing  that  for  uniaxial 

deformation,  effective  quantities  are  proportional  to  the 

uniaxial  components,  the  function  f(o-)  can  be  reduced  to 

6 


The  function  is  only  linearly  dependent  on  the  current  slope  of 
the  uniaxial  curve.  By  specifying  enough  segments,  virtually 
any  hardening  behavior  can  be  accurately  described. 

Equations  (1),  (8)  and  (9)  provide  the  fundamental 
relationships  between  stress  and  strain  rates.  The  equilibrium 
conditions  (governing  equations)  for  a  continuum  body  in  the 
absence  of  body  forces  and  inertia  effects  can  be  written  as 


Multilinear  Approximation  For  A  Uniaxial  Stress 
Strain  Curve  With  Hardening. 
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with  the  boundary  conditions 

o  n  »  t  on  S 
ij  j  i  T 

and  (18) 

u  3  u  on  S 
i  i  u 

where  are  the  specified  loading  rates  on  the  boundary 

experiencing  applied  tractions  (S-p)  and  u^  are  the 

velocities  specified  on  the  remainder  of  the  boundary  (Sy). 
Equation  (13)  provides  the  fundamental  relation  between 

the  stress  state  and  the  deformation  gradients.  For  many 
problems  in  application  the  assumption  of  "small  strain" 
introduces  minimal  error  (mathematically,  this  means  assuming 
infinitesimal  displacements  and  strains).  If  this  assumption 
is  made,  the  strain  rates  are  related  to  the  velocity 
gradients  by 


This  simplification  also  means  that  the  reference  coordinate 
system  and  the  material  coordinate  system  are  coincidental 
throughout  the  deformation.  In  the  computer  code  described, 
the  option  of  finite  or  infinitesimal  strain  theory  is  left  to 
the  user.  The  use  finite  strains  slows  convergence 
considerably  for  problems  where  the  deformations  are  small. 

As  the  strains  grow,  however,  the  solutions  assuming 
infinitesimal  theory  diverge  from  the  finite  strain  results. 
Eventually  infinitesimal  solutions  will  fail  to  converge 
regardless  of  how  small  the  load  increments  are  tak.en. 

By  either  employing  the  Principle  of  Virtual  Work,  for 
increments  of  displacement  or  by  performing  the  standard 
Galerk-in  technique  on  the  governing  equations,  (17)  and  (18), 
the  finite  element  equations  governing  the  nodal  velocities, 

0,  can  be  written  in  terms  of  the  loading  rate  vector,  R,  in 
the  form 

K(U)  .  U  -  R  -  0  (20) 

The  standard  finite  element  assumptions  made  are  given  by 
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where  N  are  the  shape  functions.  The  set  of  rate 
equations  (20)  will  be  integrated  one  load  increment  (AR) 
at  a  given  time  to  determine  the  corresponding  new 
displacement  increment,  AU.  The  Newton-Raphson  or  tangent 
stiffness  solution  procedure  is  employed.  At  load  increment 

L  +  1,  the  initial  solution  AUi  is  found  from 
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K(U  )  .  AU^  >  AR  (22) 
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The  "new”  displacement  is  then  used  in  the  stiffness  matrix, 
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KfU  +  E  AU  ),  and  a  new  correction  is  obtained  from 
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where  the  integral  is  approximated  using  Simpson's  rule.  The 
procedure  is  repeated  until  two  convergence  criteria  are  met: 
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and  (24) 
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where  R  is  the  total  load  at  step  L  +  1. 
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In  this  study,  20-node  quadratic  isoparametric  elements 
were  employed  exclusively.  All  integration  was  carried  out 
utilizing  3x3x3  Gauss-Legendre  quadrature  formulae. 

Strains  were  calculated  at  the  Gauss  integration  points  in 
each  element  from  the  strain-displacement  relations  of  (19). 
Stresses  were  cumulatively  calculated  at  the  Gauss  points  from 
the  stress-strain  relations. 


Directly  calculating  strains  and  stresses  from  the  finite 
element  relations  (21)  at  points  on  element  boundaries 
inherently  yields  poor  results.  This  is  especially  true  when 
shape  functions  are  employed.  A  superior  approach  is  to 
calculate  the  stresses  and  strains  at  the  Legendre  quadrature 
points  and  to  extrapolate  or  smooth  them  to  the  boundaries. 
This  approach  has  been  shown  to  yield  very  accurate  results 
for  a  wide  variety  of  geometric  mappings.  In  this  study  the 
smoothing  technique  as  developed  in  [6]  is  employed  for  all 
stress  and  strain  evaluations. 

For  elastic-plastic  studies,  the  authors  prefer  to  model 
the  crack,  front  region  with  a  convergent  mesh  of  conventional 
elements  rather  than  to  employ  a  "singular"  element. 

Experience  with  both  elastic  and  elastic-plastic  studies  [3,7] 
demonstrates  this  approach  to  be  accurate  (although,  for 
elastic  problems,  more  costly).  Since  the  nature  of  the 
singularity  is  unknown  in  the  elastic-plastic  problem,  it  is 
presumptuous  to  employ  a  singular  element  and  may  lead  to 


erroneous  results. 


PROBLEM  DESCRIPTION 


Consider  a  panel  of  overaged  7075  (T7651)  aluminum  with  a 
central  through  the  thickness  crack..  A  typical  panel  Is  shown 
In  Figure  4.  The  panel  used  In  this  study  had  a  width  of 
8.89cm  and  a  crack  length  to  width  ratio  of  .0.5.  The  specimen 
thickness  was  0.984mm  and  the  specimen  length  was  17.78cm. 

The  uniaxial  stress-strain  curve  for  the  material  Is  shown  In 
Figure  5.  The  metallurgical  aspects  of  this  material  and  Its 
ductility  are  discussed  In  a  subsequent  section. 

Since  the  panel  was  loaded  normal  to  the  crack  only, 
symmetry  allowed  the  modeling  of  one  octant.  The  finite 
element  grid  used  In  this  study  Is  shown  in  Figures  6a,  6b  and 
6c.  The  smallest  elements  near  the  crack  front  had  planer 
dimensions  of  a/20  (where  a  Is  the  half  length  of  the  crack). 
The  convergence  of  this  grid  is  discussed  In  [3,7]. 

The  grid  shown  consists  of  96  20-node  isoparametric 
elements  with  624  total  nodes.  The  total  grid  has  1872 
degrees  of  freedom.  The  system  was  solved  by  the  frontal 
method.  The  total  storage  required  for  the  entire  program  was 
2.2  Megabytes  (for  double  precision  computations).  Total 
runtime  for  the  problem  discussed  was  48CPU  Hours  on  a  VAX 
11/780. 


Figure  S:  Uniaxial  Stress-Strain  Curve. 


EXPERIMENTAL  PROCEDURE 


The  ease  and  accuracy  of  measurement  of  the  crack,  tip 
plastic  zone  shape  will  be  high  if  the  plastic  zone  is  very 
large.  This  would  require  a  material  in  a  very  ductile 
condition.  However,  such  high  ductility  would  create  a  large 
curvature  in  the  crack,  front  during  fatigue  precracRing  and 
difficulties  in  obtaining  convergence  in  finite  element 
analysis.  Hence,  an  alloy  in  a  moderately  ductile  condition 
was  found  to  be  desirable.  These  conditions  were  obtained  in 
7075  aluminum  alloy  by  overaging  from  the  T651  condition  for 

72  hours  at  178*C  (352®F).  In  this  T7651  condition  the  alloy 

« 

had  a  yield  strength  of  307MPa  and  ultimate  tensile  strength 
of  407MPa. 

The  specimen  geometry  used  for  this  study  was  the 
center-cracRed  type  with  width,  w  *  89mm  and  cracR  length, 

2a  3  44.5mm.  Fatigue  cracRs  were  initiated  and  extended  from 
machined  notches  to  obtain  sharp  cracR  tips.  The  fatigue 
precracRing  was  performed  at  a  load  at  least  50%  lower  than 
the  load  applied  for  plastic  zone  formation.  After  the 
fatigue  cracR  was  grown,  the  specimen  was  loaded  to  a  desired 
load  value  to  produce  plastic  zones  at  cracR  tips.  The 
maximum  load  was  limited  by  the  load  needed  for  cracR  growth 
initiation,  as  the  experimental  results  were  to  be  compared 
with  the  results  of  finite  element  analysis  without  cracR 


extension.  Attempts  were  made  to  obtain  as  large  a  plastic 
zone  as  possible  without  crack,  growth;  hence,  the  selected 
load  was  very  close  to  that  needed  for  stable  crack,  growth 
initiation.  Although  initiation  of  stable  crack  growth  is 
generally  accompanied  by  a  sudden  drop  in  the  load,  in  ductile 
materials,  this  drop  is  not  easily  detectable.  Hence,  to 
assure  that  no  crack  growth  took  place  during  loading,  the 
specimen  was  fatigue  cracked  again  after  plastic  zone  size 
measurement  to  extend  the  crack  approximately  2.5  to  5mm.  The 
specimen  was  then  loaded  to  fracture.  If  crack  growth 
occurred  during  initial  loading  for  formation  of  the  plastic 
zone,  the  crack  growth  region  would  be  marked  by  a  dull 
appearance,  distinguishing  it  from  the  fatigue  crack  growth 
region  on  either  side.  The  results  from  such  specimens  were 
rejected. 

The  plastic  zone  size  was  determined  by  measuring  the 
permanent  reduction  in  thickness  after  the  specimen  was 
initially  loaded  and  unloaded.  The  contours  of  the  plastic 
zones  were  measured  using  a  surface  profile  measuring  device. 
The  sensor  of  the  device  consisted  of  a  pointer  with  a  small 
tip  radius  attached  to  one  end.  of  a  thin  hardened  titanium 
alloy  sheet  that  was  51mm  long.  19mm  wide  and  1.3mm  thick. 

The  other  end  of  the  sheet  was  rigidly  mounted  by  sandwiching 
between  two  aluminum  pieces.  Two  strain  gauges  of  resistance 
120ohms  were  mounted  on  each  face  of  the  titanium  alloy 
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sheet.  These  strain  gauges  formed  four  arms  of  a  Wheatstone 
bridge  circuit.  The  circuit  was  similar  to  those  used  in  load 
cells,  extensometers  and  clip  gauges.  The  output  from  the 
circuit  was  proportional  to  the  movement  of  the  pointer.  The 
signal  was  amplified  using  a  D.C.  conditioner.  The  use  of  the 
thin  titanium  alloy  sheet  reduced  the  pressure  on  the  specimen 
by  the  pointer,  and  scratching  of  the  surface  was  minimized. 

The  specimen  was  mounted  horizontally  on  a  table  with  two 
micrometer  screw  feeds  at  right  angles  to  each  other.  The 
specimen  was  mounted  in  such  a  way  that  the  direction  of  crack, 
(x-direction)  was  parallel  to  the  direction  of  traverse  of  one 
micrometer  screw.  When  the  sensor  was  mounted  the  pointer  was 
pressing  against  the  specimen  vertically.  The  specimen  was 
moved  underneath  the  sensor  using  the  micrometer  screw-feeds. 
The  output  from  the  sensor  was  used  to  drive  the  x-axis  of  an 
x-y  recorder.  The  y-direction  displacement  of  the  table  was 
measured  using  an  extensometer  attached  to  the  system,  and  the 
output  from  the  extensometer  was  used  to  drive  the  y-axis  of 
the  x-y  recorder.  Several  traverses  in  the  y-direction  were 
made  for  each  face  of  the  specimen  at  regular  intervals  of 
distance  from  the  crack  tip  in  the  x-direction.  The  curves 
obtained  from  these  traverses  were  used  to  establish  points 
around  the  crack  tip  corresponding  to  a  given  thickness 
reduction.  From  these,  contour  lines  for  different 
thicknesses  were  established.  This  process  was  repeated  for 
all  the  four  faces  of  the  specimen. 


EXPERIMENTAL  RESULTS 


The  contour  lines  delineate  the  size  and  shape  of  the 
plastic  zone.  The  thicRness  of  the  specimen  falls  between 
those  required  for  plane  stress  and  plane  strain  conditions. 
The  plastic  zones  obtained  had  a  shape  representative  of  this 
thickness  range. 

A  set  of  contour  lines  obtained  from  one  face  of  the 
specimen  are  shown  in  Fig.  7.  The  scatter  in  the  data  for  the 
outer  contour  lines  is  higher  than  the  inner  ones.  This  is 
because  the  rate  of  thickness  variation  decreases  with 
increasing  distance  from  the  cracked  tip,  as  can  be  seen  from 
the  differences  in  the  spacings  between  adjacent  contour 
lines.  Although  the  resolution  of  the  sensor  is  very  high  and 
is  limited  only  by  the  extent  of  the  amplification  of  the 
signal,  errors  can  be  introduced  due  to  any  nonplanarity  of 
the  initial  specimen  surface  and  slight  variations  in  the 
pressure  applied  on  the  micrometer  screws  while  advancing 
manually.  The  planarity  of  the  surface  was  checked  initially 
before  deforming  the  specimen.  The  variations-  in  pressure  can 
cause  an  error  of  approximately  0.0025mm.  The  depth  of 
0.0051mm  and  0.0102mm  represented  by  the  outer  contour  lines 
are  very  sensitive  to  these  variations.  The  scatter  in  the 
data,  also  is  produced  by  the  nonhomogeneity  of  the  material, 
caused  by  coring  and  inclusions  during  casting  and  orientation 
effects  during  subsequent  mechanical  processing. 


Crack  Tip 
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It  is  also  seen  that  the  plastic  zone  is  not  exactly 
symmetrical  to  the  initial  notch  direction,  which  also  can  be 
attributed  to  inhomogeneity.  Plastic  deformation  occurs  more 
extensively  in  the  softer  regions.  Material  inhomogeneity 
also  may  cause  change  in  the  orientation  of  the  fatigue  crack., 
which  tilts  the  zone  ahead  of  it.  Monsymmetry  can  also  be 
produced  by  misalignment  of  the  specimen  and  the  testing 
machine,  but  the  misalignment  in  the  set  up  used  was 
negligible. 

In  the  finite  element  analysis  inhomogeneity  is  not  tak.en 
into  consideration  and  hence,  the  zone  is  assumed  to  be 
symmetrical.  A  comparison  with  the  finite  element  results  can 
be  made  by  averaging  the  distance  of  each  set  of  contour  lines 
from  the  initial  notch  direction.  Such  contour  lines 
determined  from  the  four  faces  are  shown  in  Figs.  8-11.  The 
zone  sizes  are  slightly  different  for  the  four  faces.  This 
also  results  from  the  uneven  crack,  growth  during  fatigue 
cracking  due  to  inhomogeneity.  Since  the  excess  deformation 
in  one  region  is  compensated  by  the  lack,  of  it  in  another, 
averaging  the  results  minimizes  the  error  involved. 
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COMPARISON  EXPERIMENTAL  AND  FINITE  ELEMENT  RESULTS 


The  finite  element  analysis  discussed  previously  was 
performed  and  the  residual  deformation  after  unloading  the 
specimen  was  calculated.  These  results  were  plotted  as 
contraction  contours  and  are  shown  in  Fig.  12.  The  average 
results  (the  average  of  all  four  sides  as  previously 
discussed)  are  also  plotted  from  the  experiment.  The  results 
from  the  finite  element  analysis  are  in  good  agreement  with 
those  obtained  from  the  experiment. 

The  finite  element  predicts  slightly  more  plasticity  than 
the  experiment  predicts.  This  is  expected  since  the  finite 
element  formulation  assumes  all  nonlinearity  is  due  to  plastic 
(or  permanent)  deformation.  In  real  materials,  however,  there 
is  some  recoverable  nonlinear  deformation  (i.e.,  nonlinear 
elastic  deformation).  The  averaging,  material  inhomogeneity 
and  error  in  the  two  methods  more  than  account  for  the 
deviations  (less  than  4%,  maximum). 


OOSOSOmm 


ure  12:  Comparison  of  Surface  Displacements. 


CONCLUDIMG  REMARKS 


Comparison  was  made  between  the  experimental  and  finite 
element  surface  contractions  for  a  center-cracRed  panel.  The 
results  compare  favorably  indicating  the  accuracy  and 
realistic  modeling  of  this  finite  element  formulation  and 
modeling.  The  deviations  between  the  results  which  were 
observed  are  less  than  can  be  accounted  for  due  to  inherent 
error  in  the  measurements.  In  fact,  the  results  are  more 
accurate  than  one  would  expect. 

It  is  important  to  highlight  several  factors  when 
discussing  three  dimensional  finite  element  analyses  in 
general,  and  in  particular  for  nonlinear  problems.  The 
results  are  highly  dependent  on  the  grid  characteristics  and 
on  the  convergence  algorithm  employed.  Additional  degrees  of 
freedom  do  not  guarantee  a  more  accurate  solution  [1,2,3]. 
When  employing  three  dimensional  finite  element  models, 
convergence  studies  alone  are  not  sufficient.  Comparison 
between  predictions  and  true  material  behavior  is  essential. 
With  regard  to  fracture  problems,  it  is  essential  to  compare 
predictions  from  analysis  with  local  parameters  as  erroneous 
local  models  can  be  forced  to  produce  results  which  agree 
globally  (i.e.,  on  remote  quantities,  e.g.,  nonlinear 
compliance) . 


The  unique  aspect  of  the  study  presented  in  this  paper  is 


the  direct  comparison  with  experimentally  measured  local 
quantities  demonstrates  the  accuracy  of  the  modeling 
employed.  This  approach  can  now  confidently  be  applied  to 
fracture  problems  for  the  testing  of  fracture  criteria  and  the 
prediction  of  crack,  growth  and  instability.  Without  such  a 
demonstration,  numerical  solutions  and  verification  of  failure 
criteria  are  always  suspect.  Due  to  the  complicated  nature  of 
the  problem,  each  component  of  the  analysis  must  be  verified 
independently  to  guarantee  accurate  solutions  and  meaningful 
predictions.  This  component  is  lacking  in  three-dimensional 
studies  reported  to  date. 
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ABSTRACT 


The  modeling  of  fracture  specimens  and  the  application  to 
growth  problems  is  discussed  for  elasto-plastic  materials. 
Emphasis  is  placed  on  modeling  three-dimensional  effects  and 
slow  crack  growth.  Extensions  to  creep  modeling  is  discussed 
through  reference.  Directions  of  on-going  and  needed  future 
research  is  included. 

The  modeling  of  fracture  phenomena  in  materials  which 
exhibit  nonlinearities  has  been  a  topic  of  interest  for  many 
years.  As  early  as  the  1940*s,  researchers  have  recognized 
that  real  materials  used  in  engineering  practice  exhibit  both 
material  and  geometric  nonlinearities.  Since  that  time,  much 
effort  has  been  made  to  address  and  quantify  the  effect  of 
nonlinearities  on  fracture  and  ultimate  failure  of  materials. 
The  purpose  of  this  paper  is  to  summarize  many  aspects  of  the 
modeling  of  nonlinear  fracture  mechanics.  Specifically,  the 
paper  will  address  the  elasto-plastic  and  elasto-visco-plastic 
behavior  of  metallic  materials  and  the  applicability  of 
modeling  phenomena  to  fracture  problems  in  these  regimes.  The 
emphasis  will  be  on  numerical  modeling  schemes  which  can  be 
experimentally  verified  and  computationally  implemented. 
Failure  criterion  and  fracture  modeling  is  the  topic  of  a 
separate  study. 

The  problem  of  modeling  fracture  problems  in  nonlinear 
materials  can  be  subdivided  into  three  major  numerical  (or 


theoretical)  problems:  the  constitutive  modeling  of  the  funda¬ 
mental  nonlinearities  (both  geometrical  and  material),  the 
prediction  of  stable  crack,  growth  and  fracture  instability,  and 
the  numerical  modeling  of  stable  crack,  growth.  This  paper  will 
address  the  computational  approach  to  fracture  modeling  and  the 
three-dimensional  nature  of  ductile  crack,  specimens. 

Discussion  will  also  be  presented  on  the  modeling  of  stable 
crack,  growth  from  a  computational  viewpoint.  Failure  and 
fracture  prediction  can  only  be  intelligently  addressed  through 
studies  employing  accurate  modeling  of  the  nonlinear  problem 
and  are  in  the  realm  of  ongoing  and  incomplete  research.  The 
methodology  for  time  dependent  visco-plastic  behavior  will  be 
discussed  as  extensions  and  modifications  of  these  problems. 

COMPUTATIONAL  APPROACH  TO 
ELASTIC-PLASTIC  FRACTURE  PROBLEMS 

Many  authors  have  demonstrated  the  three-dimensional 
nature  of  ductile  fracture  phenomena  (for  a  review  of  the 
literature,  see  [1]).  To  address  the  accuracy  of  the  Finite 
Element  Method  and  standard  plasticity  analyses  for  analyzing 
fracture  specimens  it  is  necessary  to  examine  the'  local 
deformation  near  the  crack..  Several  authors  have  demonstrated 
that  quantitative  agreement  with  gross  specimen  quantities 
such  as  compliance,  global  energy  (including  "Energy 
Integrals"),  gauge  displacements,  etc.  do  not  guarantee 
accurate  local  predictions  (see,  for  example,  [2,3,4]).  To 


address  this  problem,  a  local  deformation  prediction  is 
compared  with  experiment. 

A  panel  of  overaged  7075  (T7651)  aluminum  was  tested 
experimentally.  The  specimen  was  a  center-cracked  panel  with 
width  8.89cm  and  crack  length  to  width  ratio  of  0.5.  The 
specimen  was  loaded  to  the  highest  applied  load  precluding 
crack  growth  and  subsequently  unloaded.  The  residual  surface 
deformations  were  measured  using  a  special  LVDT  probe  and 
measuring  technique  described  in  [5]. 

The  same  panel  was  analyzed  using  the  Finite  Element 
Method.  J-flow  theory  plasticity  was  employed  with  a 
multilinear  uniaxial  hardening  curve.  A  mixed  kinematic  and 
isotropic  hardening  law  with  balanced  weights  was  employed 
together  with  an  adaptive  load  incrementation  scheme  [6]. 
Geometric  nonlinearity  was  modeled  using  the  Updated-Lagrangian 
approach  [6,1].  The  converged  idealization  and  complete 
formulation  of  the  problem  can  be  found  in  [1].  The  panel  was 
loaded  incrementally  and  unloaded  to  zero  applied  load.  As 
expected,  much  residual  stress  and  deformation  remained.  The 
residual  displacements  on  the  specimen  surface  were  compared 
with  experiment.  The  comparative  results  are  summarized  in 
Figure  1.  The  solid  lines  are  the  numerically  predicted 
surface  constrictions  and  the  discrete  points  are  from  the 
experimental  data. 

The  average  of  the  four  sides  was  compared  as  the  local 
inhomogeneity  and  lack  of  symmetry  is  virtually  eliminated  in 


the  experimental  average.  The  results  demonstrate  the 
accuracy  of  the  finite  element  modeling  and  solution  procedure 
being  employed.  The  experimental  data  and  predicted  contours 
differ  by  less  than  3%.  This  is  more  accurate  than  was 
expected  comparing  with  claims  made  for  the  finite  element 
solution  of  crack,  problems  in  two  dimensions. 

The  comparison  presented  above  and  described  fully  in  [5] 
demonstrates  the  accuracy  of  the  current  numerical  approach 
and  modeling.  The  results  predicted  subsequently  would  model 
fairly  well  the  true  deformation  in  engineering  fracture 
specimens.  The  remainder  of  the  problems  to  be  discussed  will 
qualitatively  examine  some  of  the  more  important  aspects  of 
fracture  specimens  and  modeling. 

To  investigate  the  effect  of  material  hardening  model,  a 
center-crached  panel  8.89cm  wide  was  studied.  The  crack, 
length  to  width  ratio  is  0.5  and  the  specimen  length  to  width 
ratio  is  2.0.  The  idealization  is  the  same  as  discussed 
previously  [1].  The  loading  is  normal  to  the  crack,  direction 
and  reaches  a  maximum  load  equal  to  1/3  for  the  material  yield 
stress.  The  material  is  modeled  as  either  exhibiting 
kinematic  hardening  behavior,  isotropic  hardening  behavior  or 
a  mixed  hardening  behavior  (as  described  previously).  Since 
the  global  applied  loading  is  monotonic,  the  three  hardening 
models  would  predict  identical  response  if  the  local 
deformation  were  truly  proportional  in  nature  (assuming  no 
crack,  growth). 


Figures  2,  3  and  4  show  the  yield  zones  on  the  surface  of 
the  specimen  at  maximum  load  for  each  of  the  hardening  models. 


The  local  response  is  definitely  nonproportional.  The 
isotropic  model  predicts  more  yielding  on  the  crack,  extent  line 
than  either  of  the  other  two  models.  This  would  suggest  a  more 
ductile  response  ahead  of  the  crack,  (implying  a  greater 
tendency  toward  stable  crack,  growth  prior  to  final  failure). 

The  results  with  mixed  hardening  have  proven  to  be  the  closest 
to  what  is  observed  experimentally.  The  others,  therefore, 
should  be  viewed  with  that  fact  in  mind. 

The  largest  affect  of  hardening  model  is  the  yield 
characteristics  on  the  crack,  extent  line.  The  maximum  yield 
radius  and  the  ’’skewing"  of  the  yield  zones  is  fairly 
independent  of  the  hardening  model.  It  is  important  to 
emphasize,  however,  that  the  differences  between  the 
predictions  are  significant  and  the  local  response  is  highly 
nonproportional . 

Figures  5,  6  and  7  show  the  yield  zones  on  the  midplane  of 
the  specimen.  These  zones  demonstrate  the  same  hardening 
effects  as  do  the  surface  zones.  The  greatest  influence  of 
hardening  model  is  seen  on  the  line  of  crack  extent.  On  the 
midplane  of  the  specimen,  a  plane-strain  type  of  zone  would  be 
expected  (i.e.,  similar  to  zones  predicted  with  a  2-dimensional 
plane-strain  analysis).  The  mixed  hardening  model  demonstrates 
the  most  realistic  results  (which  is  consistent  with  the 
surface  observations). 


• 

^  ^  • 


A*, 


The  results  demonstrate  that  the  local  response  near  a 
crack,  is  not  of  the  proportional  type.  This  has  significant 
implications  with  respect  to  valid  failure  criteria  and 
analysis  models.  The  results  presented  also  demonstrate  that  a 
mixed  hardening  rule  is  the  most  realistic  for  modeling  the 
aluminum  alloys  investigated  so  far  and  is  probably  best  for 
most  engineering  metals.  A  complete  discussion  of  the 
hardening  modeling  effects  and  the  implications  of  these 
findings  can  be  found  in  [7]. 

To  study  thickness  effects,  the  specimen  thickness  is 
varied  and  the  different  zone  sizes  and  shapes  are  reported 
(the  same  center-cracked  panel  is  used  with  mixed  hardening 
assumptions).  Figures  8a,  8b,  8c  and  8d  show  the  surface  yield 
zones  as  a  function  of  thickness.  The  thickness  is  varied  from 
1.5  times  the  ASTM  plane  strain  requirement  to  a  very  thin 
panel  dimension  (total  thickness  of  about  3mm).  Even  for  the 
thinnest  specimens,  the  classical  ’’plane  stress”  zones  are  not 
recovered  demonstrating  that  the  assumptions  of  2-D  plane 
stress  are  not  valid  for  this  specimen.  Also,  even  for  very 
thick  specimens,  the  predicted  zone  does  not  conform  to 
classical  plane  strain  zones  (the  surface  deformation  will 
never  conform  to  2-D  assumptions  since  not  only  is  the  surface 
normal  stress  zero  but  the  strains  required  to  produce  plastic 
incompressibility  require  a  nonuniform  normal  strain  through 
the  thickness  which  is  incompatible  with  2-D  plane  strain).  It 
is  important  to  remember  that  a  state  of  stress  with  zero 


normal  and  antiplane  shear  does  not  need  to  reduce  to  2-D 
plane  stress. 

Figures  9a,  9b,  9c  and  9d  show  the  midplane  zones  for  the 
four  thicknesses.  As  the  thickness  increases,  it  is  obvious 
that  the  zones  approach  plane  strain  zones  (which  they 
should).  Even  at  the  largest  thickness,  however,  they  show 
some  skewing  due  to  the  finite  geometry.  From  these  results, 
the  ASTM  criterion  may  not  be  good  enough  for  ductile 
materials . 

The  yield  radii  and  extent  of  yielding  ahead  of  the  crack 
tip  are  summarized  in  the  tables.  It  is  important  to  note 
that  the  yield  radius  changes  by  about  10%  with  thickness, 
however,  the  yielding  ahead  of  the  tip  changes  drastically 
and,  thus,  the  plastic  area  changes.  Since  more  energy  is 
being  dissipated  with  larger  areas,  the  ductility  and  fracture 
properties  are  obviously  dependent. 

This  study  demonstrates  the  thickness  effects  on  local 
yield  characteristics  and  also  mandates  3-D  analysis  for 
accurate  quantitative  predictions.  A  complete  discussion  on 
the  effects  of  specimen  thickness  can  be  found  in  [6]. 

MODELING  STABLE  CRACK  GROWTH 

The  modeling  of  stable  crack  growth  is  important  for  the 
prediction  of  ductile  fracture  phenomena.  The  processes  of 
stable  crack  growth  and  plasticity  are  almost  always  present 
in  application  and,  indeed,  are 


Table  la  -  Yield  radii  as  a  function  of 
thick.ness  for  surface  yield  zones 


Thickness 


^max/a 


^o/a 


T 

T 

T 

T 


2 . 54cin 
1 .27cm 
6 .35mm 
3.175mm 

^max 


0.307 
0.327 
0.331 
0.343 

maximum  yield  radius 

yield  radius  along  crack,  line 


0.045 

0.075 

0.205 

0.296 


Table  lb  -  Yield  radii  as  a  function  of 
thickness  for  midplane  yield  zones 


Thickness 


^max/a 


'^o/a 


T 

T 

T 

T 


2 . 54cm 
1.27cm 
6 . 35mm 
3.175mm 


0.260 

0.280 

0.283 

0.299 


0.039 

0.071 

0.193 

0.288 


extremely  interrelated.  The  modeling  of  this  process, 
however,  must  be  numerically  accurate  independent  of  any 
particular  failure  criterion. 

Ductile  fracture  is  inherently  a  three-dimensional 
problem  as  has  been  demonstrated  both  in  this  paper  and 
elsewhere  ([1]  provides  a  review  of  the  topic).  No  studies  to 
date,  however,  have  successfully  modeled  a  slow  growth  process 
in  three-dimensions.  A  major  reason  for  this  is  the  extreme 
number  of  degrees  of  freedom  needed  to  model  such  a  process  if 
standard  crack  growth  (fixed  grid)  modeling  is  employed.  A 


better  approach  is  to  employ  the  "Mesh  Adaptive"  technique 
introduced  in  [8]. 

The  standard  method  used  to  simulate  crack,  growth  consists  of 
employing  a  grid  with  extremely  small  elements  along  the  line  of 
crack,  extent.  Nodal  forces  (or  displacements  through  a  spring 
release  method  £9])  are  relaxed  at  certain  load  levels  to  create  new 
free  surface.  If  a  criterion  independent  study  is  done,  this  method 
will  yield  very  good  agreement  with  experiment  (see,  for  example, 
[10],  for  a  typical  example).  Convergence  studies  have 
demonstrated,  however,  that  a  typical  Mode  I  center-cracked  panel 
which  exhibits  15%  total  crack,  growth  can  require  upwards  of  1500 
DOF  to  achieve  good  accuracy.  The  extension  of  this  method  to 
three-dimensions,  therefore,  is  not  realistic. 

In  the  "Mesh  Adaptive"  approach,  a  standard  stress  analysis  is 
performed  until  the  onset  of  stable  growth  is  predicted  (either 
directly  from  experiment  or  from  a  theory).  The  load  is  then 
incremented  a  small  amount  (the  amount  is  influenced  only  by 
numerical  convergence,  however,  an  amount  correspondent  to  2  percent 
crack,  growth  or  less  is  accurate  for  most  aluminum  alloys)  and  the 
amount  of  crack  growth  is  predicted  (or  calculated  from  experiment). 
The  near  crack  mesh  is  convected  to  the  new  location  of  the  crack 
tip.  The  stress  along  the  new  free  surface  is  relaxed  to  zero  and 
the  new  resultant  plastic  state  is  calculated  in  the  entire 
specimen.  The  new  stress  state  is  then  extrapolated  to  the  new 
geometry  and  the  process  is  repeated  for  each  increment  of  crack 
growth. 


approach  to  yield  accurate  predictions  for  several  2024 
aluminum  specimens.  Convergence  studies  on  a  previously 
studied  problem  [8]  have  demonstrated  a  DOF  reduction  of  about 
200%  and  a  computation  time  reduction  of  about  100%,  For 
three-dimensional  problems,  the  "Mesh  Adaptive"  method  should 
yield  even  greater  improvements.  Computation  time,  however, 
will  still  be  extreme  if  accurate  predictions  are  desired. 

MODELING  VISCO-PLASTIC  BEHAVIOR 

To  model  visco-plastic  or  creep  behavior  the  rate 
equations  of  stress  and  strain  are  formulated  to  satisfy 
equilibrium  and  geometric  compatibility.  The  formulation  (an 
outline  of  which  can  be  found  in  [11]  in  addition  to  many 
other  standard  texts)  then  requires  only  the  constitutive 
behavior  and  flow  rule.  From  the  standpoint  of  fracture 
mechanics,  the  material  modeling  is  the  only  unsolved 
problem.  The  fracture  modeling  is  done  exactly  as  in 
elastic-plastic  problems.  Fracture  prediction,  however,  is 
still  an  open  question. 

The  constitutive  modeling  of  visco-plastic  behavior  is 
still  an  open  issue.  Generally  the  models  can  be  divided  into 
two  categories:  those  that  depend  on  internal  variables  and 
those  which  utilize  only  stress  and  strain  measures  for  their 
flow  rules.  A  good  (although  somewhat  biased  review  can  be 


found  in  [11]).  From  a  computational  standpoint,  a 
visco-plastic  flow  rule  should  be  only  implicitly  dependent  on 
time  and  contain  parameters  which  are  true  material 
constants.  The  last  requirement  is  the  most  difficult  to 
achieve. 

An  additional  problem  involved  is  the  material 
variability  or  "scatter"  in  the  visco-plastic  regime.  Many  of 
the  materials  involved  show  uniaxial  scatter  of  100%  or 
greater  in  their  primary  creep  behavior  [12].  A  theory, 
however,  will  not  predict  this  scatter  if  it  is  deter¬ 
ministic.  Two  approaches  are  open:  choose  a  statistical 
theory  or  attempt  to  qualitatively  model  applications  by 
quantification  of  a  single  material  sample.  If  the  latter 
approach  is  chosen,  consistent  results  can  be  obtained.  This 
approach  is  the  one  currently  employed  by  most  authors  and 
reasonable  results  are  obtained  (for  a  review,  see  [13]). 

Studies  are  underway  to  address  these  problems,  however, 
visco-plastic  modeling  of  engineering  fracture  problems  is 
still  an  unsolved  and  open  question.  The  problem,  however,  is 
one  of  constitutive  modeling.  Fracture  modeling  is  identical 
to  that  for  elasto-plastic  problems. 


SUMMARY 


The  emphasis  of  this  paper  has  been  on  the  three-dimen¬ 
sional  nature  of  ductile  fracture  specimens,  a  short  review  of 
the  state-of-the-art  modeling  and  an  emphasis  on  the 
importance  of  careful  modeling  and  the  importance  of  con¬ 
vergence  studies  for  nonlinear  problems.  The  three-dimension¬ 
al  nature  of  fracture  must  be  addressed  before  meaningful 
prediction  can  be  made  on  the  failure  of  real  engineering 
materials  in  application.  Two-dimensional  studies  are  a 
necessary  first  step  in  the  research  process,  however,  great 
caution  needs  to  be  exercised  in  the  interpretation  of  re¬ 
sults  as  the  criteria  which  appear  to  predict  reasonable 
answers  to  a  few,  idealized  example  problems  may  not  have  any 
application  to  real  situations  (e.g.,  [14]). 


ACKNOWLEDGMENT 

This  work,  was  sponsored  under  the  Office  of  Naval 
Research  Contract  #N00014-84-K-0027  to  The  George  Washington 
University. 


13 


REFERENCES 


1.  Moyer,  E.  T.,  Jr.  and  Liebowitz,  H.,  Finite  Element 
Methodology  for  Elastic-Plastic  Fracture  Problems  in 
3-Dimensions.  To  appear  in  the  Commemorative  Issue  of  J. 
Numerical  Methods  in  Engineering  in  Memory  of  Bruce  Irons 
(1985). 

2.  Jolles,  M. ,  Schroedl,  M.  and  Watson,  T.  J. ,  On  the 
Geometry  Dependence  of  R-Curves,  21st  Annual  Mtg.  of  SES, 
Virginia  Polytechnic  Institute  and  State  University, 
Blacksburg,  Virginia,  October  15-17  (1984). 

3.  Kakreti,  A.  R. ,  Khan,  A.  S.  and  Iranmanesh,  A.,  "Fracture 
Prediction  in  Elastic-Plastic  Problems  Using  Finite 
Element  Method"  in  Applications  of  Fracture  Mechanics  to 
Materials  and  Structures  (eds.  6.  C.  Sih,  E.  Sommer  and  W. 
Dahl,  Nijhoff,  Hague)  (1983). 

4.  Bakker,  A.,  "On  the  Numerical  Evaluation  of  J  in 
3-Dimensions"  in  Applications  of  Fracture  Mechanics  to 
Materials  and  Structures  (eds.  G.  C.  Sih,  E.  Sommer  and  W. 
Dahl,  Nijhoff,  Hague) (1983) . 

5.  Moyer,  E.  T.,  Jr.,  Poulose,  P.  K.  and  Liebowitz,  H., 
Prediction  of  Plasticity  Characteristics  for 
Three-Dimensional  Fracture  Specimens  Comparison  with 
Experiment.  Accepted  for  publication  in  International  J. 
Engineering  Fracture  Mechanics,  October  (1985). 

6.  Moyer,  E.  Thomas,  Jr.  and  Liebowitz,  Harold,  Effect  of 
Specimen  Thickness  on  Crack  Front  Plasticity 
Characteristics  in  Three-Dimensions.  Proceedings  of  the 
6th  International  Conference  on  Fracture,  New  Delhi, 

India,  Dec.  4-10  (1984). 

7.  Moyer,  E.  Thomas,  Jr.  and  Liebowitz,  H.,  Plastic 
Deformation  and  Hardening  Characteristics  in 
Three-Dimensional  Fracture  Specimens.  Presented  at  ICF 
International  Symposium  on  Fracture  Mechanics  (Proceedings 
published  in  China  Science  Press),  Beijing,  China, 

Nov.  22-25  (1983). 


14 


8.  Moyer,  E.  Thomas,  Jr.  and  Liebowitz,  Harold,  A  Mesh 
Adaptive  Method  for  Modeling  Slow  Crack,  Growth.  Presented 
at  21st  Annual  Mtg.  of  SES,  Virginia  Polytechnic  Institute 
and  State  University,  Blacksburg,  Virginia,  October  15-17 
(1984) 

9.  Lee,  James  D.  and  Liebowitz,  Harold,  Considerations  of 
Crack.  Growth  and  Plasticity  in  Finite  Element  Analysis,  J. 
of  Computers  and  Structures,  No.  3/4,  pp.  401-403 
(1978) . 

10.  Gifford,  L.  N.  and  Hilton,  P.  D. ,  Preliminary 
Documentation  of  PAPST  -  Nonlinear  Fracture  and  Stress 
Analysis  by  Finite  Elements,  NDW-DTNSRDC  3960/43b,  Feb. 
(1981). 

11.  Chaboche,  J.  L.  and  Rousselier,  G. ,  On  the  Plastic  and 
Viscoplastic  Constitutive  Equations,  J.  Pressure  Vessel 
Tech.,  1^,  p.  153  (1983). 

12.  Garwood,  S.  J. ,  "Effect  of  Temperature  of  Orientation  and 
Constraint  on  the  Toughness  of  A533B  Class  I  Steel"  in 
Applications  of  Fracture  Mechanics  to  Materials  and 
Structures  (eds.  G.  C.  Sih,  E.  Sommer  and  W.  Dahl, 

Nijhoff,  Hague)(1983). 

13.  Wilson,  R.  E.  and  Palazotto,  A.  N. ,  Viscoplasticity  in  a 
Superalloy  at  Elevated  Temperatures  Considering  Tension 
and  Compression  Loading,  International  J.  Engineering 
Fracture  Mechanics,  22.,  p.  927  (1985). 

14.  Jolles,  M.  I.,  Schroedl,  M.  A.  and  Watson,  T.  J. , 
Three-Dimensional  Aspects  of  Fracture  Mechanics. 

Presented  at  22nd  Annual  Mtg.  of  SES,  Pennsylvania  State 
University,  University  Park.,  Pennsylvania,  October  6-9 
(1985). 


Comparison  of  surface  displacements 


.  2  -  von  Mises  stress  contours  near  the  crack  tip  on  specimen  surface  at 
imum  load  -  isotropic  hardening 


.  4  -  von  Mises  stress  contours  near  the  crack  tip  on  specimen  surface  at 
imum  load  -  mixed  hardening 


maximum  load  -  kinematic  hardening 


Fig.  8a  -  Surface  yield  zones  for  specimen  with 
2T  =  2.54cm 


Fig.  8b  -  Surface  yield  zones  for  specimen  with 
2T  *  1.27cm 


8d  -  Surface  yield  zones  for  specimen  wi 
3.175mm 


The  first  step  in  characterizing  creep  crack  growth  is  the 
establishment  of  valid  constitutive  parameters  for  modeling  the 
deformation  process.  To  this  end,  several  uniaxial  tests  were 
conducted.  These  results  are  then  processed  and  correlated 
with  postulated  constitutive  models. 

Testing  Apparatus 

An  "ATS"  Lever  Arm  Tester  Series  2410  was  used  for  the 
creep  testing.  The  test  stand  is  designed  to  apply  static 
stress  to  a  test  specimen  for  an  extended  period  of  time  at  a 
constant  elevated  temperature.  The  apparatus  consists  of  a 
balance  beam  that  connects  the  test  specimen  to  a  weight  pan. 
Ratios  of  3  to  1  and  20  to  1  are  designated  between  the  weight 
pan  and  the  specimen.  The  lower  ratio  is  used  to  provide 
optimum  accuracy  at  lower  loads.  The  weight  pan  is  part  of  the 
overall  weights  and  is  suspended  with  a  chain  to  prevent 
bending  movements  to  the  load  trains. 

On  the  specimen  side  of  the  machine,  a  balance  beam 
leveling  motor  compensates  for  the  elongation  of  the  test 
specimen.  If  this  is  not  available,  the  balance  beam  may 
become  unlevel,  thus- changing  the  calibration  of  the  weight 
system.  This  is  important  in  maintaining  the  load  accuracy 
within  the  ASTM  requirements. 

A  tabular  furnace  with  an  electrical-resistance  winding 


heats  the  test  specimen  through  radiation  in  an  air 
atmosphere.  The  furnace  has  three  heating  zones.  The  tube  is 


located  in  a  vertical  position,  with  pull  rods  connected  to  the 
specimen.  The  temperature  is  controlled  by  means  of  a 
thermocouple,  located  in  the  vicinity  of  the  test  specimen, 
within  ±1®C. 

When  creep  data  are  required,  the  specimen  strain  must  be 
measured  as  a  function  of  time.  This  is  difficult  because  use 
of  strain-measuring  transducers  is  not  practical  at  the  test 
temperature.  A  mechanical  linkage  must  be  attached  to  the 
specimen  to  transmit  the  strain  to  the  strain-measuring 
equipment  outside  the  high-temperature  environment.  A  linear 
variable  differential  transformer  was  used  to  measure  the 
strain.  It  consists  of  a  movable  metal  core  that  changes  the 
electrical  characteristics  with  small  motion  associated  with 
strain  measurements.  The  linkage  is  made  of  an  alloy  that  can 
withstand  the  test  temperature  encountered. 


Specimen  Preparation 

The  test  specimens  were  cut  from  a  1/2"  thick  sheet  of 
Alloy  718  with  composition  listed  in  Table  1.  The  tensile  axis 
of  the  specimen  were  perpendicular  to  the  rolling  direction. 

The  specimens  were  heat  treated  and  machined  to  the  ASTM 
specifications.  The  heat  treatment  procedure  is  shown  in  Table 
2.  The  specimens  were  then  surface  treated  in  order  to 
eliminate  any  scratches  and  other  stress  raisers  that  may  have 
remained  from  the  machining. 


Specimen  Loading 

Attention  was  paid  to  avoid  straining  the  specimens  when 
mounting  them  in  the  adapters  and  load  train.  With  the 
specimen  in  place,  the  load  train  was  examined  carefully  for 
any  misalignment  that  may  have  caused  bending  of  the  specimen 
under  load. 

The  upper  load  train  was  suspended  from  the  lever  arm,  and 
the  compensating  weight  adjusted  so  that  the  lever  arm 
balanced.  Strain-measuring  clamps  and  the  extensometer  were 
attached  to  the  specimen,  and  the  load  train  was  inserted  into 
the  furnace  with  the  specimen  centered.  The  specimen  was 
stabilized  at  the  temperature  before  loading  at  least  for  one 
hour . 

Loading  the  weight  pan  was  done  smoothly  and  without 
excessive  shock..  This  was  done  by  lowering  the  supporting  jack, 
under  the  load  pan.  In  the  case  of  step  loading,  the  weight 
was  placed  on  the  weight  pan  in  measured  increments,  and  the 
strain  corresponding  to  each  step  of  loading  was  recorded. 

Temperature  Control 

The  specimen  should  not  be  overheated  while  brought  to 
temperature.  To  do  this,  the  specimen  was  brought  to  10“C 
below  the  desired  temperature  in  about  1-1/2  hours.  Then,  over 
a  longer  period,  the  specimen  was  brought  to  the  desired 


temperature. 


Table  1.  Chemical  Composition  of  Alloy  718. 


Composition.  Weight  Percent 


0.04  C 

Bal.  Fe 

1.00  Ti 

0.19  Mn 

51.14  Ni 

0.52  A1 

0.05  Si 

17.96  Cr 

0.43  Co 

0.005  P 

3.12  Mo 

0.02  Cu 

0.005  S 

5.19  Cb  and  Ta 

0.003  B 

Table  2.  Heat  Treatment  of  Alloy  718. 

Anneal  at  955®C  for  1  hour  -  air  cool 

Age  at  720®C  for  8  hours,  furnace  cool  to  620*C 

Hold  at  620‘*C  for  total  aging  time  of  18  hours 

Results 

Two  series  of  uniaxial  tests  were  conducted:  one  at  550®C 
and  one  at  650®C.  The  results  are  summarized  in  Tables  3  and 
4.  The  strain  vs  time  curves  are  shown  in  Figures  1-9. 

Evident  from  these  tests  is  that  the  specimen  exhibits  moderate 
primary  creep  at  550®C  and  virtually  none  at  650®C.  This  is 
evident  at  all  the  load  levels  tested.  For  the  550®C  data,  at 
time  explicit  model  was  fit  to  the  data.  The  results  are  given 
in  Table  5.  These  results  can  either  be  used  directly  as  a 


constitutive  law,  or,  can  be  correlated  with  a  time-implicit 
differential  model  allowing  for  more  general  hardening.  Many 
candidate  models  exist  and  are  currently  being  investigated. 

The  650“C  data  exhibits  virtually  no  primary  creep. 
Unfortunately  (from  a  modeling  standpoint),  much  tertiary 
behavior  is  evident  prior  to  failure.  Since  tertiary  creep 
involves  many  deformation  and  failure  modes  simultaneously, 
modeling  is  difficult.  The  approach  to  be  used  in  fracture 
studies  is  to  assume  tertiary  behavior  does  not  occur  and  that 
the  fracture  mechanics  models  the  damage  which  is  not  accounted 
for  by  the  hardening  laws  of  primary  and  secondary  creep. 

Fracture  tests  are  currently  underway.  Compact  tension 
specimens  will  be  loaded  at  constant  load  and  temperature. 

Crack,  growth,  crack  opening,  gauge  displacement,  mouth 
displacement,  etc.,  will  be  measured.  Finite  element  studies 
are  also  underway  to  model  the  growth.  The  results  will  be 
presented  when  they  are  completed. 


Table  3.  Tensile  Properties  of  Alloy  718  at  550®C 

Yield  Strength  =  144.08  KSI 
Ultimate  Strength  =  176.78  KSI 
Elongation  =  20.72% 

Specimen 

Number  Stress  (KSI)  Time  to  Failure 


5 

152.9 

210 

hours 

4 

153.4 

87 

hours 

9 

155.1 

189 

hours 

6 

155.2 

212 

hours 

7 

157.2 

190 

hours 

Table  4. 

Tensile  Properties 

of  Alloy 

718  at  650'’C. 

Yield  Strength  = 

133.89  KSI 

Ultimate  Strength  = 

163.67  KSI 

Elongation  = 

25.7% 

Specimen 

Number 

Stress  (KSI) 

Time  to 

Failure 

13 

143.2 

90 

minutes 

14 

132.9 

240 

minutes 

15 

110.1 

41 

hours,  30  minutes 

16 

99.9 

227 

hours,  30  ninutes 

(Minute 


This  section  outlines  the  formulation  being  employed  in 
the  creep-fracture  studies.  The  total  strain  rate  is 
decomposed  into  an  elastic  and  visco-plastic  portion  in  the 
form  (uniaxially ) • 

•  •  • 

e=e+c  (1) 

e  vp 

where  the  subscript  e  will  refer  to  elastic  quantities  and  the 
subscript  vp  will  refer  to  plastic  quantities.  The  total 
stress  rate  is  related  to  the  elastic  strain  rate  in  the  form 

e  =  D  c  (2) 

e 

All  deformation  is  assumed  to  be  elastic  except  when  the 
flow  function,  F,  first  reaches  the  condition 

F(a,  e  )  -  F^  =  0  (3) 

vp  0 

where  Fq  is  the  initial  flow  stress.  The  total  uniaxial  stress 
can  be  written  as 

o  =  o  +  H'  e  (4) 

y  vp 

where  Oy  is  the  uniaxial  yield  stress.  The  elastic  (or 
initial  strain)  is  related  to  the  zero  time  stress  as 


This  relation  also  holds  at  other  times  relating  the  total 
stress  to  the  elastic  strain.  Integrating  the  above,  the 
visco-plastic  strain  is  given  in  the  form 

•  -rt 

e  =  c  t  +  c  (1-e  )  ( 

vp  0  A 


•  •  -rt 

€  =  e  +  e  re 

vp  0  A 


( 


where  c  ,  e  ,  r  are  material  parameters. 
0  A 


For  multi-dimensional  problems,  the  visco-plastic  strain 


rate  generalizes  to  the  form 


e  =  Y  «|)(F)>  ^  ( 

~vp  Bo 


where  y  is  the  fluidity  parameter  and  is  the  yield 


function.  The  total  strain  can  be  represented  as 


e  =  e  +  e  (I 

~e  ~vp 


and  the  visco-plastic  strain  can  be  simplified  to 


c 


=  Y  «I>>  a 


(1 


For  many  practical  applications,  the  choice  of  a  flow  function 
in  the  form 


<I'(F)  =  (F  -  Fn) 
TO 


N 


F  >  F 


(12) 


=  0  otherwise 

is  appropriate  (especially  for  metallic  materials) - 

The  equations  formulated  above  are  solved  incrementally 
using  an  explicit  time  integrator.  The  increments  of 


visco-plastic  strain  are  given  by  (at  time  t^^) 


n 

A  e 
~vp 


n 


At  e 
n  ~vp 


The  strain  rates  are  incremented  by  the  relation 


•  n+1 
c 

~vp 


•  n  n  n 
c  +  H  A  o 
~vp 


where 


n 


H 


,n  n  n 
a  c  \  =  H  (<J  ) 

d  a 


n  n  n  n 

A  a  =  D  A  c  =  D(A  c  -  b,  c  ) 

"  ~  ~e  ~  ~vp 


The  strain  increments  can  be  written  as  (in  terms  of  the 
incremental  displacement  additions) 


n  n  n 
A  e  =  B  Ad 


(13) 


(14) 


(15) 


(16) 


fl71 


yielding 


n  n  n  •  n 

A  cx  =  D  [B  d  d  -  e  A  t]  (18 

~  ~  ~  ~vp 

relation  for  the  stress  increment. 

The  B  matrix  relating  the  strain  to  the  displacements  can 
be  decomposed  into  a  linear  and  nonlinear  portion  as 


n  n 

B  =  B  +  B  (19 

~0  ~NL 

The  special  case  of  linear  deformations  renders  B  =0. 

~NL 

The  equilibrium  conditions  for  a  given  time  increment  can 
be  written  as 


n  T  n  n 

[B  ]  A  0  d  V  +  A  f 


(20 


where  A  f  are  the  incremental  forces  applied  during  the 
previous  time  increment.  The  solution  is 


n  n  -1  n 
A  d  =  [K^]  A  V 


(21 


where 


n  n  T  •  n  n 

AV  =rrB]  De  AtdV+Af 

and  the  stiffness  matrix  is 

n  n  T  n 

K  =  /  (B  )  D  B 


(22 


(23 


The  stresses,  displacements,  strain  increments  and  strains  are 
then  updated  as 


n+1  n  n 

o  =  a  +  A  a 


n+1  n  n 

d  =  d  +  A  d 


n  n  n  -1  n 
Ae  =BAd-D  Aa 
~vp  ~  ~ 

n+1  n  n 

c  =  c  +  A  e 
~vp  ~vp  ~vp 


For  most  problems  of  interest,  large  strains  can  occur. 
Total  Lagrangian  formulation  is  adopted  here.  The  B  matrix  is 
decomposed  as 


B  =  B  +  B 
-0  ~NL 


The  stiffness  matrix  can  then  be  written  as 


K  =  K  +  K 
~0  ~NL 


where  the  linear  and  nonlinear  portions  are  given  as 


K=JBDBdV 
~0  V  ~0  ~  ~0 


K  =  J  (B^D  B  +  B^  D  B  +  B^  D  B  )d  V 
~NL  V  ~0 - NL  ~NL - NL  ~NL - 0 


The  displacement  increments  can  then  be  written  as 


rAWtl 


(32 


Ad=K  AV+K  AV 
~0  ~  ~NL 


where 


n  n  n  n 

AV=AV  +  AV  +Af 
~0  -ML 


n  T  n 

AV=J‘BDe  AtdV 
~0  V  -0  -  ~vp 


n  T  •  n 

AV  =JB  De  AtdV 
ML  V  ML  vp  n 


The  total  stress  increment  becomes 


(33 

(34 

(35 


n  n  •  n  n 

A  o  =  D(B  A  d  -  e  A  t)  +  D  B  Ad  (36 

-  -0  ~  ~vp  -  -ML 

In  this  formulation,  the  strain  is  decomposed  into  a 
linear  and  nonlinear  component.  For  the  special  case  of  plane 
problems,  the  specific  forms  are 


e 


c 

-L 


+  e 
-ML 


(37 


The  B  matrices  are 


a  u  a  N 


a  V  a  N 


B 


NL 


Tic  Tir 

a  u  a  N 


Tic  d  X 

a  V  a  N 


TT  'S-y 
a  u  a  N 


a  u  a  M 

i _ i 

X  Tic  6  y 


avaN  +avaM 

_ i _ i 

Ty  d  X  Tic  d  y 


(39) 


The  strain  components  can  then  be  written  as 


e 

~L 


a  u 
Tic 


a  V 

s-y 


a  u 
s-y 


a  V 

Tic 


(40) 


NL 


2  2 
1  (a  u)  +  1  (3  V) 

7  7  Tic 


1  (|_Ji)^  >  1  (|_v)^ 
7  TT  7 


(3u)  (|_u)  +  (|_v.)  (|_v) 
L  Tic  Ty  Tic  Ty  . 


(41) 


The  formulation  summarized  above  has  been  implemented  in  a 
computer  code  and  rigorously  tested.  The  next  stage  of  the 
research  is  centering  on  modeling  slow  crack,  growth  and 


alternative  constitutive  formulations. 
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ABSTRACT 


The  governing  finite  element  system  for  elastic- 
plastic  analysis  of  fracture  specimens  in  three -dimen¬ 
sions  is  formulated.  The  formulation  accounts  for  mixed 
material  hardening,  finite  strains,  finite  rotations  and 
plastic  incompressibility.  The  implementation  of  these 
aspects  into  a  computational  formula  is  presented  and 
alternative  formulations  are  compared.  Small  strain  theory 
is  recovered  as  a  special  case  to  the  present  formulation. 

Analysis  is  performed  on  a  finite  thickness  center 
cracked  specimen.  The  grid  characteristics  required  for 
converged  solutions  are  discussed.  The  effect  of  material 
hardening  model  and  specimen  thickness  are  studied.  The 
local  yield  state  is  examined  as  a  gauge  of  the  local 
deformation  processes.  The  implications  on  the  fracture 
behavior  of  the  specimen  is  discussed. 

Local  surface  displacements  are  compared  to  experi¬ 
mentally  measured  yield  surfaces.  The  formulation  is 
shovsTi  to  predict  extremely  accurate  local  deformation  in 
the  neighborhool  of  the  crack  front.  Contrary  to  the  few 
three-dimensional  fracture  studies  carried  out  to  date, 
this  analysis  concentrates  on  the  local  deformation  be¬ 
havior  which  ultimately  controls  fracture.  Accurate 
resolution  of  this  behavior  is  essential  before  meaningful 
fracture  criteria  in  three-dimensions  can  be  developed. 


INTRODUCTION 


The  field  of  fracture  mechanics  has  developed, 
primarily  over  the  past  35  years,  as  an  important  and  useful 
tool  for  the  design  of  engineering  components  and  struc¬ 
tures.  Much  of  the  success  of  fracture  analysis  is  due  to 
the  advancement  of  the  finite  element  method  as  an  accurate 
approach  to  numerical  stress  analysis  for  problems  which 
involve  geometric  discontinuities  and  singularities. 
Originally,  fracture  studies  were  aimed  at  predicting 
brittle  failure  of  materials.  IVhile  a  necessary  first 
step,  this  approach  did  not  give  meaningful  results  to  many 
problems  of  practical  importance.  Much  of  the  later  work 
in  fracture  mechanics,  therefore,  is  geared  toward 
addressing  the  problem  of  ductile  fracture. 

Ductile  fracture  problems  involve  much  local  nonlinear 
deformation.  Analytical  approaches  to  the  concerned 
boundary  value  problem  have  proven  largely  futile.  To 
obtain  analytical  solutions,  far  too  many  unrealistic  and 
inappropriate  assumptions  are  required.  Numerical  solutions 
are,  therefore,  imperative.  Much  work  has  been  done  in  the 
area  of  elastic-plastic  stress  analysis  of  fracture  speci¬ 
mens.  Very  little,  however,  has  been  carried  out  in 
three-dimensions.  In  application,  most  components  are  not 
thick  enough  to  be  considered  in  plane  strain.  Indeed, 
most  true  plane  strain  fractures  are  brittle  in  nature. 

Plane  stress  analysis  on  fracture  specimens  is  rarely 


applicable.  The  assumption  of  plane  stress  violates  the 
compatibility  equations  which  are  essential  to  the  formula¬ 
tion  of  the  full  boundary  value  problem.  While  this  is 
unimportant  far  from  the  crack  tip,  elastic  solutions 
demonstrate  the  breakdown  of  the  assimption  near  the  crack 
tip  [1].  In  plastic  analysis,  the  assumption  is  further 
compounded  by  the  fact  that  four  general  deformation  regions 
can  be  identified  near  the  crack  front:  a)  the  near  tip, 
elastic,  plane  strain  region;  b)  a  near  tip  region  governed 
by  an  asymptotic  "plane  stress  like"  region  (which,  in  fact 
shows  contributions  from  both  elastic  plane  strain,  plastic 
plane  stress  and  possibly  boundary  influence) ;  c)  a  transi¬ 
tion  region  in  which  the  deformation  is  nonlinear,  not 
dominated  but  influenced  by  the  singular  stress  field  (also 
a  fully  three-dimensional  region);  d)  the  exterior  region 
which  is  either  plane  stress  dominated,  plane  strain 
dominated  or  fully  three-dimensional  depending  on  specimen 
geometry  and  thickness  [2,3]. 

While  two-dimensional  elastic-plastic  studies  may 
provide  some  insight  into  the  qualitative  aspects  of 
ductile  fracture,  they  cannot  be  expected  to  yield  quanti¬ 
tatively  useful  results  for  problems  which  are  not  dominantly 
brittle.  This  is  not  to  imply  that  two-dimensional  studies 
are  not  important.  They  must  be  done  to  provide  qualita¬ 
tive  and  preliminary  information  which  would  be  far  too 
costly  to  obtain  with  three-dimensional  analysis.  For  the 


3 


quantitative  prediction  of  ductile  fracture  phenomena,  how¬ 
ever,  two-dimensional  studies  cannot  be  relied  upon.  The 
success  reported  for  problems  with  reasonable  ductility  in 
the  literature  can  more  than  be  accounted  for  by  the  use  of 
global  fracture  criterion  which  are  theoretically  (as  well 
as  computationally)  unsound.  This  will  be  expounded  upon 
in  a  later  section. 

The  work  discussed  in  this  paper  has  evolved  over  the 
past  three  years  in  an  effort  to  systematically,  consistently 
and  accurately  address  the  problem  of  ductile  fracture.  The 
work  is  by  no  means  complete.  The  full,  incremental  elastic - 
plastic  finite  element  formulation  is  presented.  Much  debate 
exists  in  the  literature  with  regard  to  valid  formulations 
and  yield  criterion  for  ductile  fracture  specimens.  Specific 
choices  were  made  in  the  formulation  based  on  experience  with 
two-dimensional  studies.  It  is  believed  that  the  current 
approach  will  be  adequate  for  most  engineering  metals  at 
room  temperature. 

Results  are  presented  for  three  different  aspects  of 
ductile  fracture  specimens.  An  important  aspect  of  any 
numerical  study  is  the  validation  of  the  accuracy  of  the 
numerical  method  being  employed.  With  regard  to  finite 
element  modeling  of  ductile  fracture,  this  is  extremely 
important  (not  only  is  it  necessary  to  theoretically  formu¬ 
late  a  valid  approach,  it  is  critical  to  demonstrate  the 
accuracy  of  the  discretization  employed).  Verification  of 
the  finite  element  approach  is  demonstrated  in  this  work. 


A  comparison  is  made  with  experimentally  measured  local 
deformations.  This  approach  eliminates  the  problem  of 
criteria  bias  evident  in  most  ductile  fracture  studies. 
Results  are  also  presented  demonstrating  the  effect  of 
hardening  modeling  and  specimen  thickness  on  local 


plasticity  characteristics. 
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PLASTICITY  FORMULATION 


Several  approaches  to  continuum  plasticity  have  been 
proposed  over  the  past  several  decades.  Most  are  based  on 
either  stress  or  strain  yield  criterion  and  various 
hardening  assumptions.  No  clearly  "BEST"  approach  has 
emerged  to  date.  For  most  engineering  metals,  however, 
the  J2  flow  theory  of  plasticity  is  widely  accepted  for 
practical  applications.  This  approach  is  reviewed  and 
employed  in  the  current  investigations  [4]. 

The  J2  flow  theory  of  plasticity  assumes  that  the 
material  in  question  yield,  or  starts  deforming  plastically 
when  the  "effective  stress"  (or  von-Mises  stress)  reaches  a 
critical  value  (called  the  yield  stress) .  Prior  to  the 
onset  of  plasticity,  the  material  is  assumed  to  behave 
linear  elastically.  Subsequently,  the  deviatoric  stress 
components  are  related  to  the  deviatoric  strain  rate  through 
the  tensor  relation 


1  +  Vi  ^  c' 

— E—  Slj  *  7  ^  '“e’  Sij  c, 


■  ''yd 

A  >  0 

e 


1  +  V 


(otherwise) 


where  e^j  are  the  deviatoric  strain  rates  given  by 


®ij  '  ^ij  ■  T  ^pp  °ij 


e__  6s ^ 
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are  the  current  deviatoric  stress  components  given  by 


IJ  IJ  7  pp  IJ 


are  the  deviatoric  stress  components  measured  relative 
to  the  current  yield  surface  center  given  by 


S .  .  *  S .  .  -  a .  - 
11  ij  11 


a^j  are  the  coordinates  in  stress  space  of  the  current  yield 
surface  center,  are  the  Cauchy  stress  components,  e^^- 
are  the  "true"  strain  components  (discussed  in  a  subsequent 
section),  is  the  effective  stress  .given  by 


^  7  Slj 


and  a  is  the  effective  stress  measured  relative  to  the 


current  yield  surface  center 


’  /'I — ' - 

%  = 


The  function  f  (cjg)  is  derived  from  the  uniaxial  stress - 
strain  curve  and  is  consistent  with  the  Associated  Plas¬ 
ticity  theory  (a  complete  discussion  is  given  in 
Reference  [5],  Derivation  of  f  (o^)  for  a  multilinear 
representation  of  the  stress-strain  curve  will  be  discussed 
subsequently. 
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For  plastic  strains  which  are  incompressible,  the 
hydrostatic  plastic  strain  rate  is  zero.  The  total  hydro¬ 
static  strain  rate,  therefore,  is  related  to  the  hydrostatic 


stress  rate  by 

_  1  -  2v  1 

"PP  ■  -  ^pp 

Engineering  materials  exhibit  different  types  of 
uniaxial  hardening  behavior  when  subsequently  unloaded 
after  being  plastically  deformed.  Generally,  the  behavior 
falls  between  two  extremes  called  kinematic  and  isotropic 
hardening.  The  uniaxial  representation  of  these  behaviors 
for  a  bilinear  material  are  shown  in  Figure  1.  To  allow 
for  various  hardening  behaviors  in  the  multiaxial  formula¬ 
tion,  the  yield  surface  is  permitted  to  move  and  expand 
under  certain  constraints.  These  motions  are  controlled  by 
a  single  parameter,  6  ,  which  can  be  varied  from  0  to  1 .  A 
value  of  zero  represents  isotropic  behavior  and  a  value  of 
1  represents  kinematic  behavior.  The  resulting  yield  sur¬ 
faces  in  a  three-dimensional  principle  stress  space  are 
shown  in  Figure  2.  The  yield  surface  center  moves  at  a 
rate  governed  by 


I 


i 


Figure  2:  Hardening  Models  in  Principle  Stress  Space. 
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0  (otherwise) 
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6  *  0  -►  Isotropic  Hardening 
6  *  1  -►  Kinematic  Hardening 


The  radius  expands  according  to  the  relation 


0  =20  +26(0  -0) 

y  y  max  y' 


To  allow  for  finite  strains  and  rotations,  Updated 
Lagrangian  approach  is  adopted  [6).  The  coordinate  system 
is  convected  with  the  deformation.  In  this  coordinate 
system,  the  "true  strain"  rate  is  related  to  the  deformation 
rate  (or  velocity)  through 


'^ij  “  ^ijkl  ^kl 


In  the  absence  of  rotation  the  stress  tensor  is  related  to 
the  strain  rate  tensor  in  the  classical  manner,  i.e., 


1  8U. 

7  I -sTj  ■  Tir 


This  stress  measure  is  the  Cauchy  stress.  Under  finite 
rotations,  the  stress  tensor  is  not  invariant.  At  zero 


strain  rate,  the  stress  rate  is  given  by  where  W  are  the 
rotation  rates 


=  W.  a  . 
ip  PJ 


The  total  stress -deformation  relation  is,  therefore, 
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Equations  (1)  -  (13)  form  a  complete  incremental 
representation  of  finite  plastic  deformation.  It  only 
remains,  therefore,  to  quantify  the  uniaxial  behavior 
through  the  function  f(0g).  There  are  many  functional 
ways  to  represent  uniaxial  loading  behavior.  From  a 
computation  standpoint,  a  multilinear  representation  is 
easily  implemented  and,  by  allowing  for  enough  segments, 
can  be  arbitrarily  accurate.  Consider  the  multilinear 
representation  of  a  true  stress-true  strain  curve  shown  in 
Figure  3.  The  functional  relationship  between  the  stress 
and  strain  are  given  by 


e  =  w+  ^(o-,-o,)  + 


m 


^  (^3  -  02)  + . .  .p- 


(0  - 


The  plastic  strain  rate,  therefore,  is  given  by 


Multilinear  Approximation  For  A  Uniaxial  Stress 
Strain  Curve  With  Hardening. 


Using  equation  (1)  and  recognizing  that  for  uniaxial  de¬ 
formation,  effective  quantities  are  proportional  to  the 
uniaxial  components,  the  function  can  be  reduced  to 

fCo^)  - 


The  function  is  only  linearly  dependent  on  the  current 
slope  of  the  uniaxial  curve.  By  specifying  enough  segments, 
virtually  any  hardening  behavior  can  be  accurately  described 
The  above  formulation  is  chosen  for  implementation  vith 
the  Finite  Element  Method  to  be  described  subsequently.  It 
is  important  to  recognize  that  the  formulation  is 
representative  of  a  wide  class  of  materials  and  can  model 
many  different  variations  of  material  hardening.  This 
formulation  has  been  employed  with  considerable  success  in 
many  investigation  of  material  deformation  and  fracture 
mechanics.  Several  well  known  two-dimensional  codes  were 
developed  and  have  been  widely  employed  for  many  years  (see, 
for  example,  [7]  and  [8]). 


FINITE  ELEMENT  FORMULATION 


Equations  (1),  (8)  and  (9)  provide  the  fundamental 
relationships  between  stress  and  strain  rates.  The  equilib¬ 
rium  conditions  (governing  equations)  for  a  continuum  body 
in  the  absence  of  body  forces  and  inertia  effects  can  be 
written  as 

3o^j/3x^.  *  0  (1 

with  the  boundary  conditions 


where  t.  are  the  specified  loading  rates  on  the  boundary 
experiencing  applied  tractions  (S^)  and  u^  are  the  velocities 
specified  on  the  remainder  of  the  boundary  (S^) . 

Equation  (13)  provides  the  fundamental  relation  between 
the  stress  state  and  the  deformation  gradients.  For  many 
problems  in  application  the  assumption  of  "small  strain" 
introduces  minimal  error  (mathematically,  this  means 
assuming  infinitesimal  displacements  and  strains).  If  this 
assumption  is  made,  the  strain  rates  are  related  to  the 
velocity  gradients  by 
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Eij  »  ^Ou./aXj  3Uj/3x.)  (19) 

This  simplification  also  means  that  the  reference  coordinate 
system  and  the  material  coordinate  system  are  coincidental 
throughout  the  deformation.  In  the  computer  code  described, 
the  option  of  finite  or  infinitesimal  strain  theory  is  left 
to  the  user.  The  use  finite  strains  slows  convergence 
considerably  for  problems  where  the  deformations  are  small. 

As  the  strains  grow,  however,  the  solutions  assuming 
infinitesimal  theory  diverge  from  the  finite  strain  results. 
Eventually  infinitesimal  solutions  will  fail  to  converge 
regardless  of  how  small  the  load  increments  are  taken. 

By  either  employing  the  Principle  of  Virtual  Kork  for 
increments  of  displacement  or  by  performing  the  standard 
Galerkin  technique  on  the  governing  equations,  (17)  and  (18), 
the  finite  element  equations  governing  the  nodal  velocities, 

0,  can  be  written  in  terms  of  the  loading  rate  vector,  ft, 
in  the  form 

K(U)  •  ?  -  ft  =  0  (20) 

The  standard  finite  element  assumptions  made  are  given  by 


B‘D(U)  B  dV 


elements  element  volume 


where  N  are  the  shape  functions.  The  set  of  rate 
equations  (20)  will  be  integrated  one  load  increment 
(AR)  at  a  given  time  to  determine  the  corresponding  new 
displacement  increment,  AU.  The  Newton-Raphson  or  tangent 
stiffness  solution  procedure  is  employed.  At  load  increment 
L  +  1,  the  initial  solution  is  found  from 


The  "new”  displacement  is  then  used  in  the  stiffness  matrix, 


•  T.  and  a  new  correction  is  obtained  from 
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where  the  integral  is  approximated  using  Simpson's  rule.  Th 
procedure  is  repeated  until  two  convergence  criteria  are  met 
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'/  R 

tl+l  '  ~L-*-l 


where  R^+i  total  load  at  step  L  +  1. 

In  this  study,  20-node  quadratic  isoparametric  elements 
were  employed  exclusively.  All  integration  was  carried  out 
utilizing  3x3x3  Gauss-Legendre  quadrature  formulae. 
Strains  were  calculated  at  the  Gauss  integration  points  in 
each  element  from  the  strain-displacement  relations  of  (19)  . 
Stresses  were  cximulatively  calculated  at  the  Gauss  points 


from  the  stress  -  strain  relations. 

Directly  calculating  strains  and  stresses  from  the  finite 
element  relations  (21)  at  points  on  element  boundaries  in¬ 
herently  yields  poor  results.  This  is  especially  true  when 
shape  functions  are  employed.  A  superior  approach  is  to 
calculate  the  stresses  and  strains  at  the  Legendre  quadrature 
points  and  to  extrapolate  or  smooth  them  to  the  boundaries. 
This  approach  has  been  shown  to  yield  very  accurate  results 
for  a  wide  variety  of  geometric  mappings.  In  this  study  the 
smoothing  technique  as  developed  in  [9]  is  employed  for  all 
stress  and  strain  evaluations. 

Currently,  four  methods  of  accounting  for  the  crack  tip 
singularity  are  widely  employed.  Each  of  these  methods  is 
based  on  an  established  technique  in  LEFM  (Linear  Elastic 
Fracture  Mechanics).  The  first  method,  the  enriched  element 
approach  (where  the  shape  functions  are  modified  with  the 
asymptotic  crack  solution  vanishing  at  the  nodes)  has  been 
employed  both  for  the  multilinear  stress -strain  models  and 
for  power  law  hardening  models  [10].  Enriched  elements  based 
on  the  power  law  hardening  model  assume  that  the  enriched 
element  is  fully  yielded.  This  assumption  is  physically 
unrealistic,  especially  behind  the  crack  tip.  The  singular 
solution  employed  for  the  power  lav,'  hardening  case  also 
assumes  a  circular  yield  zone  which  is  far  from  realistic. 

The  solutions  generated  using  enriched  elements  and  a  multi¬ 
linear  stress  -  strain  assumption  are  reasonably  accurate  pro¬ 
viding  a  judicious  choice  of  enriched  element  size  and  sur- 


rounding  grid  characteristics  is  made.  The  major  drawback 
to  the  use  of  enriched  elements  is  the  computation  time 
required  to  obtain  convergence  due  to  element  incompatibility. 
The  second  method,  the  most  basic  approach,  uses  a  very  fine 
mesh  near  the  crack  tip  and  employs  only  conventional  elements. 
This  method  produces  reasonable  results  far  from  the  crack 
region  but  questionable  local  results.  Convergence  is  usually 
rapid,  therefore,  gross  specimen  behavior  can  be  obtained 
quickly.  V.'ith  unrealistically  fine  grids,  good  local  results 
can  be  obtained  (except  in  the  elements  bordering  the  crack 
tip)  but  only  at  the  expense  of  computer  time  [11] .  The 
third  method  is  based  on  the  fact  that  if  isoparametric 
elements  are  chosen  with  midside  nodes,  jucicious  choice, of 
the  placement  of  these  nodes  results  in  the  inducement  of  a 
/r  term  in  the  displacement  shape  functions  [12,  13].  These 
elements  are  essentially  equivalent  to  enriching  the  shape 
functions,  however,  element  compatibility  is  preserved  re¬ 
sulting  in  faster  convergence.  The  fourth  technique  of 
modeling  crack  tip  behavior  is  through  the  use  of  hybrid 
elements  where  elements  bordering  a  surface  with  traction 
boundary  conditions  are  forced  to  satisfy  those  conditions 
exactly  and  the  elements  bordering  a  surface  with  displac- 
ment  boundary  conditions  are  also  forced  exactly.  The 
element  boundaries  are  then  matched  by  using  Lagrange  multi¬ 
pliers  in  the  variational  equations  to  ensure  element 
equilibrium  and  continuity  in  an  approximate  sense.  Little 
work  has  been  done  on  comparisons  of  hybrid  methods  to 


conventional  methods  in  elastic-plastic  crack  problems, 
however,  the  technique  was  applied  with  questionable  success 
in  [14].  The  preferred  method  in  the  literature  is  still 
to  use  a  very  fine  mesh  and  standard  elements.  Complete 
discussions  of  the  above  methods  can  be  found  in  [15-17]. 


FINITE  ELEMENT  DISCRETIZATIONS 


In  this  paper,  three  basic  aspects  of  the  finite 
element  solution  of  ductile  fracture  problems  are  examined. 
All  three  involve  a  rectangular  panel  with  finite  thickness 
and  a  center  crack  through  the  thickness.  Due  to  the 
unknovTT  singular  nature  near  a  crack  front  and  the  problem 
of  the  intersection  of  a  crack  front  with  a  free  surface, 
only  conventional  (quadratic)  20-node  isoparametric  elements 
are  employed  in  this  work.  Studies  on  linear  elastic 
through  crack  specimens  has  demonstrated  the  accuracy  of 
this  approach  for  predicting  local  deformation  and  stress 
responses  [18]. 

For  all  the  problems  discussed  in  this  paper,  the  grid 
shown  in  Figures  4a,  4b  and  4c  is  employed.  This  grid  was 
used  to  predict  stress  intensity  factors  for  linear  elastic 
problems.  The  results  (discussed  in  [18])  demonstrate  the 
accuracy  is  on  the  same  scale  with  other  approaches.  The 
advantage  of  this  modeling  is,  since  the  singularity  is  not 
accounted  for  directly,  it  should  be  directly  employable  in 
an  elastic  plastic  study.  To  verify  this,  results  are 
compared  with  experimental  measurements  in  a  subsequent 
section. 

The  grid  shown  consists  of  96  isoparametric  elements 
with  624  total  nodes  (each  element  has  20  n«des).  The  total 
number  of  degrees  of  freedom  is  1872.  Runtimes  for  the 
.studies  to  be  presented  were  rather  lengthy  on  a  VAX  11/780. 


Typical  runtimes  for  panels  loaded  to  about  1/3  of  the 
material  yield  stress  converged  in  approximately  48  CPU 
hours.  This  is  fairly  consistent  with  existing  benchmarks 
of  computer  codes  which  employ  the  frontal  solution  method 
for  solution  of  the  stiffness  equations.  Obviously  it  is 
preferable  to  run  these  types  of  problems  on  larger  and 
faster  computers  if  they  are  available.  It  should  be 
appreciated,  however,  that  even  large  problems  can  be  run  on 
fairly  small  machines  using  the  frontal  approach. 

The  elastic  results  presented  in  [18]  demonstrate  the 
convergence  of  the  grid.  It  should  be  recognized  at  the 
outset  that  the  grid  employed  in  this  study  is  the  minimum 
required  for  convergence  to  the  elastic  solution.  The  next 
section  demonstrates  its  ability  to  model  elastic-plastic 
deformation  realistically.  While  it  is  never  possible  for 
nonlinear  problems  without  analytic  solutions  to  establish 
optimal  grid  requirements,  it  is  believed  that  little 
reduction  in  the  discretization  could  be  made  without 
seriously  compromising  the  accuracy  of  the  solution. 

Studies  which  employ  fewer  degrees  of  freedom  must  be 
considered  suspect. 


COMPARISON  WITH  EXPERIMENT 


A  panel  of  overaged  7075  CT7651)  aluminum  was  tested 
experimentally.  The  specimen  was  a  center-cracked  panel 
with  width  8.89cm  and  crack  length  to  width  ratio  of  0.5. 

The  stress-strain  curve  for  the  material  is  shown  in 
Figure  5.  The  specimen  was  loaded  to  the  highest  applied 
load  precluding  crack  growth  and  subsequently  unloaded. 

The  residual  surface  deformations  were  measured  using  a 
special  LVDT  probe  and  measuring  technique  described  in  [19] 

The  same  panel  was  analyzed  using  the  finite  element 
method  using  the  procedure  and  discretization  described 
previously.  The  panel  was  loaded  incrementally  and  unloaded 
to  a  zero  applied  load  state  (the  local  plastic  state,  how¬ 
ever,  showed  considerable  residual  deformation  and  stress  as 
would  be  expected).  Displacement  contours  near  the  inter¬ 
section  of  the  crack  with  the  free  surface  were  generated 
from  the  solution  and  compared  with  the  average  experimental 
results.  The  comparison  is  summarized  in  Figure  6.  The 
solid  lines  are  the  numerically  predicted  contours  and  the 
discrete  points  are  from  the  experimental  data. 

The  average  of  the  four  sides  was  compared  as  the  local 
inhomogeneity  and  lack  of  symmetry  is  virtually  eliminated 
in  the  experimental  average.  The  results  demonstrate  the 
accuracy  of  the  finite  element  modeling  and  solution  pro¬ 
cedure  being  employed.  The  experimental  data  and  predicted 
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ure  6:  Comparison  of  Surface  Displacements 


contours  differ  by  less  than  31^,  This  is  more  accurate  than 
vas  expected  comparing  with  claims  made  for  the  finite  ele¬ 
ment  solution  of  crack  problems  in  two  dimensions. 

The  comparison  presented  above  and  described  fully  in 
[19]  demonstrates  the  accuracy  of  the  current  numerical 
approach  and  modeling.  The  results  predicted  subsequently 
would  model  fairly  well  the  true  deformation  in  engineering 
fracture  specimens.  The  remainder  of  the  problems  to  be 
discussed  will  qualitatively  examine  some  of  the  more 
important  aspect  of  fracture  specimens  and  modeling. 
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MATERIAL  HARDENING  EFFECTS 


To  investigate  the  effect  of  material  hardening  model, 
a  center  cracked  panel  8.89cm  wide  was  studied.  The  crack 
length  to  width  ratio  is  0.5  and  the  specimen  length  to 
width  ratio  is  2.0.  The  idealization  is  the  same  as  dis¬ 
cussed  previously.  The  loading  is  normal  to  the  crack 
direction  and  reaches  a  maximinn  load  equal  to  1/3  for  the 
material  yield  stress.  The  material  is  modeled  as  either 
exhibiting  kinematic  hardening  behavior,  isotropic  hardening 
behavior  or  a  mixed  hardening  behavior  (as  described  pre¬ 
viously).  Since  the  global  applied  loading  is  monotonic, 
the  three  hardening  modeles  would  predict  identical  response 
if  the  local  deformation  were  truly  proportional  in  nature 
(assuming  no  crack  growth). 

Figures  7,  8  and  9  show  the  yield  zones  on  the  surface 
of  the  specimen  at  maximum  load  for  each  of  the  hardening 
models.  The  local  response,  is  definitely  non-proportional. 
The  isotropic  model  predicts  more  yielding  on  the  crack 
extent  line  than  either  or  the  other  two  models.  This  would 
suggest  a  more  ductile  response  ahead  of  the  crack  (implying 
a  greater  tendency  toward  stable  crack  growth  prior  to  final 
failure).  The  results  with  mixed  hardening  have  proven  to 
be  the  closest  to  what  is  observed  experimentally.  The 


others,  therefore,  should  be  viewed  with  that  fact  in  mind. 
The  largest  affect  of  hardening  model  is  the  yield 
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characteristics  on  the  crack  extent  line.  The  maximum 
yield  radius  and  the  "skeving"  of  the  yield  zones  is  fairly 
independent  of  the  hardening  model.  It  is  important  to 
emphasize,  however,  that  the  differences  between  the  pre¬ 
dictions  are  significant  and  the  local  response  is  highly 
non-proportional . 

Figures  10,  11  and  12  show  the  yield  zones  on  the  mid¬ 
plane  of  the  specimen.  These  zones  demonstrate  the  same 
hardening  effects  as  do  the  surface  zones.  The  greatest 
influence  of  hardening  model  is  seen  on  the  line  of  crack 
extent.  On  the  midplane  of  the  specimen,  a  plane -strain 
type  of  zone,  would  be  expected  (i.e.,  similar  to  zones  pre¬ 
dicted  with  a  2-dimensional  plane-strain  analysis).  The 
mixed  hardening  model  demonstrates  the  most  realistic  re¬ 
sults  (which  is  consistent  with  the  surface  observations) . 

The  results  demonstrate  that  the  local  response  near 
a  crack  is  not  of  the  proportional  type.  This  has  signifi¬ 
cant  implications  with  respect  to  valid  failure  criteria 
and  analysis  models  (which  will  be  discussed  in  the  con¬ 
cluding  section).  The  results  presented  also  demonstrate 
that  a  mixed  hardening  rule  is  the  most  realistic  for  modeling 
the  aluminum  alloys  investigated  so  far  and  is  probably  best 
for  most  engineering  metals.  A  complete  discussion  of  the 
hardening  modeling  effects  and  the  implications  of  these 
findings  can  be  found  in  [2], 
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SPECIMEN  THICKNESS  EFFECTS 


In  this  study,  the  specimen  thickness  is  varied  and  the 
different  zone  sizes  and  shapes  are  reported  (the  same  center 
cracked  panel  is  used  with  mixed  hardening  assumptions) . 
Figures  13a,  13b,  13c  and  13d  show  the  surface  yield  zones  as 
a  function  of  thickness.  The  thickness  is  varied  from  1.5 
times  the  ASTM  plane  strain  requirement  to  a  very  thin  panel 
dimensions  (total  thickness  of  about  3mm).  Even  for  the 
thinnest  specimens,  the  classical  "plane  stress"  zones  are 
not  recovered  demonstrating  that  the  assumptions  of  2-D 
plane  stress  are  not  valid  for  this  specimen.  Also,  even 
for  very  thick  specimens,  the  predicted  zone  does  not  con¬ 
form  to  classical  plane  strain  zones  (the  surface  deforma¬ 
tion  will  never  conform  to  2-D  assumptions  since  not  only  is 
the  surface  normal  stress  zero  but  the  strains  required  to 
produce  plastic  incompressibility  require  a  nonuniform  nor¬ 
mal  strain  through  the  thickness  which  is  incompatible  with 
2-D  plane  strain).  It  is  important  to  remember  that  a  state 
of  stress  with  zero  normal  and  anti-plane  shear  does  not 
need  to  reduce  to  2-D  plane  stress. 

Figures  14a,  14b,  14c  and  14d  show  the  midplane  zones 
for  the  four  thicknesses.  As  the  thickness  increases,  it 
is  obvious  that  the  zones  approach  plane  strain  zones  (which 
they  should).  Even  at  the  largest  thickness,  however,  they 
show  some  skewing  due  to  the  finite  geometry.  From  these 
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results,  the  ASTM  criterion  may  not  be  good  enough  for 
ductile  materials. 

The  yield  radii  and  extent  of  yielding  ahead  of  the 
crack  tip  are  summarized  in  the  tables.  It  is  important 
to  note  that  the  yield  radius  changes  by  about  10%  with 
thickness,  however,  the  yielding  ahead  of  the  tip  changes 
drastically  and,  thus,  the  plastic  area  changes.  Since 
more  energy  is  being  dissipated  with  larger  areas,  the 
ductility  and  fracture  properties  are  obviously  dependent. 

This  study  demonstrates  the  thickness  effects  on  local 
yield  characteristics  and  also  mandates  3-D  analysis  for 
accurate  quantitative  predictions.  A  complete  discussion 
on  the  effects  of  specimen  thickness  can  be  found  in  [3]. 
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Thickness 


r  , 
max/< 


r 


o/a 


T  » 

2.54  cm 

0.307 

0.045 

T  - 

1.27  cm 

0.327 

0.075 

T  ■ 

6.35  mm 

0.331 

0.205 

T  - 

3.175  mm 

0.343 

0.296 

r  maximum  yield 

max 

radius 

-*■  yield  radius  along  crack  line 


Table  la:  Yield  Radii  As  A  Function  Of  Thickness  For 
Surface  Yield  Zones. 


Thickness 

*^max/a 

^o/a 

T  =  2.54  cm 

0.260 

0.039 

T  =  1.27  cm 

0.280 

0.071 

T  »  6.35  mm 

0.283 

0.193 

T  “  3.175  trm 

0.299 

0.288 

Table  lb  :  Yield  Radii  As  A  Function  Of  Thickness  For 
Midplane  Yield  Zones. 
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CONCLUSIONS 


This  paper  presents  a  unified  and  complete  computational 
approach  to  solving  elastic-plastic  engineering  problems. 

The  application  of  this  approach  to  fracture  specimens  is 
presented  including  the  mesh  discretization  necessary  to 
produce  accurate  results.  Several  of  the  more  important 
aspects  of  cracked  specimens  (i.e.,  material  modeling, 
specimen  thickness,  accuracy  requirements)  are  examined. 

The  results  of  the  experimental  study  demonstrate  that 
the  current  approach  and  discretization  are  extremely  accurate 
for  predicting  detailed  local  deformations.  This  comparison 
is  essential  for  establishing  the  validity  of  the  computa¬ 
tional  procedure  and  discretization  employed.  Comparison 
of  global  parameters  far  from  the  crack  (such  as  compliance, 
gauge  displacement,  mouth  opening,  etc.)  are  not  sufficient 
to  establish  a  given  approach  or  modeling  as  accurate.  The 
results  of  this  study  also  demonstrate  that  the  mixed 
hardening  model  is  the  most  accurate  for  the  material  studied. 

Initially,  the  effect  of  hardening  parameter  vas 
examined  to  delineate  the  differences  in  local  deformation 
during  unloading.  It  was  discovered,  however,  that  significant 
local  response  differences  are  predicted  during  the  loading 
segment  also.  This  is  in  direct  opposition  to  the  assump¬ 
tion  of  proportional  loading  often  made  in  the  literature. 
Indeed,  the  assumption  is  imperative  for  "deformation" 
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theory  plasticity  studies  to  have  any  validity.  The  results 
presented  (already  established  through  experimental  compari¬ 
son)  clearly  demonstrate  that  the  local  response  is  dis¬ 
tinctly  non-proportional.  This  result  also  demonstrated 
that  local  fracture  and  failure  criteria  are  the  only  criteria 
which  are  theoretically  plausible  (as  proportional  loading  is 
essential  for  the  validity  of  global  criteria)  are  local 
criteria.  This  has  been  suggested  by  many  authors  in  the 
past  and  has  often  been  a  point  of  argument  in  the  fracture 
community.  Global  criteria  such  as  J-integral,  nonlinear 
compliance,  crack  mouth  opening,  etc.,  all  require  propor¬ 
tional  local  loading.  These  criteria,  therefore,  cannot 
validly  be  employed  for  nonlinear  fracture  problems. 

The  effect  of  specimen  thickness  was  investigated  to 
determine  the  range  of  validity  for  the  assumptions  of  plane 
stress  and  plane  strain  near  a  crack  tip.  The  results 
presented  demonstrate  that  even  for  extremely  thin  specimens  , 
the  local  response  does  not  behave  in  a  plane  stress  manner. 
These  results  also  demonstrate  that  the  local  deformation 
does  not  approach  plane  strain  behavior  until  the  thickness 
is  far  beyond  the  ASTM  requirement.  ^Tiile  two  dimensional, 
studies  need  to  be  performed  as  preliminary  qualitative  in¬ 
vestigations,  they  cannot  be  expected  to  produce  results 
which  predict  local  deformation  quantitatively. 


The  extensions  of  the  work  presented  needed  to  fully 
unravel  the  mysteries  of  fracture  mechanics  are  endless.  Tc 


begin,  however,  several  directions  are  clear.  Computational! 
it  is  important  to  extend  the  algorithms  to  account  for  crack 
growth  in  a  consistent  and  accurate  manner.  The  procedure, 
however,  must  be  independent  of  any  particular  failure 
criterion  and  not  require  excessive  computational  require¬ 
ments.  Recent  work  has  proposed  one  approach  to  the  problem 
and  demonstrated  its  efficiency  and  accuracy  in  two-dimen¬ 
sions.  Further  study  is  needed  for  the  extension  to  three - 
dimensions . 

The  major  theoretical  area  in  which  more  research  is 
needed  is  in  the  area  of  failure  criteria.  It  is  beyond  the 
scope  of  this  work  to  review  all  the  criteria  which  have 
been  proposed  to  address  nonlinear  fracture  problems.  The 
results  of  the  work  presented  in  this  report,  however, 
demonstrate  that  the  global  criteria  presented  must  be 
abandoned  as  they  rely  on  the  assumption  of  propositional 
loading  which  is  incorrect.  The  local  criteria  presented 
to  date  have  not  been  sufficiently  tested  to  establish  any 
one  as  a  truly  valid  fracture  criteria.  The  most  widely 
tested  have  all  failed  the  critical  requirements  of  geometry 
and  specimen  independence  and  consistency  over  a  vide  range 
of  applications.  Perhaps  a  new  approach  is  needed.  It  is 
important  to  continue  pursuing  local  fracture  criteria. 

The  ability  to  do  large  deformation,  three-dimensional 
analysis  without  the  a  priori  bias  of  a  failure  criterion 
should  help  in  the  development  of  a  consistent  and  accurate 
approach  to  elastic-plastic  fracture  problems. 
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ABSTRACT 

A  finite  element  investigation  of  the  effect  of  thickness  on  plastic  defor¬ 
mation  and  yielding  characteristics  in  three-dimensional  cracked  bodies  is 
presented.  It  is  shovn  that  the  fundamental  deformation  modes  and  extent  of 
plastic  deformation  are  significantly  influenced  by  the  specimen  thickness. 
The  results  show  the  transition  from  a  local  plane  strain  to  plane  stress 
response  near  the  crack  front  as  the  specimen  thickness  is  decreased.  While 
the  results  are  generated  for  a  specific  aluminum  alloy  (7075-T7651) ,  the 
predictions  for  other  hardening  materials  would  be  ouailtatively  the  same. 

KEYICORD-S 

Nonlinear  finite-element  calculations,  plastic  deformation,  three-dimension¬ 
al  crack  specimens,  incremental  analysis, 

INTRODUCTION 

Of  fundamental  importance  to  the  accurate  fracture  assessment  of  components 
and  structures  made  of  metals  is  the  study  of  ductile  fracture  processes  and 
the  plastic  response  near  a  crack.  The  basic  deformation  response  near  the 
crack  front  must  be  resolved  accurately  for  reliable  predictions.  Fracture 
criteria  have  been  proposed  based  on  many  controlling  quantities  (e.g., 
stress,  strain,  energy,  displacements,  etc.)  both  on  global  and  local  scale 
levels.  Without  exception,  all  of  these  criteria  require  accurate  local 
deformation  modeling. 

To  undf.rstand  the  scale  shifting  effects  from  the  laboratory  specimen  to  the 
structural  component,  it  is  imperative  to  discover  the  effects  of  specimen 
thickness  on  the  deformation  response.  This  problem  is  an  essentially 
three-dimensional  one  and  must  be  investigated  accordingly. 

The  purpose  of  this  Investigation  is  to  delineate  the  effect  of  specimen 
thickness  on  local  crack  front  yielding  characteristics  in  a  cracked  speci¬ 
men.  The  three-dimensional  elastic  plastic  finite  element  code  developed  in 
[1]  is  employed  for  the  analysis.  Specimen  thicknesses  investigated  range 


from  well  beyond  ASTM  plane-strain  requirements  to  thin  sheet  dimensions. 

The  yield  rones  calculated  in  this  work  demonstrate  the  transition  from  di- 
latational  to  distortional  dominance  ahead  of  the  crack  tip  as  a  function  of 
thickness  (equivalent  to  a  transition  from  plane  strain  to  plane  stress). 

The  magnitude  as  well  as  the  extent  of  yielding  is  shown  to  be  highly  thick¬ 
ness  dependent.  The  results  of  this  study  also  demonstrate  that  two-dimen¬ 
sional  analysis  based  on  plane  strain  (for  thick  specimens)  or  plane  stress 
(for  thin  specimens)  can  fall  to  accurately  model  the  local  response  when 
simple  standards  would  dictate  otherwise. 

PLASTICITY  FOBMULATION 

The  Incremental  theory  of  plasticity  employed  In  this  work  Is  based  on  the 
classical  rate  proportionality  assumptions  and  J-  flow  theory.  While  the 
mathematical  details  vary  with  the  choice  of  yield  criteria,  the  salient 
features  of  all  Incremental  theories  are  the  same.  This  discussion  will, 
therefore,  be  confined  to  the  specific  theory  employed  In  this  work. 

Assuming  stress  strain  rate  proportionality  and  J-  flow  theory  (which 
assumes  the  plastic  deformations  are  Incompressible)  the  stress-strain  rate 
relations  can  be  written  as  [2] 


4^  4  f  (Oe) 


c  •  a  ;  a  >0 
e  ye 


1  +  V 


Otherwise 


where : 


are  the  devlatorlc  stress  components. 


e. .  *  -  -r  £  _  6..  are  the  devlatorlc  strain  rate  components, 

ij  ij  3  pp  Ij 

V  Is  Poisson's  ratio, 

E  Is  Young's  modulus, 

S.  .  ■a.,  -  IT  a  _  i5..  are  the  devlatorlc  stress  components, 

Ij  Ij  J  PP  ij 

are  the  coordinates  in  stress  space'  of  the  yield  surface  center 

s'  *  S. .  -  a.,  are  the  devlatorlc  stress  components  measured  relative 
^  ^ ^  to  the  current  yield  center, 

Y  ^Ij^lJ  effective  stress. 


/  Y  ^ii^ll  effective  stress  relative  to  the  current  yield 

^  ^  center. 

Is  the  current  yield  stress,  and 

*  denotes  time  differentiation. 

Due  to  the  incompressibility  condition,  the  hydrostatic  strain  rate  Is  pro¬ 
portional  to  the  mean  stress  rate  and  Is  given  by 

•  1  -  2v  1 


The  function  f(<T^)  is  dependent  on  the  uniaxial  stress-strain  curve  and  will 
be  discussed  subsequently.  For  a  von  Mises  (J.)  material,  the  center  of  the 
yield  surface  moves  at  a  rate  proportional  to  the  projection  of  the  stress 
rate  vector  onto  the  local  normal  to  the  current  yield  surface  and  can  be 
written  as 


(3) 
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c  >  0 
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■ 

0  Otherwise 

where  6  varying  from  0  to  1  will  model  hardening  behavior  from  kinematic 
(B  “  0)  to  isotropic  (B  ■  1). 


The  function  f(<7  }  is  derived  from  Che  uniaxial  stress-strain  curve.  For  a 
uniaxial  speclmeS,  ec^uacion  Cl)  reduces  to 
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in  Che  plastic  range.  Thus, 
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f(a  )  can  be  written  as 
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If  the  uniaxial  stress-strain  curve  is  expressed  in  a  multilinear  fashion 
Che  stress-strain  relation  is 

where  a  ,  <  o  <  a  and  a  is  given  by 
m-1  m  m  “  ^ 

EAe  -  Aa 
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From  equation  (7) ,  the 


a  a 
m  e 


plastic  E 
and  thus  from  (6) 


“  eT- 
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plastic  strain  rate  is  given  by 


(9) 


CIO) 


Equations  CD  t  (2) ,  C3)  and  CIO)  provide  a  complete  set  of  elastic-plastic 
constitutive  relations.  Together  with  the  equilibrium  equations  and  Che 
strain-displacement  relations,  a  governing  system  will  be  formed.  It  is 
important  to  note  that  the  constitutive  formulation  outlined  above  is 
acceptable  for  finite  aa  well  as  infinitesimal  sti’ains.  Also  of  Importance 
is  the  fact  that  this  formulation  is  strain-rate  independent.  This  assump¬ 
tion  appears  to  be  realistic  for  most  engineering  metals  at  room  temperature 
Cor  cooler).  For  high  temperature  problems  a  rate-independent  formulation 
is  dubious. 


Equations  Cl) .  C2) ,  (3)  and  CIO)  provide  the  fundamental  relationships  be¬ 
tween  stress  and  strain  rates.  The  equilibrium  conditions  Cgoveming  equa¬ 
tions)  for  a  continuum  body  in  the  absence  of  body  forces  and  inertia 
effects  can  be  written  as 


with  the  boundary  conditions 
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where  T  are  the  specified  loading  races  on  the  boundary  experiencing  applied 

X 

tractions  CS„)  and  u  are  the  velocities  specified  on  the  remainder  of  the 
boundary  (S  Utilizing  Che  standard  infinitesimal  strain-displacement 
relations 
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(13) 


and  either  employing  Che  Principle  of  Virtual  Work  for  Increments  of  dis¬ 
placement  or  by  performing  Che  standard  Galerkln  technique  on  the  governing 
equations,  (^1)  and  (12),  Che  finite  element  equations  governing  the  nodal 
velocities,  U  can  be  written  In  terms  of  the  loading  rate  vector,  R,  in 
the  form 


K(U)  •  U  -  R  -  0 

The  standard  finite  element  assumptions  made  are  given  by 


(14) 
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where  N  are  Che  shape  functions.  The  set  of  rate  equations  (14)  will  be  In- 
tegracid  one  load  increment  (AR)  at  a  given  time  to  determine  the  corre¬ 
sponding  new  displacement  Increment,  AU.  The  Newcon-Raphson  or  tangent 
stiffness  solution  procedure  Is  employed  as  described  In  [3J. 


PROBLEM  DESCRIPTION 


To  study  the  effects  of  specimen  thickness  on  the  yielding  characteristics 
of  typical  fracture  specimens,  a  finite  center-cracked  plate  was  chosen  for 
Investigation.  The  standard  mode  I  configuration  shown  In  Pig.  1  was 
analyzed  for  total  thicknesses  of 

2T  •  2.54  cm  2T  -  1.27  cm 

2T  -  6.35  mm  2T  -  3.175  mm 

The  material  Investigated  was  a  7075-T7651  alumlntim  alloy  with  elastic  pro¬ 
perties 

E  «  7.24  E+04  MPa 
V  -  0.3 

Oy  -  4.07  E+02  MPa 

The  uniaxial  streaa-straln  curve  Is  shown  In  Fig.  2. 


z/ 


w  «  S.89  CB 


Fig.  1.  Through  crack  geometry  and  loading. 


Fig.  2.  Uniaxial  stress-strain  curve  for  7075-16731  ftluminum. 

The  finite  element  discretization  employed  in  the  analysis  utilizes  20-Node 
quadratic  isoparametric  elements  exclusively.  A  fine  mesh  near  the  crack 
front  is  employed  for  accurate  modeling.  The  grid  characteristics  and  con¬ 
vergence  properties  are  discussed  in  [1,4].  The  maximum  load  applied  was 

a  -  1.77  £+02  MPa 
max 

A  hardening  parameter  of  0  ■  0.5  was  also  assumed  in  the  analysis. 

RESULTS  AND  DISCUSSION 

The  yield  zones  predicted  at  the  maximum  load  for  each  of  the  four  thick¬ 
nesses  studied  were  calculated  and  plotted  both  on  the  surface  and  midplane 
of  the  specimen.  The  results  demonstrate  the  significant  influence  thick¬ 
ness  has  both  on  the  nature  and  extent  of  the  yielding. 


Figure  3a  is  a  plot  of  the  von  Mises  stress  contour  corresponding  to  the 
specimen  yield  stress  calculated  at  the  maximum  load  on  the  surface  of  the 
2.54  cm  thick  specimen.  As  expected  for  a  thick  specimen,  this  zone  has  the 
characteristic  form  of  a  plane  strain  yield  zone  (i.e.,  minimal  yielding 
ahead  of  the  crack  tip  and  a  very  upright  yield  zone) .  The  maximum  extent 
of  yielding  is  30. 7Z  of  the  half  crack  length  which  is  consistent  with  the 
small  strain  assumptions  made  in  the  analysis  requiring  contained  yielding. 
Figure  3b  is  a  plot  of  the  surface  zones  for  a  specimen  with  total  thickness 
of  1.27  cm.  The  yield  zone  is  slightly  wider  (more  rounded)  with  this 
thickness.  The  maximum  radius  is  now  32.7%  of  the  half  crack  length  and  the 
yielding  ahead  of  the  tip  has  increased  (though  it  is  still  small) .  The 
zone  still  maintains  the  basic  plane  strain  characteristics  at  this  thick- 


Figure  3c  shows  the  surface  yield  ■ 

zone  for  a  specimen  with  thickness  of  | 

6.35  mm.  The  zone  is  now  much  wider  '  / 

with  a  larger  maximum  radius  and  /  /  \ 

yield  extent  ahead  of  the  tip.  The  /  \  ^  ) 

zone  no  longer  exhibits  the  plane  /  ^  / 

strain  characteristics  but  is  in  j  / 

transition  between  plane  strain  and  I  / 

plane  stress.  Figure  3d  is  a  plot  of  .  \  /  / 

the  surface  yield  zone  for  a  specimen  \  /  j/ 

with  total  thickness  of  3.175  mm.  _ _  \  / 

The  zone  is  significantly  more  * 

rounded  than  any  of  the  previous  r/»  «  .ois 

zones  with  a  larger  maximim  radius  ^  ^  , 

and  yield  extent.  The  maximum  yield  yield  zones  for 

radii  and  extent  of  yielding  ahead  specimen  with  2T  -  2.54  cm. 

of  the  crack  tip  for  the  four  thick¬ 
ness  surface  zones  are  given  in  '  ^ 

Table  la.  These  yield  parameters  /N 

both  increase  with  decreasing  thick-  ^  /  \ 

ness  as  was  expected.  The  final  ,  /  y  \ 

zone  at  a  thickness  of  3.175  mm  has  /  •$'.  I 

the  roui^.ded  characteristic  of  a  /  ■>  '  J 

plane  stress  yield  zone.  The  direc-  j  y 

tion  of  maximum  yielding,  however,  I  ^  jT 

is  still  a  fairly  large  angle  rela-  \  ^ 

tive  to  the  crack  line  suggesting  \ ^ ^ - 

some  influence  of  dilatation.  I  1 

Though  for  this  problem  (with  a 
relatively  small  amount  of  plastic  „ 

deformation  present)  the  difference  Surface  yield  zones  for 

between  the  maximum  radii  is  not  specimen  with  2T  -  1.27 

large,  the  nature  and  extent  of 

yielding  ahead  of  the  crack  tip  show  a  large  dependence  on  the  specimeu 
thickness. 


Fig.  3b. 


Surface  yield  zones  for 
specimen  with  2T  ■  1.27  cm. 


Figure  4a  is  a  plot  of  the  von  Mises  stress  contour  corresponding  to  the 
material  yield  stress  on  the  midplane  of  the  2.54  com  thick  specimen.  The 
zona  is  typical  of  plana  strain  zones  and  is  smaller  than  the  surface  zone 
for  the  same  thickness  specimen.  The  shape  of  the  zone  with  a  minimal  ex¬ 
tent  of  yielding  ahead  of  the  crack  tip  suggests  high  dilatation  in  that 
region.  The  midplane  zone  for  the  1.27  cm  thick  specimen  is  shown  in 


Fig.  4b.  The  zone  is  larger  than  that 
of  the  thicker  specimen,  however, 
there  is  still  minimal  yielding  ahead 
of  the  tip.  The  angle  of  maximum 
yielding  is  more  acute  than  in  the 
thicker  specimen.  The  stress  state, 
however,  would  still  be  characterized 
by  plane  strain. 

Figure  4c  shows  the  midplane  yield 
zone  for  the  6.35  mm  thick  specimen. 
The  zone  is  considerably  wider  and 
more  rounded  than  for  the  thicker 
specimens.  It  shows  characteristics 
of  both  plane  strain  and  plane 
stress  zones  suggesting  a  region  of 
transition.  Figure  4d  is  a  plot  of 
the  midplane  yield  zone  for  the 
3.175  mm  thick  specimen.  The  zone 
is  basically  a  plane  stress  zone 
and  is  larger  than  for  the  thicker 
specimens.  The  maximum  yield  radii 
and  radius  of  yielding  ahead  of 
the  crack  tip  on  the  specimen  mid~ 
planes  are  given  in  Table  lb.  Both 
increase  with  decreasing  thickness  ' 
as  was  expected.  In  all  cases,  the 
midplane  yield  zones  are  smaller 
than  the  surface  zones. 


Fig.  4a.  Midplane  yield  zones  for 

specimen  with  2T  ••  2.54  cm. 


Fig.  4c.  Midplane  yield  zones  for 

specimen  with  2T  ■  6.35  mm. 


Fig.  3c.  Surface  yield  zones  for 

specimen  with  2T  ■  6.35  mm. 


Fig.  3d.  Surface  yield  zones  for 

specimen  with  2T  ■  3.1"/b  mm 


r/a  •  0.071 

Fig.  4b.  Midplane  yield  zones  for 

specimen  with  2T  ■  1,27  cm. 


Fig.  Ad.  Midplaue  yield  zones  for 

specimen  with  2T  ■  3.175  mm 


TABLE  la  Yield  Radii  as  a  function 
of  Thickness  for  Surface 
Yield  Zones. 


TABLE  lb  Yield  Radii  as  a  Function  of 
Thickness  for  Midplane  Yield 
 Zones.  


Thickness  'nsk/s 


T  -  1.54  ea  0.307  0.045 
T  -  1.27  ea  0.327  0.075 
T  -  4.35  M  0.331  0.205 
I  •  3.175  B  0.343  0.2M 


r  maximum  yield  radius 

max  _ 


IMckMsa 

*o/. 

T  •  2.54  a 

0.260 

0.039 

T  -  3.27  a 

0.280 

0.071 

T  -  6.35  ■■ 

0.283 

0.193 

T  ■  3.375  tm 

0.299 

0.268 

r  -*•  yield  radius  along  crack  line 

o  _ _ 


CONCLUSIONS 


The  results  of  this  study  demonstrate  the  thickness  dependence  of  the  yield 
zones  near  a  crack  front  on  specimen  thickness.  It  is  shovm  that  both  the 
extent  of  plastic  deformation  and  the  dominance  of  deformation  type  (i.e., 
dilatation  or  distortion)  are  controlled  by  the  thickness.  The  nature  of 
the  deformation  is  fundamental  to  the  understanding  of  the  incipient  frac¬ 
ture  processes.  The  delineation  of  the  fundamental  deformation  response 
near  a  three-dimensional  crack  front  is  an  imperative  first  step  in  the 
understanding  and  accurate  prediction  of  ductile  fracture  processes. 


To  further  the  understanding  of  ductile  fracture,  it  is  necessary  to  compare 
theoretical  and  experimental  deformation  predictions  local  to  the  crack 
front.  Only  thx’ough  such  comparisons  can  an  assessment  be  made  of  the 
accuracy  and  rft.llability  of  the  numerical  methods  fr  -  analysis.  To¬ 

ward  this  goal,  it  is  proposed  to  measure  the  tesinuJ.i,  uef tu  ,'iation  on  the 
surface  of  the  specimen  in  the  unloaded  state.  The  theoretical  study  pre¬ 
sented  above  demonstrates  that  the  finite  element  prediction.'!  are  qualita¬ 
tively  realistic  and  sensitive  to  specimen  thlcknes.s.  Comparison  with  ex¬ 
perimental  results  will  delineate  the  grid  characteristics  and  hardening 
models  which  best  model  specific  geometric  and  material  applications.  After 
successful  "tuning"  of  the  finite  element  model,  a  complete  description  of 
the  stress  and  energy  state  in  a  cracked  body  can  be  predicted  with  confi¬ 
dence.  Once  fully  three-dimensional  stress  fields  are  predicted,  ductile 
failure  theories  can  be  tested  and  skeptically  compared  without  the  bias  of 
unrealistic  analytical  approximations. 
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